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1.  Executive  Summary 

1.1  Survey  of  effort 

This  document  summarizes  effort  conducted  under  the  contract  F49620-92-C-0006, 
“Development  of  Combined  Asymptotic  and  Numerical  Procedures  for  Trarisonic  and 
Hypersonic  Rows.”  In  the  hypersonic  portion  of  the  effort,  stability  and  transition, 
unsteady  aerodynamic  characteristics  and  stability,  strong  inviscid/viscous  interactions 
and  potential  q^plication  of  instabilities  to  scramjet  combustors  were  emphasized.  Major 
aspects  of  the  work  included  a  study  of  three-dimensional  disturbances  in  h^rsonic 
boundary  layers.  These  are  particul^y  important  to  scramjet  inlet  design.  We  have 
investigated  low  frequency  instabilities  of  viscous  hypersonic  flow  over  a  flat  plate  using 
triple-deck  methods.  Our  eigensolutions  indicate  the  existence  of  an  upstream 
propagating  stable  mode  which  is  the  unsteady  counterpart  of  the  steady  two-dimensional 
Neiland-Brown-Stewartson  solutions.  We  have  also  discovered  a  downstream- 
propagating  unstable  mode.  This  leads  to  a  three-dimensional  instability  and  is  controlled 
moderately  supersonic  flow  physics.  Calculations  of  growth  rates  and  analysis  of  the 
forced  problem  for  the  experimentally  significant  time-harmonic  vibrator  show  that  at 
sufficiently  low  frequencies,  the  downstream  far  field  transverse  projection  consists  of 
Gaussian  lobes  symmetrically  displaced  about  the  source  axis.  For  frequencies  exceeding 
a  critical  value,  only  one  lobe  occurs.  The  width  of  this  lobe  contracts,  in  contrast  to  the 
expansion  of  conventional  wakes.  For  wavelengths  of  the  order  of  the  boundary  layer 
thickness  and  phase  speeds  less  than  unity,  we  have  discovered  new  waveguide 
reflectional  inst^ility  mechanisms  that  can  be  important  for  strong,  moderate  and  weak 
interaction  between  the  shock  and  boundary  layers.  These  results  can  be  readily  extended 
to  three-dimensional  flows.  Our  asymptotic  analyses  show  that  to  first  order  in  the 
displacement  thickness,  the  upper  boundary  layer  edge  behaves  as  a  solid  wall  that 
reflects  boundary-layer  disturbances.  This  is  a  viscous  generalization  of  energy  trapping 
instability  mech^isms  discussed  by  the  Principal  Investigator  in  an  inviscid  context  in 
[Cl].  These  compete  with  other  more  well  known  instability  modes  that  produce 
premature  transition  and  are  therefore  of  importance  to  hypersonic  design  and  propulsion- 
airframe  integration.  In  addition  to  this  boundary  layer  “lower”  wave  guide,  and  “upper” 
one  exists  between  the  shock  and  the  sonic  line.  At  mode  synchronism  points  between 
upper  and  lower  wave  guides,  singular  growth  spikes  occur  that  invalidate  the  modal 
dMomposition  used  in  state-of-the-art  linear  stability/transition  prediction  theories/codes. 
The  waveguide  concepts  have  been  studied  in  connection  with  control  of  transition 
by  absorbing  walls  and  demonstrated  rcdoction  of  amplification  with  increased  wall 
absorption.  From  the  stability  asymptotics,  practical  transition  delay  concepts  have 
been  identified  using  conventional  thermal  protection  materials  having 
ultrasonlcally  selective  properties.  For  thin  shock  layers,  we  have  extended  steady 
Newtonian  asymptotic  theory  to  the  unsteady  case.  We  have  used  the  theory  to 
demonstrate  the  existence  of  a  new  finite-ampUtndc  Inflectional  Instability  that  may 
be  ttseftil  to  promote  pre>mixtng  and  enhanced  bnming  efficiency  for  scramjet 
combustors.  The  unsteady  Newtonian  model  is  useful  in  other  contexts  such  as 
hypersonic  oeroelasticity  and  flutter.  As  another  application  of  our  hypersonic 


M 


SCTlOeiFTR 


theories,  similitudes  have  been  used  to  collapse  re-entry  cone  transition  data  to 
provide  an  efficient  prediction  tool  for  the  sdentist,  engineer  and  designer.  In  the 
transonic  regime,  we  have  developed  the  basis  for  wind  tunnel  wall  interference 
predictions  for  slender  and  not-so-slendcr  bodies  in  wind  tunnels.  The  singular 
perturbation  theory  is  interacting  with  experiments  that  we  are  conducting  in  the  Russian 
TsAGI  T-128  wind  tunnel  le^ng  to  the  generalization  of  our  Transonic  Area  Rule  for 
Wind  Tunnel  Wall  Interference  (TARWI)  to  test  articles  of  length  comparable  to  the  wall 
hci^L  The  TsAGI  tests  have  validated  our  matched  asymptotics  derivation  of  the 
TARWL  The  TARWI  reduces  the  three-dimensional  wall  Interference  correction 
problem  to  mic  in  two-dimensions.  This  will  provide  test  engineers  a  basis  for 
estimation  of  transonic  blockage  corrections  and  optimization  of  model  size  for 
maximum  unit  Reynolds  number  and  minimum  wall  interference.  Also  to  be  described  in 
what  follows  are  other  theories  for  lift  interference  in  porous  wind  tunnels  as  well  as 
corrections  for  non-circular  walls  that  we  have  derived  in  the  effort.  To  our  knowledge, 
this  is  the  first  time  .such  corrections  have  been  embedded  in  systematic  asymptotic 
approximation  procedures. 
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•  Dr,  Norman  D.  Malmuih 

•  Professor  Julian  D.  Cole 

•  Dr.  Alexander  V.  Fedorov 

•  Dr.  Andii  Khokhlov 

•  Professor  Vladimir  Ya  Neiland 
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•  I  lofessor  Oleg  Ryzhov 
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Asbpis,  T.B.  Oatski  and  R.  Hirsh.  Kluwer  Academic  Publishers,  1993. 

C2.Malmuth,  N.D.,  Jafroudi,  H.,  Wu,  .C.C,  McLachlan,  R.  and  Cole,  J.D.  “Asymptotic 
methods  applied  to  transonic  wall  interference,"  AIAA  J.  31  5,  pp.  91 1-918,  May 
1993. 
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Flows,  Hotel  Amber  Baltic,  Miedzyzdroje,  Poland,  July  12-14,  1993  in  Advances  in 
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2.  Mechanism  of  Hypersonic  Flow  Stabilization  by  a  Semi- 
Transparent  Wall 


2.1  Introduction 

The  ability  to  predict  and  control  the  location  of  boundary  layer  transition  is  of  critical 
importance  in  the  design  and  optimiratior.  of  aerospace  vehicles  [1).  Transition  causes 
significant  increases  in  heat  transfer  and  skin  friction.  This  leads  to  diminished  vehicle 
performance,  primarily  because  of  the  additional  wei^t  of  thermal  protection  system 
(TPS). 

If  external  disturbances  ore  small  and  the  TPS  roughness  is  negligible,  then  the  transition 
process  is  due  to  amplification  of  unstable  modes  in  the  boundary  layer  [I],  [2].  In  this 
case,  stability  theory  and  experiment  form  a  foundation  for  the  prediction  and  control  of 
transition  [3].  [4],  [3].  Several  difierent  instability  mechanisms  may  be  responsible  for 
transition  on  the  hypersonic  vehicle  surface,  nantely 

1.  the  first  and  second  modes  which  are  dominant  in  2*D  or  quasi  2>D  boundary  layers  on 
a  flat  plate,  axisymmetric  bodies,  wings,  etc. 

2.  Gortler  vortices  which  can  play  a  major  role  on  concave  surfaces. 

3.  cross  flow  instability  which  can  be  dominant  in  3-D  boundary  layers  on  the  leading 
edge  of  a  swept  wing,  axisymmetric  bodies  at  the  angle  of  attack,  etc. 

Our  analysis  addresses  the  mechanism  in  Item  (I).  Its  relation  to  other  types  of  instability 
will  be  discussed  at  the  end  of  this  chapter. 

The  first  mode  is  an  extension  to  high  speeds  of  the  Tollmien-Schlichting  waves.  They 
represent  viscous  instability  at  low  Mach  numbers.  However  the  inviscid  nature  of  the 
first  mode  begins  to  dominate  when  the  Mach  number  increases,  since  the  compressible 
boundary  layer  profiles  contain  a  generalized  inflection  point  [3].  This  mode  may  be 
stabilized  by  wall  cooling,  suction  and  favorable  pressure  gradient 

The  second  mode  is  the  result  of  an  inviscid  instability  present  due  to  a  region  of 
supersonic  mean  flow  relative  to  the  disturbance  phase  velocity.  It  belongs  to  the  family 
of  trapped  acoustic  modes.  For  an  insulated  surface,  these  higher  modes  appear  for  Mach 
numbers  M  >22.  However,  it  is  not  until  the  Mach  number  is  of  the  order  of  4  or 
greater  that  the  second  mode  is  at  low  enough  frequency  to  have  experimental 
consequences.  Once  the  .second  mode  sets  in,  k  becomes  the  dominant  instability  since  ks 
growth  rate  tends  to  exceed  that  of  the  first  mode.  For  insulated  surfaces,  this  occurs  for 
M  >4.  For  cooled  surfaces,  the  second  mode  can  dominate  at  even  lower  Mach 
numbers.  In  contrast  to  the  first  mode,  the  higher  subsonic  modes  are  destabilized  by 
cooling.  The  existence  of  the  second  modes  was  established  by  the  experiments  of 
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Kendall  [6],  Demetriades  [7],  and  Stetson  et  al.  (8I-[10].  The  most  amplified  second 
mode  wavelengths  tend  to  be  approximately  twice  the  boundary  layer  thickness,  and  the 
second  mode  phase  velocities  tend  to  be  approximately  equal  to  the  edge  velocity  of  mean 
flow. 

The  second  mode  induces  pressure  disturbances  of  very  high  frequency.  For  example,  the 
most  amplified  waves  observed  in  the  experiment  of  Stetson  et  al.  [8]  at  Mach  8 
correspond  to  a  frequency  of  about  100  kHz.  TPS  porous  surfaces  can  absorb  these  high 
frequency  fluctuations.  It  is  assumed  that  the  absorption  can  cause  stabilization  of  the 
second  and  higher  modes  by  this  energy  extraction  mechanism.  This  assumption  has 
been  examined  in  a  theoretic^  study  that  is  summarized  in  what  follows.  The  eigenvalue 
problem  for  inviscid  disturbances  in  a  supersonic  boundary  layer  is  formulated  in  Section 
2.2.*  The  WKB  method  is  used  to  specify  the  boundary  condition  on  the  semi-transparent 
wall  for  the  acoustic  normal  modes  of  high  frequency  including  the  second  mode.  In 
Section  2.3,  the  eigenvalue  problem  is  solved  numerically  and  ft  is  shown  that  absorption 
stabilizes  the  second  mode.  A  discussion  of  the  future  effort  concludes  this  chapter. 

2Jt  Eigenvalue  problem  forlnvlacld  disturbances 

A  supersonic  boundary-layer  flov/  over  a  flat  plate  or  sharp  cone  schematically  .shown  in 
Fig.  1  is  considered.  The  coordinates  X ,  y  ire  made  nondimcnsionless  using  the 

boundary  layer  displacement  thickness  S* .  In  the  locally  parallel  approximation,  the 
mean  flow  is  characterized  by  the  profile  of  -component  velocity  U(y)  and 
temperature  T(y)  nondimensionaiized  using  their  corresponding  reference  quantities  U, 

and  at  the  upper  boundary  layer  edge.  The  inviscid  stability  equations  are  obtained 
from  the  linearized,  parallel  flow,  viscous  equations  for  a  perfect  gas  with  constant 
specific  heat  by  setting  all  viscous  and  heat-conduction  terms  equal  to  zero  [II].  Since 
the  second,  third  and  higher  modes  have  maximum  growth  rate  in  the  two-dimen.siona] 
case^  we  consider  2-D  disturbances  in  normal-mode  form  as 

[u,v,p,eY(.x,y,t)=[u,v,p,eY{y)exp[i((Xx-QX)]  (1) 

where  S  and  P  are  the  velocity  components  in  the  X  and  y  directions,  P  is  the 
pressure  normalized  by  the  d)mamic  pressure,  pJUg  ,  9  is  the  temperature,  Ct  is  the 
wavenumber,  and  (0  is  frequency. 


‘This  work  and  a  portion  of  that  described  in  the  other  chapters  was  partially  funded  by  other 
sources  in  addition  to  the  contract. 

'  This  is  also  termed  the  most  dangerous  case. 
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Fig.  1  Acoustic  mode  in  a  supersonic  boundary  layer  on  semi-transparent  wall. 


The  system  of  equations  for  v(y)  and  p(y)  ( see  for  example  [11])  is 


v'  = - v  +  ia - - — 

U-c  ,  U-c 


(2) 


where  c  =  <u  /  a  is  the  phase  speed,  and  prime  denotes  the  derivative  with  respect  to  y . 
The  boundary  conditions  for  Eqs.  (2)  arc 


v(0)  =  Ap(0) 

(3) 

o 

II 

T 

(4) 

Equation  (3)  is  a  general  form  (Darcy’s  Law)  of  the  inviscid  boundary  condition  on 
porous  or  compliant  walls.  The  coefficient  A  is  a  complex  quantity  which  depends  on 
properties  of  the  wall  material,  mean  flow  characteristics  on  the  wall  surface,  and  flow 
perturbation  parameters  such  as  a  wave  frequency  and  wavelength.  Note  that  the 
boundary  condition  (4)  requires  the  normal  mode  to  vanish  outside  the  boundary  layer. 
The  eigenvalue  problem  (2)-(4)  provides  the  dispersion  relation  F(a,0))  =  0.  For 
temporal  stability,  the  wavenumber  Cl  is  real,  and  (i){a)  is  the  complex  eigenvalue.  If 
Imfm)  >  0,  then  the  disturbance  is  unstable. 

The  system  of  equations  (2)  can  be  reduced  to  one  equation  for  the  pressure  amplitude 
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T 


with  the  following  boundary  conditions 

AO)  =  A~piO)  (6) 

p(oo)=0  (7) 

According  to  the  asymptotic  analysis  [12]  characteristics  of  the  second,  third  and  higher 
modes  can  be  obtained  by  the  WKB  method.  Without  going  into  details,  we  briefly 
describe  the  results  of  this  analysis.  For  a  subsonic  disturbance  with  phase  speed 

1-— <Re(c)<l+—  (8) 

M  M 

the  boundary  layer  behaves  as  an  acoustic  waveguide  schematically  shown  in  Fig.  I.  In 
the  waveguide  region  between  the  wall  and  the  sonic  line; 
>'  =  )  =  ReW  “  ^(y^ )  *  pressure  amplitude  piy)  oscillates  in  they  direction. 

Here  the  real  part  of  the  function  A(y)  is  large  and  positive,  and  hence  the  solution  of  Eq. 
(5)  can  be  expressed  in  the  following  WKB  form 


p(y)  ~  My) «xp  -t J  ^\+ My)  * J  '^y  0{e) 


n  U-c  .  MHU-cf  , 

where  £  =  1  /  max(|A|}  is  a  small  parameter.  The  Hrst  and  second  terms  respectively 
represent  incident  and  reflected  acoustic  waves  to  the  wall.  Near  the  sonic  line,  the 
acoustic  ray  turns  around  its  schematically  shown  in  Fig.  1.  In  the  upper  region  y>y,  , 

the  real  part  of  A(>')  is  negative  and  the  pressure  amplitude  decays  exponentially.  For  the 


neutral  modes,  the  dispersion  relation  obtained  by  the  WKB  method  has  the  following 
form  for  the  case  of  a  non-transparent  wall  ( A  s  0  ) 
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The  solution  of  Eq.  (  10),Q)^{(X),  corresponds  to  the  second  mode  for  «  =  0,  to  the 
third  mode  at  n  =  1  etc. 

To  obtain  the  explicit  form  of  the  coefficient  A  on  a  semi-transparent  wall,  we  use  the 
WKB  approximation  (9).  The  reflection  coefficient  T  is  defined  as 


P2i0) 

pAO)’ 


(11) 


If  T  =  1 ,  then  disnifbances  are  completely  reflected  by  the  wall.  If  T  <  1 ,  then  some 
portion  of  the  disturbance  energy  is  absorbed  by  the  wall  material.  Differentiating  the 
expression  (9)  we  obtain  in  the  first  order  approximation 


p'(0)  =  -(A(0)[A(0)-A(0)) 


Substituting  (12)  into  the  boundary  condition  (6)  and  using  Eqs.  (9),  (1 1),  we  get 
p'(0)  =  -/A(0)p(0)~  ( 

From  Eqs.  (6),  (13)  we  obtain  the  explicit  form  for  the  absorption  coefficient 


c  V  ^(0)  ’  1+f 


2.3  ReaultM  and  df9cut9lon 

The  eigenvalue  problem  (2H4)  is  solved  by  numerical  integration  using  the  method  of 
Mack  [11].  The  absorption  coefficient  A  is  specified  by  the  expression  (14).  The 
integration  starts  at  the  upper  freestream  and  continues  to  the  wall  along  an  indented 
contour  in  the  complex  plane.  Tlic  contour  lies  below  the  critical  point  >'c:f/{y,)  =  c. 
The  eigenvalue  a)(a)  is  found  by  the  iterative  Newton  method. 
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&)  Fig.  2,  the  second-mode  growth  rate  Ini(a>)  is  .shown  as  a  function  of  the  wavenumber 
CC,  at  various  values  of  the  reflection  coefficient  T  .  Calculations  were  performed  for  the 
flat-plate  boundary  layer  on  a  cool  wall  with  a  temperature  ratio  typical  of  a  hypersonic 
vehicle  surface.  Here,  7.  is  the  surface  temperature,  and  7^  is  the  temperature 
ratio  r,  /  7^  «>  0.2  on  an  adiabatic  wall.  The  freestream  parameters  have  the  following 
values:  Mach  number  M  =  6,  sp«:ific  heat  ratio  y  =  1.4 ,  and  Prandtl  number  Pr  =  0.72 . 
It  is  seen  that  the  semi-transparent  wall  reduces  the  growth  rate  by  the  factor  of  4  at  the 
reflection  coefficient  r  =  0.6.  A  similar  stabilization  effect  has  been  obtained  for  other 
Mach  numbers,  and  wall  temperatures.  Since  the  physics  of  the  third  and  higher  modes  is 
similar  to  that  of  the  second  mode,  it  is  assumed  that  the  absorption  of  acoustic  energy 
can  cause  similar  stabilization  of  these  higher  modes. 


Fig.  2  Effect  of  the  wall  reflection  coefficient  T  on  the  second  mode  growth  rate 
Iin((30)  for  hypersonic  boundary  layer  at  Af  =  6,  7,  /  7^  =  0.1986 ,  Pr  =  0.72, 
and  7  =  L4 . 

Z4  ConetualonM 

1.  According  to  inviscid  stability  theory,  the  absorption  of  acoustic  energy  by  a  semi¬ 
transparent  wall  can  essentially  stabilize  the  second  mode  of  a  hypersonic  boundary 
layer.  This  stabilization  may  Ic^  to  delay  of  transition  on  hypersonic  vehicle  surfaces. 

2.  Since  frequency  of  pressure  disturbances  incident  on  the  wall  by  the  second  mode  are 
very  high  (of  the  order  of  100  kHz),  TPS  porous  materials  used  on  hypersonic  vehicle 
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surfaces  can  be  semi-transparent  with  respect  to  these  disturbances.  It  is  this  property 
that  provides  an  opportunity  to  use  absorption  to  maintain  or  increase  the  laminar  flow 
mn. 


Since  the  boundary  layer  stabilization  is  due  to  the  absorption  of  disturbance  energy  by 
the  wall,  we  believe  that  a  similar  effect  may  occur  for  other  high-frequency  instabilities 
such  as  unsteady  cross  flow  vortices.  A  semi-transparent  wall  may  also  affect  transition 
past  TPS  roughness  elements.  &  is  recommended  that  these  assumptions  be  verified  by 
wind  tunnel  experiments  using  models  with  porous  surfaces. 

In  our  analysis,  we  have  used  a  simple  model  for  the  boundary  condition  on  a  semi¬ 
transparent  wall.  This  model  permits  rough  estimates  of  the  absorption  effect  and 
formulations  for  its  experimental  verification.  For  reliable  quantitative  results,  more 
comprehensive  models  of  the  absorption  mechanism  as  well  as  detailed  experimental 
.studie.s  are  required.  In  order  to  demonstrate  the  effect  of  transition  delay  on  a  semi¬ 
transparent  surface  and  apply  it  to  hyjiersonic  vehicle  design,  the  following  effort  should 
be  performed: 

•  Test  current  TPS  materials  and  obtain  their  transparency  characteristics  in  a  frequency 
band  relevant  to  hypersonic  instability 

•  Demonstrate  the  effect  of  surface  transparency  on  transition  by  experiment  on  models 
of  a  simple  shape  (such  as  a  flat  plate  or  sharp  cone)  in  hypersonic  wind  tunnels 

•  If  the  demonstration  is  successful: 

-  Develop  theoretical  models  accounting  for  viscous  effects  as  well  as  physics  of 
interaction  between  the  flow  disturbance  and  the  TPS  material  microstructurc. 

-  Conduct  detailed  experiments  in  hypersonic  wind  tunnels. 

-  Perform  flight  testing. 

-  Develop  new  TPS  materials  providing  maximum  stabilization  of  hypersonic 
boundary  layer  and  transition  delay. 
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3.  Reflectional  Instability  of  a  Hypersonic  Boundary  Layer 

3.1  Introduction 

The  hypersonic  speed  range  contains  some  of  the  most  important  problems  in  boundary 
layer  transition.  At  high  Mach  numbers,  early  transition  causes  temperatures  high  enough 
to  destroy  air  vehicle  surfaces.  Uncertainty  in  the  transition  locus  leads  to  diminished 
performance,  primarily  because  of  the  additional  weight  of  the  thermal  protection  system. 
This  has  stimulated  considerable  experimental  and  theoretical  studies.  The  transition 
problem  is  an  initial  boundary  value  problem  which  requires  the  identification  of  the 
freestream  disturbance  field,  understanding  of  the  laminar  flow  receptivity  to  external 
disturbances  including  surface  irregularities,  as  well  as  linear  and  nonlinear  amplification 
of  unstable  disturbances  [1].  Current  knowledge  is  inadequate  to  solve  and  apply  the 
initial  value  problem  to  design,  even  at  low  speeds.  However,  it  has  been  found  that 
transition  can  be  correlated  with  linear  stability  theory  if  the  external-disturbances  are 
relatively  small.  This  motivates  development  of  hypersonic  flow  stability  theory. 

Hypersonic  instability  is  more  complicated  than  for  moderate  speeds.  In  addition  to 
Tollmien-Schlichting  waves,  cross  flow  and  Gbrtler  vortices,  the  following  special 
observations  apply  at  high  Mach  numbers: 

•  In  the  region  of  the  boundary  layer  that  is  supersonic  relative  to  the  phase  velocity, 
Mack  [2]  showed  that  multiple  solutions  of  the  inviscid  stability  equations  occur. 
These  additional  solutions  are  called  higher  modes.  Physically  they  represent  acoustic 
modes  which  propagate  in  a  ‘‘waveguide”  located  between  the  wall  and  the  sonic  line. 
The  first  mode  of  this  family,  which  is  called  the  second  mode,  is  the  most  unstable. 
Its  wavelength  is  of  the  order  of  boundary  layer  thickness. 

•  Fedorov  and  Khokhlov  [3]  demonstrated  that  the  second  mode  instability  is  related  to 
synchroniaation  between  the  first  and  second  modes.  Near  the  .synchronism  point, 
where  frequencies  and  phase  speeds  of  both  modes  coincide,  the  eigenvalues  split  into 
two  branches.  In  the  vicinity  of  branch  points  the  dispersion  relation  is  singular, 
causing  a  strong  exchange  between  modes  due  to  nonparallclli.<m.  This  implies  that  a 
multi-mode  approach  is  required  for  transition  prediction  rather  than  the  conventional 
single  mode  method. 

•  If  the  shock  wave  is  close  to  the  surface,  which  is  relevant  to  the  windward  side  of  re¬ 
entry  vehicle  at  the  high  angle  of  attack,  a  new  class  of  unstable  modes  occurs  due  to 
reflection  of  disturbances  by  the  shock.  Qualitative  asymptotic  analy.ses  of  Fedorov 
and  Gushchin  [4],  [5]  showed  that,  in  addition  to  the  lower  waveguide  of  the  boundary 
layer,  there  is  an  upper  waveguide  located  between  the  shock  wave  and  the  sonic  line. 
At  certain  conditions  acoustic  waves  of  the  lower  waveguide  are  in  resonance  with 
those  of  the  upper  waveguide  and  have  identical  phase  speeds.  The  synchronization 
between  modes  causes  a  new  type  of  instability. 
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•  In  the  hypersonic  limit,  the  inflectional  instability  is  related  to  a  thin  transitional  layer 
located  near  the  upper  boundary  layer  edge.  This  inviscid  instability  is  called  the 
vorticity  mode.  It  has  been  studied  by  Blackaby,  Cowley  and  Hall  [6],  Smith  and 
Brown  [7],  Grubin  and  Trigub  [S]  using  asymptotic  methods.  They  analyzed  stability 
of  the  hypersonic  flow  over  a  flat  plate  for  weak  and  strong  viscous-inviscid 
interaction.  It  has  been  found  that  the  vorticity  mode  weakly  disturbs  the  flow  outside 
the  transitional  layer  which  separates  a  hot  flow  of  the  boundary  layer  from  a  cool 
inviscid  flow.  The  phase  speed  of  the  vorticity  mode  tends  to  unity  With  increasing 
Mach  number.  Its  wavelength  is  much  less  than  the  boundary  layerthickness. 

•  Malmuth  [9]  studied  dishirbance  flowfleld  in  a  thin  shock  layer  between  the  shock 
wave  and  the  upper  boundary  layer  edge.  He  examined  the  strong  interaction  regime 
on  a  flat  plate  at  the  Newtonian  limit  when  the  specific  heat  ratio  y  ->  1 .  His  studies 
show  that  disturbances  of  longitudinal  scale  of  the  order  of  the  length  of  the  strong 
interaction  region  can  strongly  amplify  in  a  thin  shock  layer.  This  new  instability  does 
not  relate  to  the  normal-mode  decomposition  normally  used  in  stability  theory.  It  was 
revealed  only  by  solving  the  complete  initial  boundary  value  problem  for  strongly 
nonparallel  flow. 

•  Hypersonic  flow  past  a  blunt  body  generates  a  strong  curved  shock  near  the  nose.  This 
causes  an  enttt^y  layer  in  the  inviscid  region  between  the  shock  and  the  boundary 
layer.  As  the  shock  wave  approaches  the  limit  of  a  sharp  cone  solution,  the  entropy 
layer  thickness  decreases  and  eventually  it  is  swallowed  by  the  boundary  layer.  The 
entropy  layer  can  be  unstable  due  to  presence  of  the  inflection  point  [5],  [10]. 
Reshotko  and  Khan  [10]  found  that  within  the  swallowing  region,  the  boundary  layer 
stability  is  affected  by  the  nonuniform  flow  near  the  boundary  layer  edge.  The 
interplay  between  the  entropy  modes  and  the  boundary  modes  may  strongly  affect  their 
amplification. 

Coexistence  of  several  instabilities  of  essentially  different  scale  make  it  difficult  to 
predict  transition  from  a  universal  stability  code.  Computations  of  the  hypersonic  mean 
flow  and  its  stability  from  codes  develop^  for  moderate  Mach  number  flows  are  very 
time-consuming.  Computational  time  and  intensity  grows  dramatically  with  increasing 
Mach  number  because  of  the  strong  non-uniformity  of  hypersonic  mean  flow  profiles  and 
the  high  sensitivity  of  the  related  eigenvalue  problem  to  its  first  and  second  derivatives 
[1].  Solution  proc^ures  based  on  combined  asymptotics  and  numerics  deal  effectively 
with  this  difficulty  since  they  systematically  guide  accurate  resolution  of  the  high 
gradient  regions  in  a  robust  way.  This  has  stimulated  the  research  described  here. 

Asymptotic  models  have  been  developed  to  treat  the  short-scale  vorticity  mode  of  the 
transitional  layer  in  [6]  •  [8],  and  large-scale  disturbances  in  [9]  of  the  inviscid  shock 
layer.  Both  instabilities  convect  with  a  phase  speed  close  to  the  freestream  velocity.  The 
intermediate  case  in  which  the  distuibance  scale  is  of  the  order  of  boundary  layer 
thickness  and  its  phase  speed  is  less  than  the  freestream  speed  has  not  been  studied  by 
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asymptotic  methods.  This  gap  is  targeted  in  the  analysis  that  follows.  In  Section  3.2,  the 
composite  system  of  stability  equations  for  a  strong  interaction  hypersonic  boundary  layer 
is  derived.  The  disturbance  field  in  the  inviscid  shock  layer  is  treated  in  Section  3.3.  In 
Section  3.4,  matching  between  the  boundary  layer  and  shock  layer  is  used  to  establish 
boundary  conditions  at  the  upper  boundary  layer  edge.  The  weak  interaction  case  is 
considered  in  Section  3.5.  Section  3.6  generalizes  consideration  of  the  eigenvalue 
problem  to  3-D  mean  flows  and  fluctuations.  In  Section  3.7,  a  self-similar  solution  for 
steady  strong  interaction  boundary-layer  flow  is  obtained.  These  mean  flow  profiles  are 
used  for  stability  computations  in  Section.3.8. 


3J  Stability  equathna  In  th»  strong  Interaction  boundary  layer 

Consider  the  boundary  layer  on  a  flat  plate  in  hypersonic  flow  of  a  perfect  gas  (see  Fig.  1) 
with  speed  f/ * ,  temperature  72 .  and  density  p2 .  Introduce  nondimensional  coordinates 
and  time  as 

(x.y)  =  (x,y)lS\t-.n(S‘lu:).  (1) 


where  8'  ii  of  the  order  of  the  boundary  layer  thickness.  The  pressure  p\  temperature, 
density  and  the  velocity  components  u*.v*  are  nondimensionalized  as 


v: 


(2) 


Besides  the  fast  variable  x ,  a  slow  variable  x^  corresponding  to  the  longitudinal  scale  V 
(see  Fig.  I)  is  introduced  ,  where 


x^=8•x,S^^,  (3) 

Consider  the  following  hypersonic  limit  corresponding  to  the  strong  interaction  between 
inviscid  flow  and  boundary  layer  (see  for  example  Bush  [1 1]) 

^-♦0,Af-+-,5’Af’->-.(7'-l)  =  0(l).  (4) 

where  /  is  the  specific  heat  ratio.  In  this  case,  the  shock-layer  thickness  8,  is  of  the 
order  of  the  boundary  layer  thickness  8' . 
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Fig.  1  Hypersonic  smjng  interaction  flow  over  a  flat  plate. 

The  mean  flow  and  its  unsteady  peiturbations  are  represented  asymptotically  as 

p(x,y./)  =  5*[F,(.r,)+ep(x,x„y,/)]+-- 

p(x,y./)  =  ^~5^[;?,(x„y)+ep(x.x,,y,r)]+- 

(x,  ,  y)  +  « •  7*(x,x, ,  y,  r)] + •  •  • 
y.  r)  =  (x„  >•)  +  c  •  b(x,x„  y,  0  +  •  •  • 

v(x.  y,t)  a  ^  •  V,  (X,  ,y)  +  ev(x,x,,y,/)+-* 

/2»r  a[(y-l)M*r  /t,(x,.y)+e  Ar(x,.y./) 


(5) 


In  (5),  e  is  a  small  parameter  characterizing  the  di.sturbance  amplitude,  I  ^  is  the 
normalized  viscosity  and  the  temperature-viscosity  dependence  is  approximated  by  the 
power  law  with  the  exponent  n.  With  the  exception  of  the  vertical  velocity  v,  these 
expansions  resemble  those  of  Bush  [11]  for  steady  flow.  Since  the  disturbance  scale  is  of 
the  order  of  boundary  layer  thickness,  the  pressure  perturbation  is  not  constant  across  the 
layer  and  induces  the  vertical  velocity  perturbation  v  s  0(1) . 

Substitution  of  (5)  into  Navier-Stokes  equations  and  linearization  gives  a  linear  system. 
In  particular,  the  continuity  equation  yields 
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From  the  x -momentum  equation 


The)' -momentum  equation  gives 


^o.+o(|]+0(i’). 


From  the  energy  equation 


(7) 


(8) 


In  (9),  the  similarity  parameter  R  is 


(9) 


1710  quantities  (7, ,  C7] ,  and  Gy  are  viscous  dissipation  terms  that  can  be  neglected 
eveiywhere  in  the  boundary  layer.  Viscous  terms  proportional  to  the  highest  y -derivative 
are  included  in  (7)  and  (9),  since  they  are  dominant  in  viscous  sublayers  as  will  be  shown 
below.  Terms  proportional  to  8  describe  nonparallcl  effect.s  associated  with  downstream 
boundary  layer  growth.  The  small  parameter  ^  can  be  expressed  as  (see  Bush  [11]) 


6 


(II) 
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From  (10)  and  (11) 


=  (12) 

y-i[(r-i)A#»r 

Substitution  of  (U)  into  (12)  gives5  =  l//?  and  the  result  that  the  similarity  parameter 
i?  -» «•  for  the  strong  interaction  case. 

For  stability  considerations,  the  disturbance  characteristics  are  represented  in  the  normal 
mode  form 

«(jc.jt,.y./)  =  (<?o(jCpy)+5-4,(x„y)+...lexpj^^Ja(x,)<iri-fa>fj  (13) 

Modes  with  phase  speed  c  =  a)/a  satisfying  the  condition  Re(c)  <  1  will  be  considered. 
U  c-^1.  the  disturbance  is  concentrated  in  the  transitional  layer  located  near  the 
boundary  layer  edge.  This  case  has  been  studied  in  [6]  -  [8].  Substitution  of  (13)  into  (6)  • 
(9)  gives  the  following  system  of  stability  equations  to  the  first  order  approximation  with 
re.specttod 


,  ..  ,  ■  r~i  I  ^ 


/p.(o£/,-<»)v,+Z^^  =  0 


Since  the  mean'flow  pressure  is  constant  across  the  boundary  layer,  it  is  convenient  to 
nom^ize  (he  pressure  to  P,  so  that  the  equation  of  state  can  be  expressed  as 

T,p,  ml  (18) 


Using  the  noutions 


«  ^0  ^  Po  g%  M 

f  m  U^,<p  m  p  *  2o* 


we  get  the  following  system  of  stability  equations  in  the  form  convenient  for  numerical 
integration 
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(20) 

(21) 

(22) 

(23) 

where  prime  denotei  they -derivative.  Thi.s  system  is  a  hypersonic  analog  of  the  Dunn- 
Lin  system  for  the  supersonic  boundary  layer.  Derivation  of  (20)-{23)  is  based  on  an 
asymptotic  analysis  similar  to  that  conducted  by  Cheng  [12],  and  Alekseev  [13]  for 
moderate  speeds.  The  asymptotic  structure  of  disturbances  across  the  boundary  layer  is 
schematically  shown  in  Fig.  2  for  0  <  Re(c)  <  1 . 


inviscU  sMock  laytr 
boTittdary  laytr 


vrntcidxont  2 


inviaeidzopu  I 
■  viscous  Oiblaysr 


X 

Fig.  2  Asymptotic  structure  of  the  disturbance  field. 

The  following  observations  are  relevant: 

1.  In  Zones  1  and  2  with  the  thickness  0(1),  viscous  terms  are  negligible,  and  the  system 
(20)*<23)  is  reduced  to  the  inviscid  stability  equations. 

2.  In  a  viscous  sublayer  of  thickness  Sf  ■  (a/?)’’'* ,  viscous  terms  from  (20)  and  (23)  are 
of  the  order  of  inviscid  terms.  Here,  the  solution  is  expanded  in  a  power  series  of  small 
parameter^. 
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3.  In  the  critical  layer  of  thickness  =  (cicR)"''\  viscous  terms  are  also  of  the  order  of 
inviscid  terms.  This  layer  is  located  near  the  critical  point  y,:C/(y,)  =  c.  Here,  the 
solution  is  expanded  in  a  power  series  of5  2  • 

4.  hi  the  inviscid  shock  layer  located  outside  the  boundary  layer,  viscous  terms  are 
negligible  and  the  system  (MHI?)  is  reduced  to  the  inviscid  stability  equations. 

Listead  of  solving  stability  equations  in  each  sublayer  and  performing  matching  between 
asymptotic  solutions,  we  can  solve  the  system  (20)*(23)  numerically  keeping  the  viscous 
terms  with  the  highest  y-derivative  throughout  the  entire  region  of  integration.  This 
allows  us  to  obtain  a  composite  solution  which  is  uniformly  valid  in  all  decks.  A  similar 
^proach  has  been  applied  by  Mack  [2],  Gaponov  and  Maslov  [14],  Tumin  and  Fedorov 
[15],  and  others  to  stability  analysis  of  supersonic  flows.  Note  that  the  system  (20)-(23) 
can  be  directly  derived  from  the  Dunn-Lin  system  by  applying  the  hypersonic  limit  and 
using  the  asymptotics  (5). 

Summarizing,  we  conclude  that  the  disturbance  structure  and  stability'  equations  in  the 
hypersonic  boundary  layer  are  qualitatively  similar  to  those  for  supersonic  flow. 
However,  the  behavior  of  the  disturbance  amplitude  near  the  upper  boundary  layer  edge 
and  within  the  inviscid  shock  layer  is  essentially  different  from  the  moderate  sp^  case. 
The  mean  flow  temperature  tends  to  zero  at  the  upper  boundary  layer  edge^  which 
separates  the  hot  flow  in  the  boundary  layer  from  the  cool  flow  in  the  shock  layer.  As 
will  be  shown  below,  the  large  difference  between  temperatures  of  these  flows  strongly 
affects  boundary  conditions  for  the  disturbance  amplitude. 


3.3  DIaturbaneea  In  the  Inviscid  shock  layer 


Consider  the  inviscid  flow  between  the  shock  wave  and  the  upper  boundary  layer  edge 
(see  Fig.  1).  We  specify  nondimensional  variables  of  inviscid  flow  as 


-  p'-pl  f  jf"  •  P* 


The  shock- wave  shape*  is  defined  as 


(24) 

(25) 


(26) 


Subscript  s  is  used  interchaiigeubly  herein  to  denote  shock  quantity  or  (steady)  mean  flow.  This  distinction 
should  be  obviated  by  the  useage  context. 
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The  inviscid  flow  in  the  shock  layer  is  described  by  the  system  of  equations 

dp  dpu  ■ 

A,  “ 


{ 

J do  ^  do  ^  do'' 


du  ^du  .<3m1  ^  - 


p"  ’ 

where  S  is  the  entropy  function. 

The  boundary  conditions  at  the  shock  wave  are 

2 

v(jr,.y„./,)  =— 5D,| 


p(^\*yU’^\^  ~  y  _  I 

^  d  d 
■  at/  ar, 


We  assume  that  the  shock  wave  has  the  shape 

where  x  =  x,  /S ,  and  t  =  tj6  arc  fast  variables  of  the  disturbance  flowfield. 
the  boundary  conditions  (32)-(35)  arc  written  as 

2 

v  =  -^<5l^o^+£Og,+...J 


(27) 

(28) 

(29) 

(30) 

(31) 


(32) 

(33) 

(34) 

(35) 

(36) 


(37) 

Using  these, 

(38) 

(39) 
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These  conditions  have  been  shifted  to  the  undisturbed  shock  using  the  following 
expansions  for  the  steady  flow  profiles 


F'(y,)  =  F,(go)  +  e5-^(go)g,, 
where  F  signifies  a  flow  dependent  variable. 

Since  the  coefficients  appearing  in  (46)-{49)  are  functions  of  slow  variables  only,  we  can 
apply  the  WKB  method.  Accordingly,  a  solution  is  represented  in  the  form 


.  7)  +  <^1 1  (^t  •  >)+•  •  •  1  •  }adx,+Aify-  iVuf  j 


In  the  first  order  approximation  with  respect  to  the  small  parameter^ ,  the  sy.stcm  of 
equations  for  pressure  and  vertical  velocity  are 

^[i(a-fi))  +  iAv,]-^+iAv,o  =  0  (54) 

^l/(a-<u)  +  iAvJv,„+/A-^  =  0  (55) 

'ij  Pt 

(56) 

Equations  (47),  and  (50)  for  the  entropy  function  5,  and  longitudinal  velocity  u,  can  be 
found  separately  using  solutions  of  (54)  and  (55).  At  some  point  (x,,y) ,  the  system  (54)- 
(56)  has  a  non-trivial  solution  if  the  wave  numbers  satisfy  the  dispersion  relation 


A*  - — ((a  -  (u)  +  Av,  f  *  0, 


which  gives 


a-<o 


Accordingly,  the  general  solution  is 


Pio  =  P,  C,  -exp  “  j  \dy  +  C,  •  exp  ^  J  X^dy 
\  >»  A  V  n 


C, - — exp  -rfAjiiy  +Cj — — exp  •^fAiriy 
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where  y,  corresponds  to  the  upper  boundary  layer  edge.  The  constants  C, ,  and  must 
be  found  from  the  boundary  condition  (53)  and  matching  with  the  boundary  layer 
solution. 

3.4  Matching 

In  the  inviscid  zone  2  of  the  boundary  layer  (see  Fig.  1),  the  system  of  equations  (20)- 
(23)  can  be  reduced  to  the  inviscid  equations  for  the  pressured,  and  the  vertical  velocity 
Vo  =  cup 


■■‘“■SI''" 


(59) 

(60) 


Near  the  upper  boundary  layer  edge  U,(y)  I ,  and  7)  ->0,  as  y,  -0.  the  limiting 
form  of  (59),  and  (60)  is 


/  « 

0  4-/ - ar  =  0 

y 

(y-i)7;^ 


(61) 

(62) 


The  solution  of  (61),  (62)  is  represented  in  the  WKB  form 

>  / 


a-o) 


(63) 


.  V(y^ 
- 


JO 


a~o} 


Aexra  ~if  -j''-’^'«*dy  , 

a  J 


-  Aj  exp 


.Sy  dy 

J  ./tv -nr 


liV(y-i)r,  'j 


where  A,,  A,  are  constants.  Using  the  definition 


ra~  =  (r-l)A/*T,. 


(64) 


we  obtain  the  following  explicit  forms  for  the  disturbance  amplitude,  an  y  y,  *-  0 
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Now  we  consider  the  shock  layer  solution  (58)  as  y-*y,+0.  Since  the  mean  flow 
temperature  tends  to  infinity,  and  the  mean-flow  vertical  velocity  v,  -00),  the 
asymptotic  form  of  the  solution  (58)  is 


Using  the  scaling  (44),  (45),  and  the  definition  (56),  wc  obtain  the  relation 

Accordingly,  the  vertical  velocity  amplitude  v,  ■  5^,,  and  the  pressure  amplitude  arc 
respectively 


Comparison  of  the  outer  asymptotic  limit  (65),  (66)  with  the  inner  limit  (69),  (70)  and 
accounting  for  the  relation  p,{y,)^P,  shows  that  the  pressure  and  vertical  velocity 
disturbances  match  at  the  upper  boundary-layer  edge,  y-y,. 
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On  comparison  of  the  boundaiy-layer  flow  field  representations  (5)  with  the 
corresponding  shock-layer  flow  field  expressions  (42)-(4S)  and  if  the  disturbance 
characteristics  are  norm^ized  in  such  a  way  that  the  pressure  amplitude  is  p  =  0(1)  in  the 
boundary  layer,  then; 

1.  p  =  0(1)  in  the  shock  layer. 

2.  The  temperature  amplinide  r=0(l)in  the  boundary  layer,  and  T  =  0(6^)  in  the 
shock  layer. 

3.  The  longitudinal  velocity  amplitude  is  u  =  0(1)  in  the  boundary  layer,  and  u  =  0(6^) 
in  the  shock  layer. 

4.  The  vertical  velocity  amplitude  is  v  =  0(1)  in  the  boundary  layer,  and  v  =  O(^)  in  the 
shock  layer. 

This  statement  is  schematically  illustrated  in  Fig.  3. 


Rg.  3  Orders  of  magnitude  for  the  amplitude  of  the  velocity,  pressure  and  temperature 

disturbance,  (normalized  with  respect  toe ). 

Consequently, ,  the  boundary  conditions 

Vo(y,)»a<p(y,)«0  (71) 

hold  to  00 ) . 

On  the  wall  surface  the  no-slip  conditions  are 

/(0)  =  i»(0)  =  fl(0)-0.  (72) 

Equations  (20)-<23)  with  the  boundary  conditions  (71),  and  (72)  consititute  an  eigenvalue 
problem.  To  00 ) ,  this  problem  does  not  depend  on  the  shock-layer  flow.  According  to 
(71)  the  upper  boundary-layer  edge  separating  the  hot  viscous  flow  from  the  cool  sliock- 
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layer  flow  behaves  as  a  solid  wall.  In  spite  of  the  strong  interaction  between  steady  flows 
the  disturbances  with  the  phase  speed  less  than  unity  only  weakly  interact  with  each  other. 
However,  the  eigenvalue  problem  implicitly  rather  than  explicitly  depend.s  on  the  type 
of  interaction,  because  the  stability  equations  include  the  steady  flow  profiles. 


3.S  Weak  Interaction  eaae 

For  weak  interaction,  the  inviscid  flow  parameters  tend  to  their  freestream  values  as 
y-^oo.  Then  the  representations  of  the  pressure,  density  and  temperature  in  the 
boundary  layer  are  [16] 

Hx.  y.t)  =  ^^[^1  (Xi)+€-  p(x,x, ,  >,/)]+•  •  • 

P(.x,  y,t)=  - [p,  ixj ,  y)  +  e  •  p(x,  X, .  y ,  r)]+-  •  •  (73) 

t(x,y,t)  =  iy  -l)M^[T,(x^,y)  +  e^Tix,x^,y,t)]+^ 


From  an  analysis  similar  to  the  strong  interaction  case,  we  find  that  the  stability  equations 
coincide  with  (20)-(23).  if  the  similarity  parameter  is  defined  as 


SVIp:  1 

Pi 


(74) 


From  the  expression  for  the  boundary  layer  displacement  thickness  [  16] 


(75) 


the  explicit  form  of  the  small  (interaction)  parameter  is 


•#1 


(76) 


From  (74)-(76)  we  obtain  the  relation  5  =  I  /  f? ,  which  is  similar  to  the  strong  interaction 
case.  To  0(5 ) ,  the  boundary  conditions  for  the  disturbance  at  the  upper  boundary  layer 
edge  are  given  by  (71).  Similar  to  the  strong  interaction  case,  the  upper  boundary  layer 
edge  separates  a  cool  free  stream  from  a  hot  flow  in  the  boundary  layer,  and  behaves  as  a 
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solid  wall  reflecting  the  boundary  layer  disturbance.  The  eigenvalue  problem  (20}*(23) 
with  boundary  conditions  (71)»  (72)  is  valid  for  the  weak  interaction  case,  if  we 
simply  replace  the  mean  flow  profiles  by  new  ones.  Consequently  wc  can  study  the 
boundary-layer  stability  with  respect  to  disturbances  of  phase  speed  c  <  1  for  strong, 
weak,  and  moderate  interaction  regimes  using  the  same  algorithm. 


3.6  Three-dimenalonal  case 


The  foregoing  analysis  has  been  extended  to  the  case  of  3-D  disturbances  in  a  3-D 
hypersonic  boundary  layer.  Here,  the  disturbance  is  represented  in  the  normal  mode  form 


(77) 

where  z  is  the  lateral  (spanwise)  coordinate. 

For  this  case,  the  3-D  system  of  stability  equations  is  reduced  to  a  2-D  system  by  the 
following  transformation  of  the  mean-flow  and  disturbance  characteristics 


^(x,y.2.r)  =  ?o(x,.y)cxp 


+ipz-iax  , 


e  » I  (78) 

r  1  /  ax  '’o 
/  = -(OMo  +  ^o).(p  *  J.C *  J, 


where  W,,  andwo  are  the  mean  value  and  the  disturbance  amplitude  of  z -velocity 
component  respectively. 

Equations  (20)-<23),  and  the  boundary  conditions  (71),  (72)  take  the  form 


-c)/ 

(79) 

(80) 

f  • 

(81) 

.  (y  -  iX,.' + if  )+.(£/,  -  c)-^ 

(82) 

/(0)ap(0)«d(0)=o 

(83) 
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f(.y,)  =  (p(y,)  =  B(y.)=o  (84) 


Note  that  this  system  may  be  used  to  study  cross  flow  instability  at  hypersonic  speeds.  In 
contrast  to  the  full  Navier-Stokes  stability  equations,  the  system  (79)-(82)  docs  not 
contain  second  derivatives  of  the  mean-flow  profiles.  This  saves  computational  time  in 
otherwise  time-consuming  und  intense  calculations  of  the  basic  flow  distributions. 

3.7 Interactive  steady  flow 

Strong  interaction  mean  flow  profiles  for  f/,  and  T,  can  be  obtained  using  a  self-similar 
solution  of  the  boundary  layer  equations  (see  for  example  [II],  and  [17]).  We  consider 
the  boundary  layer  on  a  flat  plate  with  constant  surface  temperature.  For  brevity  we 
assume  that  the  exponent  of  the  viscosity-temperature  law  is  n  =  1 .  On  introduction  of 
the  similarity  variables 


V,WF'mH,(ri)=jC(n)  (86) 

TArn  =  ^C-F-').  (87) 

the  boundary  layer  equations  become 

y  _  1 

F'"  +  FF''+-^(G-F'*)  =  0  (88) 

f (0)  -  f '((D .  0.  cm  .  a,  (90) 

F'(«)  =  1,G(«-)  =  I.  (91) 

For  strong  interaction,  the  pressure  is 

Vr*'*.  (92) 
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where  the  constant  5o  is  obtained  from  the  interaction  law.  With  the  Tangent-Wedge 
approximation',  this  constant  is  [17] 


The  system  of  equations  (88H91)  was  integrated  using  a  fourth  order-accurate  Runge- 
Kutta  method.  Distributions  of  the  velocity  ,  and  temperature  T,iri)  are  shown  in 
Fig.  4.  The  calculations  were  performed  for  an  adiabatic  wall;  the  specific  heat  ratio 
is  7  =  1.4 ,  and  the  Prandtl  number  is  Pr  =  0.72 .  Integrating  the  temperature  profile,  we 
obtain  from  Eq.  (93)  that  A  =  1.463 . 


Fig.  4  Self-similar  profiles  for  velocity,  [/.(t?)  and  temperature,  T.(n)  in  a  hypersonic 
boundary  layer  on  an  adiabatic  flat  plate;  y  »  U ,  Pr  ■  a72,  «  =  1 . 


Hypenonte  Small 

Djnufbwce  Tb«)^.  Numerical  .wiuliom  cun  tive  other  catimalcK.  However,  neither  will  chanse  th« 
ctientui  liability  behavion  diKuucd  here. 


3-18 


SC71062.FTR 


3.8  Stability  computatlona  and  naulta 

To  obtain  the  results  to  br;  discussed,  the  eigenvalue  problem  (79)-(84)  was  solved 
numerically  using  the  strong  interaction  mean  flow  profiles  (86).  (87).  The  system  of 
stability  equations  is  transformed  into  the  form  convenient  for  integration 

^  =  A(k.S,U„T,.R)Z  (93) 

drj 

2,(0)  =  2,(0)  =  Z,(0)  =  0 

2,(-)  =  2,(o*)»2j(«)  =  0 


where  =  the  eigenvector.  According  to  the  definition  (85)  of 

Jj ,  the  similarity  parameter .  frequency  and  wavenumber  are  rescaled  os 


li=rR,oj=r<ij,k  =  rk 


(94) 

(95) 


Fourth  order-accurate  Runge-Kutta  integration  of  the  boundary  problem  (93)  was 
performed  u.sing  an  the  onhogonalization  technique.  For  the  spatial  stability  problem,  the 
eigenvalue  k(0)  is  searched  at  fixed  real  frequency  m  by  the  Newton  iteration  method. 
Unstable  regime  corresponds  to  Im(r)  <0. 

Calculations  were  performed  for  2~D  eigenmodes  with^aO  and  the  phase  speed 
Re(c)<l.  The  mean  flow  corresponds  to  the  adiabatic  case  for  specific  heat  ratio 
y  =  IA,  and  Prandtl  number  Pr  =  0.72  (set  Seaion  3.7).  The  exponent  of  the 
temperature-viscosity  law  is  n^l. 

In  Fig.  5,  distributions  of  Rcf^)  are  shown  as  a  function  of  the  frequency  ^  for 
Reynolds  number  7 «  2795 .  It  is  seen  that  spectrum  contains  two  families  of  normal 
denumerated  as  1,2,...,  and  r,2',....  At  the  intersection  points, 
Re(f],)  =  Re(f]^),  the  phase  speeds  of  two  modes  coincide  implying  their 
sytrehronization.  This  pattern  is  quantitatively  similar  to  those  obtained  by  Fedorov  and 
Gushchin  [4],  (5]  from  the  inviKid  stability  analysis  of  the  short-scale  disturbances  in  a 
shock  layer.  According  to  this  model,  the  family  1,2, ...  in  Fig.  S  corresponds  to  the  upper 
waveguide  located  between  the  boundary  layer  edge  and  the  second  turning  point 

where  a(y)»^(/-l)7,(y)  it  the  normalired  local  sonic 
speed.  The  family  r,2',...cone8ponds  to  the  lower  waveguide  located  between  the  wall 
and  the  first  turning  point  y„:U,(y„)»Re(c)-a(y„).  The  disturbance  pattern  is 
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schematically  shown  in  Fig.  6  for  the  short-wave  limit.  Large  number  modes  behave  as 
acoustic  waves  with  small  damping  due  to  viscous  dissipation. 


Fig.  5  Distributions  for  the  lower  waveguide  modes  m»  VX . and  for  the 

upper  waveguide  modes  m  »  lA— ;  if  »  2795 . 


Fig.  6  Schematic  pattern  of  disturbance  flowfield  in  the  boundary  layer. 

Near  synchronism  points  however,  the  upper  waveguide  wave  Interacts  with  the 
conesp^ing  lower-waveguide  wave  throu^  the  quiet  zone  depicted  in  Fig.  6.  In 
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accordance  with  asymptotic  model  [4],  [5],  this  interaction  causes  a  splitting  of  the 
dispersion  curves  k(0))  that  can  lead  to  instability. 


Fig.  7  The  growth  rate  a  =  -  Im ?,(<»)  for  the  mode  1  in  the  vicinity  of  the  synclwonwm 
point  1-r;  ^*2795. 

As  an  example.  Fig.  7  shows  a  spike  of  the  growth  rate  <7  =  ~Inik^(af)  near  the 
synchronism  point  1-1'.  For  this  case,  the  eigenfunction  of  the  unstable  mode  I  is 
shown  in  Figures  8a-8c.  It  is  seen  that  the  pressure  waves  are  excited  in  both  waveguides 
and  interact  with  each  other  through  the  quiet  zone. 

The  temperature  distribution  clearly  indicates  the  presence  of  critical  and  viscous 
sublayen  (see  Fig.  8c).  In  accordance  with  qualitative  pauem  of  disturbance  structure 
shown  in  Fig.  2,  viscous  sublayers  are  separated  by  a  relatively  thick  inviscid  zone.  In  the 
critical  layer,  viscous  effects  cause  strong  temperature  perturbations  and  very  high 
gradients.  The  velocity  and  temperature  amplitude  quickly  vanish  as  rj  ,  whereas 
the  ptvsiure  disturbance  is  of  the  order  of  0(1)  at  the  upper  boundary  layer  edge.  This 
feature  looks  more  graphic  and  rralistic  if  the  phj’sicaJ  coordinate  y  is  used  instead  of  the 
coordinate  7J  which  is  highly  stretched  near  the  boundary  layer  edge. 
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uppe^  wavefiuida 


g.  8b  Continued;  the  vertical  velocity  v(t7) . 
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Im(k) 


Fig.  9  The  Mode- 1  increment  Imife,  (m)  as  a  function  of  frequency;  R  =  2795 . 


Fig.  lOa  Eigenfunction  of  Mode  I  at  the  frequency  ar  =  0.39,  iR=2795, 
i  »  2i{  1 3  - 1  ■  0.02 ;  the  pressure  >r  (»]) . 
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Fig.  lOd  Gsntinued;  the  temperature  ^(17) . 


For  the  modes  depicted  in  Fig.  5,  the  most  unstable  is  Mode  t.  Its  increment  lmiF,(a))  is 
shown  in  Fig.  9  for  the  Reynolds  number  X  » 2795 .  It  is  seen  that  the  growth  rate 
achieves  its  maximum  at  the  frequency  oT  =  039 .  The  corresponding  eigenfunctions  are 
given  in  Figures  10a  -  lOd.  For  this  case,  the  pressure  amplitude  is  relatively  large  in  the 
upper  waveguide.  The  lower  waveguide  does  not  exist  because  its  upper  boundary  is 
y„  <0.  Since  the  phase  speed  is  small,  the  critical  layer  is  merged  with  the  wall  sublayer 
and  blends  to  one  deck  as  clearly  seen  in  Figure.s  lOc,  and  lOd.  We  believe  that  the  Mode* 
1  instability  is  due  to  reflection  of  the  pressure  disturbance  by  the  upper  boundary  layer 
edge.  For  this  reason,  we  refer  to  this  mode  as  a  “rcflectional  instability."  Distributions  of 
its  wavenumber  Re^fR)  and  growth  rate  ff(^)  =  at  the  frequency  ST  *  039 

are  shown  in  Figures  1  la  and  I  lb  respectively.  The  rcflectional  instability  has  a  relatively 
small  phase  speed  Rc(c)»0.14,  which  is  typical  for  acoustic  modes.  For  reference  we 
recall  that  the  vorticity  mode  located  in  the  transitional  layer  has  the  phase  speed 
Re<c)  I  [6].  From  d^  shown  in  Figures  9,  and  1  lb  we  conclude  that  the  rcflectional 
instability  has  a  relatively  high  growth  rate,  of  the  order  of  2  10'*.  Such  a  strong 
instability  may  achieve  the  critical  amplification  necessary  for  transition. 
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Since  the  eigenvalue  problem  (79)  -  (84)  is  isolated  from  the  outer  shock  layer,  we  can 
assume  that  similar  mechanism  of  the  reflectional'  instability  exists  in  hypersonic 
boundary  layers  at  the  moderate  and  weak  interaction.  Analysis  of  these  cases  is  in 
progress. 

3.9  Dfseusalon 

We  have  investigated  the  instability  of  the  flat-plate  hypersonic  boundtiry  layer  to  the 
normal  modes  with  wavelength  of  the  order  of  boundary  layer  thickness  and  phase  speed 
less  than  1.  These  inodes  are  associated  with  acoustic  waveguides  located  in  the  bounds 
layer.  Asymptotic  analysis  showed  that  in  the  first  order  approximation  with  respect  to  the 
small  parameter^,  the  upper  boundary  layer  edge  behaves  as  a  solid  wall  which  reflects 
the  boundary-layer  disturbance.  In  spite  of  strong  interaction  between  the  steady  shock 
and  boundary  layer  mean  flows,  the  unsteady  disturbances  under  consideration  only 
weakly  interact  with  each  other.  From  asymptotic  analysis  for  the  weak  interaction  case 
we  have  found  that  the  eigenvalue  problem  has  the  same  form  for  both  weak  and  strong 
interaction.  This  allows  as  to  study  the  boundary  layer  stability  for  the  strong,  moderate 
and  weak  interaction  regimes  using  a  single  approach.  The  eigenvalue  problem  is  easily 
extended  to  the  3-D  case  and  can  be  ased  for  analysis  of  the  cross  flow  instability  in  three 
dimensional  hypersonic  boundary  layers. 

From  our  eigenvalue  problem  solution  for  the  strong  interaction  case,  we  conclude  that 
the  reflection  of  disturbance  by  the  upper  boundary  layer  edge  leads  to  instability  of  the 
upper-waveguide  modes.  We  refer  to  this  mechanism  as  the  reflectional  instability. 
Similar  waveguide  mechanisms  have  been  ascribed  to  the  dramatic  effect  of  thinning  the 
shock  layer  fory-»lin  [9].  There,  thinning  produces  downstream  amplification  of 
otherwise  long  wavelength  fluctuations  that  decay  downstream  for  thicker  .shock  layers. 

Calculations  have  shown  that  the  reflectional-instability  mode  has  a  relatively  small  phase 
speed  c*O.I4.  This  is  consistent  with  recent  experiments  of  Maslov  and  his  group, 
performed  on  a  flat  plate  in  hypersonic  wind  tunnel  at  Mwh  20  (private  communication). 
They  observed  slow  waves  propagating  in  the  strong  and  moderate  interaction  region  with 
a  phase  speed  of  the  order  of  0.2,  and  growing  downstream.  Since  the  reflectional- 
instability  mode  has  a  large  growth  rate,  it  can  compete  with  other  unstable  modes  such, 
as  the  vorticity  mode.  We  expect  that  downstream  amplification  of  this  mode  can  lead  to 
early  transition  on  hypersonic  vehicle  surfaces.  These  findings  motivate  us  to  study  the 
reflectional  instability  mechanism  further. 

We  also  found  that  the  upper  waveguide  modes  can  synchronize  with  the  lower 
waveguide  modes.  Near  synchronism  points,  additional  growth  rate  spikes  have  been 
found.  This  instability  mechanism  is  similar  to  those  obtained  from  qualitative  analysis  cf 
a  simplified  flow  in  a  shock  layer  [4],  [5].  In  the  vicinity  of  synchronism  points  the 
normal  mode  decomposition  is  not  valid  due  to  the  singular  nature  of  the  dispersion 
relation.  This  situation  has  been  studied  by  Fedorov  and  Khokhlov  [3]  for  the  cose  of 
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moderate  Mach  numbers.  Using  their  results,  we  assume  that  the  nonparallel  effects*  may 
cause  strong  exchange  between  modes  near  synchronism  point. 
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4.  Low-frequency  instability  of  hypersonic  Interactive  boundary  layers 
on  a  cooled  surface 


4.1  Introduction 

This  chapter  reviews  research  supported  by  the  contract  and  other  funding  concerning 
propagation  of  unsteady  three-dimensional  perturbations  in  hypersonic  interactive 
boundary  layers  on  a  cooled  surface.  Such  high  speed  boundary  layers  that  are  subject  to 
the  combined  effect  of  the  interaction  with  a  shock  layer  and  cooling  of  the  surface  are 
typical  of  realistic  flight  conditions  and  lead  to  substantial  heat  transfer  and  drag.  The 
geometry  of  a  flat  plate  is  considered  herein.  For  this  simple  shape,  the  stability  of  even 
the  steady  flow  is  complicated  by  a  multiplicity  of  different  mechanisms  not  exhibited  at 
moderate  supersonic  Mach  numbers.  Two  classical  limits  naturally  evolve  at  high  Mach 
numbers  in  this  flow.  Near  the  leading  edge,  a  strong  interaction  regime  occurs  in  which 
the  shock  wave  is  not  far  from  the  surface  and  the  viscous  boundary  layer  and  inviscid 
.shock  layer  are  strongly  coupled.  Further  downstream,  the  effect  of  the  shock  wave 
disappears  and  the  boundary  layer  develops  under  zero  pressure  gradient  (weak 
interaction).  Strong  and  weak  interactions  have  been  studied  extensively  (see  for 
example,  [1-3].  Transition  to  turbulence  has  been  observed  in  the  last  region  only. 
However,  as  the  Mach  number  increases,  the  strong  interaction  region  expands  rapidly. 
This  could  significantly  influence  the  early  stages  of  transition,  since  it  will  affect 
important  receptivity  and  linear  amplification  processes. 

In  the  limit  of  infinitely  large  Mach  number,  the  temperanire  dLstfibution  has  increasing 
nonunifoimity  and  cannot  be  described  in  closed  form  independent  of  the  Mach  number. 
Additional  distinctions  are  made  in(4].  These  include  the  observations  that  the  domain  of 
importance  is  not  semi-infinite  as  in  lower  speed  regimes,  the  shock  is  strong  and  curved 
and  non  parallelism  can  be  more  important  and  typical.  Accordingly,  the  Orr-Sommerfeld 
approach  cannot  be  applied  a  priori,  and  a  suitable  hypersonic  limit  involving  the  Mach 
and  Reynolds  numbeii  needs  to  be  applied.  Different  types  of  instabilities  can  be  studied 
using  appropriate  scalings  relevant  to  the  characteristic  wave  lengths  present. 

An  interesting  type  of  instability  is  inviscid  which  is  due  to  a  generalized  inflection  point 
within  a  thin  transitional  layer  near  the  boundary  layer  upper  edge.  It  was  analyzed  for  the 
moderate  interaction  regime  in  [S]  where  the  transition^  layer  between  the  external  flow 
and  a  shear  layer  was  considered  and  in  [6]  for  a  gas  with  a  Chapman  viscosity  law  and  a 
unit  Prandtl  number.  The  main  conclusion  of  these  studies  is  that  there  is  a  numerous  set 
of  relatively  weakly  unstable  acoustic  modes  and  a  most  unstable  vorticity  mode.  The 
acoustic  modes  were  also  analyzed  in  [7]  and  [8].  The  phase  speed  of  all  the  modes  is 
close  to  the  external  flow  velocity  so  that  the  relative  growth  is  small.  The  results  are  in 
good  qualitative  and  quantitative  agreement  with  the  computations  for  finite  Mach 
numbers  in  [9].  The  more  general  case  of  arbitrary  Prandtl  number  and  power  law 
viscosity  was  considered  in  (lOJ  where  the  effect  of  the  viscosity  on  the  stability 
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characteristics  was  included.  Continuation  into  the  strong  interaction  zone  was  treated 
in[ll]  assuming  the  Sutherland  formula  for  the  viscosity  and  Prandtl  number  equal  to 
unity.  Weak  and  strong  interaction  regimes  were  considered  and  it  was  conclude  that 
although  the  scaling  is  quite  different,  the  qualitative  features  of  the  vorticity  mode 
instability  are  the  same.  The  interaction  leads  to  differences  in  the  unperturbed  flow 
proHle  only  and  does  not  change  the  physical  mechanism  of  the  instability.  No  influence 
of  the  shock  wave  arises  because  the  perturbations  are  mostly  concentrated  in  the 
transitional  region  between  the  shock  and  boundary  layers.  However,  as  shown  in  [8]  for 
the  non-interactive  boundary  layer  on  a  wedge  such  influence  may  occur  in  some 
situations  when  there  is  a  proper  correlation  between  the  wavelength  scale  and  the 
thickness  of  the  shock  layer.  Another  example  of  ampliflcation  of  the  small  time- 
harmonic  perturbations  in  ^e  curved  shock  layer  was  obtained  in  [4]  by  analysis  of  the 
initial  boundary-value  problem  in  the  downstream  direction. 

Other  types  of  unstable  waves  are  small  phase  speed  viscous  modes  described  within  the 
framework  of  triple-deck  theory  in  [12]  which  showed  that  in  supersonic  flows,  unstable 
modes  of  this  type  are  essentially  three-dimensional,  where  the  direction  of  propagation 
lies  outside  the  Mach  cone.  They  dominate  at  moderate  Mach  numbers,  but  have 
negligibly  small  growth  rate  at  hypersonic  speeds  and  inviscid  instability  prevails.  The 
last  conclusion  is  correct  if  the  temperamre  of  the  wall  is  of  the  order  of  the  stagnation 
temperature.  The  effect  of  surface  cooling  at  moderate  Mach  numbers  was  estimated  in 
[13]  in  which  it  was  found  that  it  increases  the  growth  rates  and  compresses  the  scales  of 
the  eigenmodes.  For  sufficiently  strong  cooling,  the  wavelength  becomes  comparable 
with  the  boundary  layer  thickness  and  a  new  moderate  cooling  regime  occurs.  Surface 
cooling  is  typical  for  hypersonic  flight  conditions.  Its  influence  on  the  viscous  mcxles  may 
be  significant  Special  study  of  the  strong  interaction  case  is  particularly  warranted. 

Analysis  of  the  viscous  eigensolutions  in  the  strong  interaction  boundary  layer  was 
initiated  in  [14]  for  the  steady  two-dimensional  case.  It  showed  that  their  streamwise 
evolution  is  controlled  by  a  power  rather  than  exponential  law  in  the  weak  interaction 
case.  According  to  [15],  the  exponent  in  this  power  law  grows  to  infinity  as  the  specific 
heat  ratio  tends  to  unity.  Therein,  an  analytic  formula  for  the  exponent  was  obtained  and 
it  was  shown  that  the  disturbance  field  is  governed  by  the  interaction  between  the  viscous 
region  and  the  main  part  of  the  boundary  layer,  while  the  upper  region  plays  a  passive 
role.  A  relationship  between  this  mode  and  the  Ughthill  m^e  in  the  weak  interaction 
region  was  postulated.  This  conjecture  was  subsequently  verified  by  Brown  et  al.  (1975). 
Another  important  result  is  that  the  exponent  exhibits  a  sharp  growth  when  the  surface  of 
the  plate  is  cooled.  A  similar  conclusion  was  obtained  in  [17]. 

This  effort  focused  on  three-dimensional  unsteady  disturbances.  From  the  foregoing,  it 
can  be  assumed  that  cooling  of  the  surface  will  lead  to  appearance  of  short  wavelength 
eigensolutions  in  the  unsteady  three-dimensional  us  well  as  the  steady  two-dimensional 
case.  Moreover,  these  perturbations  can  be  effectively  confined  within  the  boundary  layer. 
Because  the  local  Mach  number  within  the  boundary  layer  is  of  the  order  of  unity,  a  three- 
dimensional  instability  similar  to  that  considered  in  [12]  for  the  moderate  Mach  numbers 
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may  -be  anticipated.  In  Section  2,  the  proper  scaling  for  the  perturbations  field  is 
developed.  The  asymptotic  expansions,  analysis  of  the  each  region  and  matching  leading 
to  the  complete  formulation,  which  is  a  type  of  the  triple-deck  problem  with  a  new 
interaction  law,  are  discussed  in  Section  3.  As  a  simple  framework  to  exhibit  some  of  the 
high  Mach  number  physics  and  relate  to  the  earlier  two-dimensional  steady  effort,  the 
linearized  forced  uansieui  problem  of  a  harmonically  oscillating  hump  starting  its  motion 
at  some  instant  of  time  is  considered  in  Section  4. 

In  Section  4,  the  dispersion  relation  is  analyzed  and  all  eigenmodes  with  respect  to 
upstream  and  downstream  propagation  are  classified.  In  particular,  the  mode  of  Neiland  • 
Brown  and  Stewartson  is  shown  to  propagate  upstream.  As  a  confirmation  of 
expectations,  another  three-dimensional  mode  appears  that  is  unstable  in  the  downstream 
direction.  Its  growth  rate  is  not  bounded  in  the  short  wave  limit  so  that  the  forced  problem 
is  well-posed  on  the  set  of  functions  with  compact  support  in  wavenumber  space.  A 
limiting  process  to  a  wider  class  of  functions  demonstrates  that  the  part  of  the  solution 
containing  the  fast  growing  modes  is  convected  downstream.  Hence,  in  the  large-time 
limit  the  perturbation  field  consists  of  the  time-harmonic  component  without  an  infinitely 
growing  mode.  The  far-fteld  asymptotic  of  this  component  is  e.stimatcd  and  the  direction 
of  the  maxima]  amplification  versus  frequency  is  determined 

Other  possible  mechanisms  of  instability  and  correspondence  with  experimental 
observations  in  [18]  and  [19]  are  discus.sed  in  Section  5. 

4.2  FormuttUon 

We  consider  three-dimensional  short-scalc  unsteady  perturbations  in  a  hypersonic  flow 
past  a  semi-infinite  flat  plate  developing  at  a  distance  L  from  the  plate  leading  edge  in  a 
region  where  strong  interaction  regime  occurs.  The  Cartesian  coordinates  (x,y,z)  and 
corresponding  velocities  («,w,v)  ure  chosen  in  streamwise,  vertical  and  spanwisc 
directions,  respectively.  Freestream  pressure,  density,  temperature  and  viscosity  arc 
denoted  asp„,  p.,  r,,p.and  the  external  flow  velocity  as  I/.. 

The  specific  heat  ratio  y  and  a  Prandtl  number  a  ore  assumed  to  be  constant,  and  the 
dependence  of  the  viscosity  on  the  temperature  is  taken  in  a  power  law  form 

p/p.=(r/r.r  (2.1) 

The  mean  flow  is  basically  controlled  by  interaction  between  the  boundary  and  shock 
layers  and  allows  self-similarity.  Its  profile  shown  schematically  in  Fig.  1  includes  three 
characteristic  regions:  viscous  high-temperature  boundary  layer,  inviscid  shock  layer  and 
transitional  layer  appearing  due  to  singularity  in  temperature  behavior  near  upper  edge  of 
the  boundary  layer  [2].  As  the  plate  temperature  decreases,  the  boundary  layer  solution 
becomes  singular  at  the  will  because  the  viscosity  vanishes  there,  and  according  to 
[13,17],  an  additional  fourth  buffer  layer  occurs  near  the  surface  where  heat  and 
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momentum  fluxes  are  constant  while  the  temperature  is  of  the  order  of  the  wall 
temperature  and  changes  rapidly  together  with  the  velocity. 


Rg.]  Qualitative  distribution  of  mean  flow  temperatures  and  velocity  across  boundary 
layer. 


For  convenience,  the  following  two  small  parameters  characterizing  mean  flow  are 
introduced 


/»=/?e-‘'^W^,£  =  r^/(WX)  .  (2.2) 

where  M  is  the  Mach  number,  and  Re  is  the  Reynolds  number.  Both  of  these  parameters 
are  based  on  freestxeam  conditions  and  the  distance  L.  The  quantity  T^is  the  wall 
temperature.  For  the  strong  interaction  approximation  to  be  applicable,  the 
parameter/i  which  is  the  ..^lative  thickness  of  the  boundary  layer  must  meet  the  additional 
constraint 
hM 


4.2. 1  Preliminary  considerations  and  scales 

In  this  unsteady  framework  and  consistent  with  the  foregoing  discussion,  the 
characteristic  length  of  the  disturbances  are  assumed  to  be  small  in  comparison  with  the 
global  scale  L,  and  their  structure  in  they  •direction  consists  of  multiple  layers.  There  are 
at  least  two  different  regions:  the  main  part  of  the  boundary  layer  and  a  thin  viscous  wall 
sublayer,  hi  what  follows,  the  characteristic  length  and  time  of  the  perturbation  are 
respectively  denoted  as  Dx  and  Dt  ;velocity  and  pressure  variations  as  Du  and  Dp  all 
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values  referring  to  the  viscous  sublayer  and  to  the  main  part  of  the  boundary  layer  by 
subscripts  v  and  b ,  respectively.  For  ordcr-of-magniiude  estimations  it  is  a.ssumed  that: 

(a)  the  pressure  variation  across  the  boundary  layer  is  negligibly  small. 

(b)  streamwise  and  spanwise  scalings  are  the  same. 

(c)  the  viscous  sublayer  lies  at  the  bottom  of  the  buffer  layer  so  that  the  undisturbed 
flow  in  that  region  has  a  linear  shear  velocity  profile  whose  slope  can  be  evaluated 
by  comparison  of  the  momentum  flux  on  the  surface  with  that  of  the  boundary  layer, 
so  that 


(d)  Within  this  sublayer,  the  viscous  forces  are  in  balance  with  the  inertia  and  pressure 
gradient  terms  so  that 


!  b>x-l^pl  Px  - . 


(e)  in  the  main  part  of  the  boundary  layer  the  perturbation  field  is  mainly  inviscid  so 
that  only  inertia  and  pressure  gradient  terms  must  be  kept  in  balance.  This  implies 

Ac  -  A/»/Ar, 

(f)  interaction  occurs  basically  between  the  main  part  of  the  boundary  layer  and  the 
viscous  sublayer  and  perturbations  of  the  displacement  thickness  have  the  same 
order  in  these  layers 

(This  follow.s  f  15]  for  the  steady  case.) 

(g)  flow  in  the  viscous  sublayer  is  un.steady  and  nonlinear  with  respect  to  the  lime  and 
amplitude  scaling  of  the  disturbance 

Ax/ Ar  - 

For  completeness,  these  esu'mations  must  be  .supplemented  by  the  order  of  the  mean 

boundary  layer  values 


y»  -  Uir*  -  t/..p»  -  -  T,M\ 


as  well  as  the  equation  of  state 
p/(pT)  ■  Const. 
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These  relations  lead  to  a  unique  scaling  for  the  perturbation  field  in  the  viscous  sublayer 
and  the  main  part  of  the  boundary  layer.  To  provide  matching  between  these  regions,  an 
intermediate  layer  of  the  constant  momentum  and  heat  fluxes  is  considered.  Its  thickness 
may  be  evaluated  by  comparison  of  the  heat  flux  near  the  surface  and  with  that  of  the 

boundary  layer.  Thus, 

Also  for  matching,  a  fourth  upper  layer  located  at  the  bottom  of  the  shock  layer  is 
introduced.  In  that  region,  pressure  variations  in  the  vertical  direction  are  assumed  to  be 
essential.  Neglecting  viscous  forces,  the  thickness  of  this  layer,  Dy  is  given  by  the 
condition  that  the  local  slope  of  the  characteristics  has  the  same  order  as  the  relative 
thickness  of  the  region, 

Ay/Ax-l/Mfc,, 

where  the  local  Mach  number,  M  loc  is  connected  with  the  local  temperature  T  loc 
within  this  layer 

which,  in  turn,  may  be  explicitly  written  following  the  asymptotic  formulas  (3.02),  (3.13) 
in  [2]  reexpressed  in  terms  of  physical  variables  as 

Ts^-T,MW^y/Ui). 

This  completes  the  estimates.  As  shown  in  1 15]  for  the  upper  layer,  [17]  and  [13]  for  the 
buffer  layer,  both  layers  play  a  secondary  role  in  the  interaction  but  are  necessary  for 
rational  formulation  of  the  problem. 

Coupling  of  all  of  the  above  estimations  gives  us  the  appropriate  scaling,  which  may  be 
written  in  terms  of  the  small  parameters  h  and  e  in  the  form 

x  =  +  (2.3) 

z=lj£  Zi  , 

r  =  (L/f/.)£‘*'‘’<i, 

y=  iMe'^^^y,,  y,.y„  + 
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The  perturbation  field  is  schematically  shown  in  Fig.  2.  It  is  decomposed  into  four 
characteristic  regions:  the  lower  viscous  region  y,  =  (?(1) ,  the  buffer  layer  =  0(1) ,  the 
mean  part  of  the  boundary  layer  y,  =  0(1),  and  the  upper  layer  y^  =0(1)  located  at  the 
bottom  of  the  shock  layer  above  the  boundary  layer  upper  edge  y,  =  y^ .  Assumption  (a) 
that  the  pressure  variation  across  the  boundary  layer  is  negligible  will  be  correct  if  the 
following  restriction  on  the  small  parameters  is  valid: 

(2.4) 

i.e.  ,  the  streamwise  scale  is  greater  than  the  thickness  of  the  boundary  layer.  With  such  a 
restriction,  the  upper  deck  will  always  be  thicker  than  the  transitional  layer  at  any  realistic 
(tr.y.d)). 


Fig-  2  Asymptotic  structure  of  the  perturbation  field,  ^  =  2(1  +  2fit))(3y  -  2)(3y  - 1). 
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4.3  AsymptoUe  approach. 

From  condition  (g),  a  typical  amplitude  of  the  perturbation  within  each  layer  can  be 
formulated.  For  these,  the  independent  variables  for  the  region  arc 
respectively,  with/I  =  l,..-,4.  These  will  be  studied  in  the  next  subsections. 

4.3.1  The  viscous  sublayer  y^  *0(1) 

The  solution  in  this  layer  is  represented  in  the  form 

p  =  +  (3-') 

u  =  U..e'*“un+... 

V  *  C/-A£‘'“v„+... 

r  =  r.W'(£7’ot+-) 

P  =  P-/«*(«‘'Ao+-) 

In  the  absence  of  perturbations,  the  solution  is  a  uniform  shear  flow 

Foi  =itPio*Pio 

where  the  undisturbed  pressure  Po  and  shear  stress  C,  are  determined  by  matching  with 
the  remainder  of  the  mean  flow.  Substituting  the  expansions  in  the  cquation.s  and 
retaining  like  orders  gives 


PjCau,,  /dr, +«„*<„ /ax,  +v„au„  /ay, +w„au„  /az,« 

y*'ap„ /dx, =a^«„/ay,i 

(3.3a) 

Po(dw„  /  dr,  ^  +  v„aw„  /  ^,  +  H'„dw„  /  dz,  ♦ 

r‘‘9pii  /9«|  =  /  ay,, 

(3.3b) 

Pii  “ 

(3.3c) 

d«„  /  dxj  +  a>Vj,  /  dz,  +  dv,,  /  dx,  ■  0 

(3.3d) 
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The  boundary  conditions  on  the  surface  arc  the  no-sIip  conditions 

“ii  =  =  Wii  =  0  ai  y-j  =  0.  (3.3e) 

At  the  outer  edge  of  the  layer,  the  asymptotic  behavior  of  the  solution  is 

Uii=Ciyi+A(Xi,Zi,tt)+..., 

»*'ii  = -<9PoC,y, )■' J(9ph • 

=  -{ypoCiyxY' i(9p,i  /3z,V£x,+...  as  y,  ->  «>  (3.3f) 

To  complete  the  problem  (3.3)  a  relation  between  the  prcsjuirc  perturbations  p,,  (x,  ) 

and  the  displacement  thickness  -A(x,,2,,/,)must  be  applied.  This  relation  (interaction 
law)  will  be  obtained  after  analysis  of  the  other  decks. 

4.3.2  The  bt^er  regiony.^  =  0(1) 

Altliough  this  layer  is  still  thin  in  comparison  with  the  boundary  layer,  variations  of  the 
mean  flow  temperature  and  density  within  it  are  not  small.  Requirement  of  continuity  of 
the  pressure  and  the  appearance  of  the  asymptotics  (3.30  rewritten  in  terms  the  variable 
y2  implies  the  following  expansions  in  the  buffer  region 

psp_A*A/*(Pjo+f:'*^"pj,+..,)  (3.4) 

tt  =  0.(6tra(yj)  +  e'*"Mj,+...) 

r  =  r,Af»(£rjo(y,)  +  £'*T„+...) 

P-P^h^-'PM 

The  orders  of  the  temperature  and  density  disturbances  in  the  above  expansions  are 
estimated  from  the  energy  and  state  equations.  Substituting  these  expansions  gives  the 
problem  for  the  mean  flow  and  unsteady  perturbations.  The  mean  flow  is  characterized 
by  constant  fluxes  of  heat  and  momentum,  i.e. 

d(iT^  )“du^  /  dyj )  /  Jy,  =  0 

daT^rdT^/dy,)/dy^^O 
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Pm^2q  =  /’2o  =  Const. 

The  boundary  conditions  are 

^2o  =  *’“2o=Oatyj=0. 

This  system  leads  to  the  solution 
Tjn  =  (1  +  Q(1  +  (o)yi 

u^=(Q/C,){T^-l) 
Pia  “  Po 
Pm  ~  P>  /  T20 . 


(3.6) 


(3.7) 


The  pressure  p^,  wall  shear  stress  Q  and  heat  flux  Cj  cannot  be  specified  until  the  global 
mean  flow  solution  is  obtained.  The  next  approximation  gives  us  the  perturbation  field 
equations 


«m3«ji  f  /  dyj  =  0 

(3.8a) 

(3.8b) 

(3.8c) 

(3.8d) 

Pjo9«2,  /  +  «no9pj,  /  3.x, + 3(P2j  Vj, )  /  3y2  =  0 . 

(3.8e) 

PjI  ^  p»  ^31  ^  ® 

(3.80 

This  is  a  system  of  linear  invisctd  equations  with  the  absence  of  the  pressure  gradient  in 
the  X  -momentum  and  energy  equations.  Note  also  that  this  deck  and  all  decks  above  are 
quasi-steady  in  contrast  to  the  viscous  sublayer  as  for  other  speed  regimes.  Formally,  the 
problem  (3.8)  is  second-order  with  respect  to  the  vertical  variable  y,  (cqs.  (3.8c), (3.8c)). 
To  obtain,  a  unique  solution  it  is  therefore  sufficient  to  specify  continuity  of  the  pressure 
and  matching  of  the  y-component  of  the  velocity  with  asymptotic  (3,30  at  the  outer  edge 
of  the  viscous  sublayer.  By  direct  substitution,  it  can  be  verified  that  a  .solution  of  the 
system  (3.8)  meeting  these  conditions  is 

Pji  **  Pii('*i’^i*^^  (3.9) 
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Mjj  —  i4(Xj ,  Zj  ,/j  )C^du^  f  ^2 

^11  =  -'^JSPoU2o  r'  I  ('^1  ,/*, 

Vjj  =  — MjqCj  j  3i4  /  oX| 

r„=AC7'^/r„/^, 


•^.i.  J  A/a/n  parr  a/ the  boundary  layer  >•,  =  0(  I)  J 

Asymptotic  expansions  within  this  region  may  be  obtained  by  analysis  of  the  limiting 
form  of  the  buffer-layer  solution  (3.9)  for  large  in  tenns  of  the  variable  yj,  or  by  direct 

estimation  from  the  equations  of  motion,  assuming  that  the  pressure  perturbation  is  of  the 
same  order  as  before.  This  gives 

p  =  pyM\p2o+e'‘^‘^Pn-t-...)  (3.10) 

w  =  C/.£r'**-w„+... 

V  =  +...+V3o(yj)+...) 

r=r^*(r„(y,)+£'''‘'r3,+...) 

p=p-/»Vj(,(j's). 


where  the  out-of-order  term  Vj^  included  in  the  expansions  for  the  vertical  velocity  i.?  not 
involved  in  the  leading-order  problem  for  the  disturbances  but  is  necessary  for  a 
consistent  description  of  the  mean  flow.  The  unperturbed  boundary  layer  flow  profile  is 
locally  uniform  (no  dependence  on  the  short  scale  streamwise  variable  x, )  and  obeys  the 
system  of  ordinary  differential  equations  that  result  from  the  usual  three  quarter  power 
similarity  (sec  for  example  [20]).  A  complete  fomiulation  also  includes  the  problem  for 
the  inviscid  shock  layer.  This  complete  system  wa,s  investigated  in  detail  in  [2],  and  in 
von  Mises  variables  in  [3].  A  brief  analysi.s  of  the  problem  in  terms  of  the  physical 
independent  variables  is  given  in  the  Appendix.  Here,  we  note  only  that  the  undisturbed 
pressure  is  con.stant  and  the  temperature  and  density  distributions  are  connected  by  the 
state  et]uation 
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Pya  “  Po’Pya^y)  ^  Po  (3. 1 1 ) 

The  leading-order  problem  for  the  perturbations  has  the  form 

/  2t,  +  v,|<f«3o  / dy^  +  /''clpj,  /  3jf,  =  0  (3.12a) 

PvMy^^ix  /  2c,  +  Y~'^P^l  / 2z,  -  0  (3. 1 2b) 

3p3, /3yj=0  (3.12c) 

Pjo(M3o27’3i/2*i  +  v„dr3g/<fy3-(y-l)y‘'ujoDpj,/aJc,  =  0  (3.12d) 

Pii  ^ P'M  ^31  ^  ^30  ~  Pii^ Po‘  (3.12e) 


Like  the  problem  (3.8)  for  the  buffer  layer,  this  gives  rise  to  an  inviscid  linear  system  of 
second  order  with  respect  to  the  vertical  variat^le,  but  with  non  trivial  pressure  gradient 
forces.  Due  to  the  inviscid  nature  of  the  problem,  it  is  sufficient  to  apply  conditions  of 
continuity  of  the  pre.ssure 


PnixvZfO*  ^3.13) 

and  matching  the  vertical  velocity  component  with  the  buffer  region  solution.  It  is  .shown 
in  what  follows  that  the  upper  region  is  also  mainly  inviscid  and  the  same  continuity 
conditions  should  be  appli^  for  its  matching  with  the  main  part  of  the  boundary  layer. 
Therefore,  the  analysis  can  be  limited  to  solution  for  the  two  variables  P3i,v„  only. 
Successively  excluding  all  other  variables  by  expressing  them  in  terms  of  PutVji ,  and 
their  derivatives  gives  accordingly 

PoPm^\  / 2yj9.c,  -  PoPxidu^  /  dy^^Vj,  /  dx,  + 

r''P3o“M2 V, I  /  ^ci"  -  f'Poi^Pu  I  +  2 /  *,0=  ®  (3.14) 

The  solution  of  Eq.  (3. 1 4)  matched  with  the  buffer  layer  is 

5v3,  /3x,  =«jo(“Cr’2*^/2x,,  +(yA»)"'/{"2Vii  +r„3H,‘*Ap„Ky3) 

(3.15) 

where  the  operator  ^  is  the  two-dimensional  Laplacian,  i.e. 

Ap„  *3*Ph/2c,*+3Vh/2z,*.  (3.16) 
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The  integrand  in  the  right-hand  side  cf  the  expression  (3.15)  has  a  singularity  at 
>3  =  0  which  may  be  evaluated  from  asymptotics  (A.4)  as  Oiy^  ).  Since  it  is  a  weak 
singularity  for  positive  O) ,  the  integral  exists  and  the  solution  is  regular  at  any  y,  >  0. 


4.3.4  Tlic  upper  region  y\  =  0(1) 

The  undisturbed  flow  in  this  layer  is  an  asymptotic  “tail”  of  the  shock  layer  solution 
considered  by  many  authors  (see,  for  example,  [2,3].  These  analyses  show  that  the 
leading-order  pressure  and  velocity  are  constant  in  this  layer  while  the  temperature  and 
density  have  a  power  Jaw  singularity  at  the  lower  boundary  (y,  -♦O*).  As.suming  the 
pressure  perturbation  to  be  of  the  same  order  as  that  within  the  boundary  layer,  the 
amplitudes  of  the  other  perturbation  components  can  be  estimated.  This  leads  to  solution 
expansions  in  the  form 

P  =  (3.17) 

u  =  UJl  +...) 

p  =  ’(p«  (>3 ) + +•  •• ) 

The  above  expansions  do  not  include  perturbations  of  the  sneam-  and  spanwise  velocity 
components  since,  like  the  mean  flow  problem,  the  leading-order  equations  for  the 
perturbations  decouples  and  these  velocities  are  not  involved  in  the  problem  for 
{p,P,T,v)- 


The  mean  flow  di.stribution  may  be  obtained  in  an  explicit  form  from  results  of  the 
previously  mentioned  authors,  giving  after  some  manipulation  of  the  dependent  and 
independent  variables 

P40  =  Po  (3J8) 

v^  =  (3/4)y,o 
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Here,  the  constant  B  is  expressed  in  terms  of  the  the  boundary  layer  pressure  p^and  the 
location  of  the  shock  layer  >},  as 

R  =  {3V  (7 +  I)'’'(3y  -  2)'*  y^y, . 

In  the  next  approximation,  the  equations  for  the  perturbations  take  the  form 

P4odv^,/atc, +7‘'3p4,/ay4  =  0  (3.19a) 

dy^).iy-\)y‘*^p^J^x^=0  (3.19b) 

ap4,  /  ax, +a(p4oV4,)/ 9x4=0  (3.i9c) 

Pa\^  Pii^  Po  (3.19d) 

This  system  is  reduced  to  a  second-order  equation  for  the  pressure  with  boundary 
conditions  of  continuity  of  the  pressure  at  the  lower  boundary  and  a  radiation  condition  at 
the  upper  one.  This  leads  to  a  standard  hyperbolic  problem  in  a  nonuniform  medium 

a*p4, !  ax,j  -  3(7409^4,  / 8x4) /  3x4  =  0  (3.20a) 

P4i=Pii  a*  X4*0  (3.20b) 

radiation;  {9  /  9x,  +  V  / Sy,  )/;4,  0  as  X4  .  (3,20c) 

A  forma]  solution  of  this  problem  may  be  obtained  by  application  of  the  Fourier 
transformation  with  respect  to  x, ,  solution  of  the  resulting  ordinary  differential  equation, 
and  inversion.  If  the  pressure  on  the  lower  boundary  is  decays  sufficiently  rapidly  in  the 
upstream  direction,  the  solution  exists  and  its  explicit  form  is 

=  (3.21) 

The  kernel  function  G  has  the  integral  representation 
G(5.74)  =  ^(1  - 10(1  -  ff-'rO -?)sin(2»?) . 

Jo>J'/2)V,(;t:7>  ‘Vi/it.  (3.22) 

where 

r"*(l  +  r/2)‘'yi*"*,  r«2/(37-2),  9*(r+l)/(r+2) 
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Here,  //(jr)is  the  Heaviside  function  and  !^(x)  is  the  modified  Bessel  function.  Solutions 

for  the  other  components  may  be  expressed  similarly.  The  expression  for  the  velocity  at 
the  lower  boundary  is 

"'4.  .  (3.23) 


where  the  factor  D  is 

Z)=(y«p,)*'fl'-V+2)'-*r(l-?)r(2^-l)r'(^)sin(2«g. 


At  the  outer  edge  of  this  region,  the  nearly  uniform  medium  approximation  becomes 
applicable.  Accordingly,  disturbances  will  propagate  along  characteristics.  The 
disturbances  will  enter  the  main  part  of  the  shock  layer  along  positive-slope  first  family 
characteristics  and  propagate  up  to  the  shock  wave,  lliey  then  reflect  a.s  characteristics  of 
the  opposite  family  toward  the  boundary  layer.  From  an  estimate  of  a  typical  slope  of  the 
characteristics  within  the  shock  layer  it  can  be  concluded  that  the  reflection  zone  lies 
downstream  of  a  distance  of  the  order  of  L .  Assuming  the  upstream  influence  through 
the  boundary  layer  is  weak,  the  reflection  effect  can  therefore  be  neglected. 

Another  comment  concerns  matching  with  the  solution  in  the  main  part  of  the  boundary 
layer.  To  provide  this  matching  Bush  in  [2]  has  demonstrated  that  a  thin  transitional  layer 
needs  to  interspersed  between  the  boundary  and  shock  layers  for  the  case  /  *  1  where 
the  singularity  of  the  mean  flow  is  being  smoothed  out.  Without  special  analysis,  it  is 
evident  that  the  viscous  forces  are  still  negligible  in  this  region  (for  the  perturbations 
only)  and  leading-order  equations  for  the  vertical  component  of  the  velocity  and  pressure 
are  very  simple,  since  both  variables  are  constant  in  the  vertical  direction.  This  behavior 
may  be  established  by  direct  ordcr-of-magnitude  e.stimations  and  also  follow.s  from  the 
fact  that  the  pressure  and  vertical  velocity  are  finite  at  both  edges  of  the  transition  region 
while  its  thickness  is  much  less  than  the  thicknesses  of  the  boundary  und  upper  layers.  In 
other  words,  the  transitional  layer  is  passive  with  respect  to  type  of  the  perturbations 
considered. 

On  the  other  hand,  the  typical  amplitude  of  the  vertical  velocity  in  the  boundary  layer  is 
much  greater  than  within  the  upper  deck.  Therefore,  the  necessary  matching  condition  is 
that  the  boundary  layer  velocity  is  zero  at  the  upper  edge,  Le. 

v„aOaty,sy3o.  (3.24) 

This  means  that  the  boundary  layer  disturbance.s  are  almost  fully  reflected  by  the  shock 
layer.  This  result  is  actually  the  same  as  that  obtained  in  [15]  for  rapidly  varying  steady 
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perturbations.  Applying  the  requirement  (3.24)  to  the  formula  (3.15)  a  relation  between 
the  pressure  and  displacement  thickness 

-a^A/aXi,  +C,(yiOo)JJ"{(?^  Pi\ldx\  +  r«j'^4p„)dy3  =  0  (3.25) 

is  obtained.  Equation  (3.25)  completes  the  problem  (3.3)  for  the  viscous  sublayer.  This 
gives  a  new  kind  of  triple-deck  problem  with  a  different  interaction  law  than  previously 
considered.  In  the  two^imensional  case,  the  last  term  in  (3.25)  degenerates,  and  the 
formula  is  reduced  to  a  more  simple  form 

A+constpj,  =  0 

which  was  analyzed  in[  16]  and  [21]. 

4.3.5  Final  formulation. 

The  problem  (3.3)  with  the  interaction  law  (3.25)  is  now  considered.  To  reduce  it  to  an 
invariant  form,  an  additional  transformation  of  the  dependent  and  independent  variables 

(xpz,) *PoC,'’yi'(X.Z);  (a„.w„)  = 

~ PqCI  Yt  v„  =«  p„  efyjoV 

Pu  ~  F ,  (3.26) 

is  applied.  A  new  parameter  6  depending  on  the  properties  of  the  gas  via  parameters 
(o.y.a) 

(3-27) 

is  also  introduced.  Numerical  values  for  all  factors  used  in  this  scaling  are  given  for  air 
and  helium  in  Table  2  of  the  Appendix  .  The  right-hand  side  of  the  last  expression  is  the 
average  vnlue  of  the  inverse  square  of  the  local  Mach  number  across  the  toundary  layer 
and  is  less  than  unity  for  both  these  gases  (probably  for  most  others  as  well).  After 
substitution  of  the  rescaled  variables  into  (3.3),  (3.25),  the  governing  system  of  equations 
takes  (he  form 

Ur  +  UU^  +  VUy  +  WUy  +Px=  U„ 

WrfUW„+VW,i-  WWj  *Pz‘Wrr 


(3.28a) 

(3.28b) 

(3.28c) 


and  (he  boundary  conditions  are 
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U  =  V  =  W  =  0sitY=O 

(3.29a) 

U=Y+A(X,Z,T)+...,W^OasY^«> 

(3.29b) 

Aa  +sin*  eP„  -  cos*  ePix  =0. 

(3.29c) 

Here,  subscripts  denote  corresponding  partial  derivatives  and  the  function  A  is 
renormalized  in  the  same  way  as  the  streamwisc  velocity.  That  formulation  docs  not 
include  any  effects  of  forced  disturbances  coming  from  the  surface  or  external  flow.  For  a 
perturbed  surface,  such  generalization  can  be  made  by  simply  changing  the  boundary 
condition  (3.29a),  to 

U‘W^O,V‘FratY  =  F{X,Z,  T) ,  (3.29a’) 

where  the  function  F  describes  the  scaled  shape  of  the  surface.  Reformulation  for  the 
case  of  external  flow  disturbance.s  demands  a  special  study  of  their  transformation  within 
the  shock  layer  and  is  outside  the  scope  of  this  chapter.  Instead  of  the  homogeneous  wall 
conditions  (3.29a)  the  modified  version  (3.29a')  will  be  applied  in  what  follows.  In 
keeping  with  the  primary  thru.st  of  this  study,  a  linear  stability  analysis  of  the  resulting 
problem  will  be  made. 

4.4  UnearizMthn  mnd  solution 

In  the  linear  approximation,  the  solution  is  expanded  in  a  power  series  of  a  small 
parameter  (a)  characterizing  the  amplitude  of  the  perturbation  so  that 

£/  =  y+ au+...;lV  =  aw+...;V- av+...;P  =0/7+...;^=  a/+...  (4.1) 

Here  and  below,  small  letters  denote  only  scaled  values  and  do  not  correspond  to  primary 
dimensional  variables  as  in  Sections  2,3.  Substitution  of  this  representation  into 


(3.28).(3.29)  gives 

+  V + Pj, »  u„  (4.2a) 

WT  +  Ywg-*-pj[  =  Wyr  (4,2b) 

(4.2c) 

a»-/,w=0,vx./yat  y»0  (4.2d) 

+  •'*»*  ^Pa  “  cos*  dpg  »  0,^  -♦  0  at  y  (4.2e) 
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For  the  solution  of  this  problem  to  be  unique,  an  initial  condition  at  some  instant  of  time 
must  be  specified.  This  condition  is  assumed  to  be 

«=v®>v=p  =  /=0at  r<0,  (4.3) 

implying  that  surface  variations  are  switched  on  at  r=0.  These  variations  are  assumed 
to  be  bounded  continuous  functions  vanishing  in  the  X-and  Z-  directions 

fiX.Z,T)  6  C./  -»  0  as  X*  +Z*  (4.4) 

From  the  restriction  (4.4)  it  is  plausible  that  the  solution  is  also  bounded  at  any  finite 
instant  of  time*  and  the  double  Fourier  transformation  with  respect  to  the  streamwise  and 
spanwise  variables  can  be  applied.  Based  on  tf..  standard  assumption  that  the  temporal 
growth  is  not  faster  than  the  exponential  one,  the  Laplace  transformation  with  respect  to 
the  time  can  be  applied.  These  transformations  written  symbolically  as 

UT -4  -i(o)F{Z  ip,X  -* ta)[u,w,v,p,f),ia,P)  e  (-«>o.+e«>).(-f6)) >  0  (4.5) 

reduce  the  problem  (4.2), (4.3)  to  a  system  of  ordinary  differential  equations  (the  Laplace 
transform  variable  was  rotated  in  the  complex  plane  for  a  more  obvious  analogy  with  the 
canonical  triple^leck  formulation  to  connect  with  previous  work).  The  solution  procedure 
has  been  used  by  many  authors  (for  the  three-dimensional  unsteady  case  see  for  example, 
[12,22].  Differences  with  the  present  formulation  relate  only  to  the  form  of  the 
interaction  law  used.  For  this  reason  details  of  all  intermediate  steps  including  the 
inversions  back  to  the  physical  space  will  not  be  presented.  Only  the  integral 
representation  for  pressure  is  given  which  is 


p  =  (2^)"' J  Jcxp{iaX +ipZ]dadpj^exp{-iQ)T]D'\a,fi,o))f{a,P,(o)d(ij  (4.6) 

where  the  function  f(a,p,0D)  is  the  shape  of  the  surface  under  the  transformations  (4.5). 
For  the  problem  considered  here,  the  dispersion  function  D  expressed  in  terms  of  the 
Airy  function,  Ai  is 

Aa,^,®)=  (ia)''’n+(^ /a)*](A/(-Zo))‘'J’ Ai(z)*+sin* 

03/ a)*  cos*  d  (4.7) 

Zfsioi/  (/a)*'*;-/r  <arg(/a)<  «,(»a)''’  >  0  for  la  >  0. 


*  Of  eoune.  it  does  not  exclude  (he  pOMibite  existence  of  unbounded  (unsuble)  xolutions  in  (he  lar  jc-time 
limit. 
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The  contour  C  in  (4.6)  passes  from  left  to  right  half  planes  of  the  complex  to -plane 
above  all  singularities  of  the  integrand  (due  to  causality  principle),  and  both  its  ends  are 
directed  into  the  lower  halfplanc  to  provide  convergence  of  the  inversion  integral  for  p . 
Since  the  integrand  is  a  function  of  the  wavenumber  (a,/3),  the  contour  C  also  depends 
on  the  wavenumber. 

Properties  of  the  solution  are  determined  by  zeroth  of  the  dispersion  relation  (4.7)  which 
correspond  to  eigenmodes  of  our  problem. 

4.4.1  Analysis  of  the  eigenmodes 

Since  the  square  of  the  spanwise  wavenumber  p  occurs  in  the  di.spcrsion  relation  only 
the  case  need  be  considered.  Results  of  numerical  analysis  of  the  dispersion 

relation  for  6=61.67_  corresponding  to  air  (see  the  Appendix)  arc  presented  in  Fig.  3. 
The  simplest  picture  is  observed  in  the  two-dimensional  case  There  is  a 

denumerable  set  of  eigenwavenumbers  a  with  real  positive  frequencies,  one  of  which 
lies  in  the  left  lower  quadrant  of  the  complex  a -plane,  and  others  move  from  the  zero  at 
oj  =  0  into  the  right  upper  quadrant.  In  the  steady  case,  the  first  root  may  be  written 
analytically  as 

ia  =  {-3A/(0)sin*6}*.  (4.8) 

Accounting  for  the  scalings  (2.3),  (3.26).  this  value  agrees  in  the  Newtonian  limit 
discussed  in  [23]  (noting  that  according  to  (3.27)  $in6  -»  1  as  y  1 )  with  the  exponent 

in  [16]  their  power  law  for  disturbance  propagation  is  converted  to  the  approximate 
exponential  form 


Fig.  3  Movement  of  the  roots  in  complex  plane  as  frequency  goes  from  zero  to  infinity. 
(Solid,  three-dimensional  ca.se  at  ^  =  1  -.dashed,  plane  modes  at  ^  =  0).  Arrows 
indicate  increasing  frequency. 

(x t LY  -  txplnAxt  L)asn^m,n^x  /  Lm 0(1)  (4.9) 
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This  implies  that  the  two-dimensional  steady  state  version  of  the  asymptotic  regime 
considered  here  describes  the  Neiland  •  Brown  and  Stewartson  eigensolution. 

In  the  three-dimensional  case,  all  these  modes  also  exist  and  exhibit  the  same  qualitative 
behavior,  but  in  addition,  two  more  eigenwaves  appear.  One  of  these  is  located  in  the  left 
upper  quadrant  and,  like  all  previous  ones,  does  not  intersect  the  real  axis.  Another  one 
starts  from  the  zero  in  the  right  upper  quadrant,  but  after  some  critical  value  of  the 
frequency,  intersects  the  real  axis.  The  large-frequency  limits  of  these  modes  are 
respectively 

a“’  =  -^cot0  +  (O;a‘^’ sPeotS-ZOaso-^oo.  (4.10) 

The  imaginary  part  of  the  second  three-dimensional  mode  versus  the  frequency  is  shown 
in  Fig.  4  at  different  positive  spanwise  wavenumbers  j3 .  As  is  increased,  the  depth  of 
the  minimum  also  increases,  i.e. 

min  {Im  ( a  )}->-*»  as  -»+»,  (4.11a) 

or,  similarly  for  the  real  wavenumbers  and  complex  frequency 

max„,,{Im(o)))  +««as/J -» +«.  {4.11b) 

Equations  (4.1  la)  and  (4.11b)  imply  that  the  contour  C  in  the  integral  (4.6)  actually 
depends  only  on  the  spanwise  wavenumber  ,  and  hence  the  order  of  integration  with 
respect  to  frequency  and  streamwise  wavenumber  can  be  changed.  Accordingly,  Briggs’ 
[24]  criterion  can  be  used  to  specity  the  direction  of  propagation  of  each  eigenmode. 
According  to  this  criterion,  the  streamwise  wavenumbw  of  the  upstream  (downstream) 
propagating  eigenmode  should  belong  to  the  lower  (upper)  halfplane  for  a  sufficiently 
large  imaginary  part  of  the  frequency.  Specially  conductetl  computations  show  that  the 
root  lying  entirely  in  the  lower  half  plane  at  real  frequencies  correspond  to  the  upstream 
propagating  wave,  and  all  other  modes  propagate  downstream  with  one  of  them  being 
unstable.  In  particular,  the  stationary  mode  of  Neiland-Brown  and  Stewartson  propagates 
in  the  upstream  direction. 

Another  important  consequence  of  (4.11)  is  that  the  growth  rate  of  the  unstable 
eigenmode  goes  to  infinity  as  the  wavelength  decreases.  In  fact,  realistic  physical  systems 
do  not  contain  infinitely  fast  growing  modes,  and  the  present  situation  is  a  result  of 
simplifications  made  in  the  asymptotic  limit  considered.  Thus,  the  possibility  of  well- 
posing  of  the  forced  problem  requires  additional  investigation. 

4.4.2.  Amlysis  of  the  initial  botmdary-value  problem 

Li  typical  formulations  of  the  forced  problem,  the  large  time  behavior  of  the  boundary 
surface  is  a.ssumed  to  be  harmonic.  This  implies  that  the  function  f(a,p,(o)  is  analytic 
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with  a  simple  pole  at  some  real  frequency  6)  =  >  0 .  To  provide  the  existence  of  the 

solution  of  the  forced  problem  which  is  equivalent  to  convergence  of  the  integral  (4.6), 
the  class  of  the  surface  .shape  functions  is  restricted  by  imposing  limitations  on  their 
decay  in  wavenumber  space.  Properties  of  the  solution  for  this  set  are  investigated  in 
what  follows.  Implications  of  relaxing  these  limitations  are  then  considered.  For 
simplicity,  the  special  case  of  shape  functions  of  compact  support  in  the  h*direction 

0 for l/3l> /J,  “♦+*» as (4.12) 

is  studied.  Within  this  clasi  of  functions,  the  infinitely  fast  growing  waves  are  excluded 
and  the  solution  doe.s  exist.  The  region  far  downstream  from  the  source  of  the 
perturbations  where  the  unstable  eigenmode  dominates  is  now  considered.  Accordingly, 
only  the  contribution  of  the  unstable  eigenmode  is  important.  For  convenience,  the  order 
of  integration  is  reversed  with  respect  to  or  undo),  and  the  integrals  with  respect  to  each 
variable  are  considered  successively.  The  first  integral  may  be  estimated  as 

[‘y‘^f(a,p,a))D~\a,p,o})da*» 

27[iexl^iia^X)f(a^,p,Cl>){^,ia^,p,(o)y'  (4.13) 

Here  or,  =ar,(d,o))  is  the  root  of  the  dispersion  relation  corresponding  to  the  unstable 
mode.  The  next  integral  with  respect  to  the  frequency  contains  contributions  of  the  shape 
function  pole  at  (U=  tUoand  the  branch  points  dD/  da^O.  Dividing  it  into  the  time- 
harmonic  contribution  written  in  an  explicit  form  and  a  remainder  term  /, ,  gives 

J^exp(-iC!)T+ ior,X)/(a,./8,a))|D„(flr|,^,<u))''d(u  » 
-2nfexp(-toor+ior,o^resl/(or,o,/J,tUo)I{D„(a,o,/?,o)o}’'  +  4.  ('^•I^) 

where  ar,o  =  Ci^(^,o>o).  The  first  term  of  the  representation  (4.14)  includes  only  the 
modes  with  bounded  growth  rates  so  that  its  behavior  is  regular  us  j3  --» <«.  The  second 
term  contains  the  modes  growing  infinitely  fast  as  /J  ->  <» . 

At  finite  X  and  T,  /,  in  (4.14)  can  be  estimated  approximately  by  the  method  of 
steepest  de.scent  assuming  the  spanwisc  wavenumber  b  is  large.  In  that  limit  the 
eigenvalue  is  expanded  into  the  descending  power  series  in  /3 

<u  =  ^;  a,  =  ^Ao(fl)+/?"*A,+..., 

(4.15) 

where  the  coefficients  A*  are  described  by  the  implicit  equations  obtained  after 
substitution  of  the  expansion  (4.15)  into  the  dispersion  relation  (4.7) 
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(1  +  Ag,  )Ao  +  flCsin  ^  OA^,  -  cos  *  6)  =  0 ,  (4.1 6a) 

A,{(A^  - 1)  /  ft  +  2cos'  e  /  Ao)-Al(-iCl)-^''‘(Al  + 1)  =  0.  (4. 16b) 

There  are  three  branches  of  (4.16a).  The  branch  of  interest  here  is  confined  within  the 
finite  region  (<^  (4.10)) 

-cot0<Aj<cot0,IA,l<«>  in-»<Q<oo,  (4.17) 

The  location  of  the  saddle  point  in  the  ©-plane  i.s  determined  from 
d(a^X-0}T)JdQ)  =  0. 

This  gives 

(4.18) 

where  the  expansion  terms  are  obtained  from  conditions 

dAJda=TlX,a,  =-dAJdSl[d^AJd^iy ,  (4.19) 

with  all  derivatives  taken  at  (2=  Elcmentaiy  analysis  of  the  conditions  (4.19)  with 
(4.16)  shows  that  the  saddle  point  exists  at  TIX<co^9,  and  its  contribution  is 
exponentially  large,  i.e., 

=  0(expCC^''^ D.  as  /J  -»  oo;X  >  r  /  cos’  0,  C  >  0.  (4.20) 

In  the  region  X<T / cos’ 0,  the  saddle  point  contribution  is  absent  at  large  b  and  the 
remainder  part  /,  of  the  integral  defined  by  (4.14)  should  be  evaluated  at  finite  values  of 
P,  For  Xand  T  large,  /,  is  evaluated  from  the  contribution  of  the  branch  points 
(3D/  da  =  0)  which  arc  closest  to  the  real  axes  in  a-  and  ©-planes.  If  position  of  these 
points  is  denoted  as 


oc^(P).a>„{p) 


the  value  of  /,  is  roughly  expressed  by 

/,  =0{cxp[lm(©*,(;3)r-a*,(/3)X]).  (4.21) 

This  contribution  represents  a  vorlicity  spot  appearing  when  the  oscillator  starts  its 
motion.  This  effect  was  considered  in  detail  in  [25]  for  two-dimensional  subsonic 
boundary  layers  and  in  [22]  for  the  three-dimensional  case.  All  branches  occur  in  the 
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lower  half  of  the  a-  and  (o -planes  and  do  not  intersect  the  real  axis,  so  that  the  vorticity 
spot  is  convected  downstream  with  Hnite  velocity:  The  present  case  differs  mainly 
because  some  of  the  branch  points  approach  zero  from  below  as  jS  -» 0.  However,  they 
do  not  intersect  the  real  axes  so  that 

ImtWi,  S0,Ima4,  <0,ImQ)3,/Ima*,^c„,.,  >0  (4.22) 

Behavior  of  and  versus  fi  for  the  first  branching  point  of  this 

kind  is  shown  in  Fig.  5.  Relations  (4.22). imply  that  /, represents  disturbances  moving 
downstream  with  (he  finite  velocity.  On  a  ray  X IT  =c<\lcos^ 0,  (Case(i))  the 
disturbances  are  growing  not  faster  than  exponentially  aind'are  bounded  at  any  P,  while 
for  X/T>1  /cos^  0,  (Case  (ii)),  they  arc  unbounded  as  ^  in  agreement  with  the 
formula  (4.20). 

The  last  integration  in  (4.6)  with  respect  to  the  spanwise  wavenumber  gives  a  finite  limit 
as  ->oo  in  Case(i)  and  generally  has  no  limit  in  Casc(ii)*.  Therefore,  the  asymptotic 
regime  considered  is  applicable  for  the  region  Jf/r<  l/cos*fl  only.  The  perturbation 
field  further  downstream  requires  a  special  analysis  of  the  high  frequency  eigenmodes. 

For  increasing  time,  the  vorticity  spot  displays  far  downstream  and  the  perturbations  field 
consists  of  time-harmonic  oscillations  only.  Thus 

P  •  J7_exp(i)3Z+ia,oX-i®or]  res  (/(o, )](/)„(«, o,)3,£Oo))‘‘d^ 
Z~»+«,X/r<c;^  (4.23) 

The  distribution  of  the  Imaginary  part  of  0,^  versus  p  is  shown  in  Fig.  6  for  different 
frequencies.  When  the  frequency  is  small,  there  is  a  minimum  of  Im  at  finite  P .  If  the 
frequency  exceeds  some  critical  value  (O*  (»0.30  for  air,  -0.78  for  helium),  the 
minimum  disappears  and  Imo,o  tends  monotonically  to  its  lower  boundary  as/3  -» 
Accordingly,  the  solution  behavior  is  different  for  both  cases. 

Estimating  the  integral  (4.23)  by  the  method  of  steepest  descent  on  a  ray  Z  =  X  tan  ,  the 
position  of  the  saddle  point  is  determined  by  a  condition 

d(a^o+P^an(f>)^dP=^0,  (4.24a) 

which  is  equivalent  to  Values  of  cc,,P^  corre.sponding  to  the  saddle  point  can  be 

calculated  from  (4.24b),  and  the  growth  rate  cr  can  be  determined  along  the  ray 
considered  to  be 


*  Even  If  (he  shape  function  is  rapidly  decaying  in  P  -direction  so  that  (he  integral  exists  as  /3„  ,  the 

temporal  growth  of  the  solution  is  fnaterthan  exponential  and  the  application  of  the  Laplace  transformation 
is  not  correct 
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Fig.6  Imcr  vs.  j3  for  unstable  mode  at  different  frequencies ;  £» =.05, ,  1 5,  .30. 

/<//)  + tan  ip  =  0.  (4.24b) 

If  there  is  an  angle  (p^  of  maximum  amplification,  the  relation 

da/  d<p  =  hn[dP^  /  d<^da^^  /dfi-i-  tan  <p] + /  cos*  9}  =  0  at  (p  =  <p^.  (4.26) 

holds.  Combination  of  (4.26)  with  (4.24b)  leads  to  tile  conclusion  that  on  the  ray  of 
maximum  amplification  the  saddle  point  coincides  with  minimum  of  the  function 
Ima,Q(^)  at  real  /3.  Results  of  calculations  of  the  angle  9^  and  the  corresponding 
growth  rates  are  given  in  Fig.  7.  For  frequencies  below  the  critical  value  O)  *,  this  angle  is 
different  from  zero  so  that  the  far  field  disturbance  distribution  consists  of  two  Gaussian 

humps  in  the  spanwise  direction  at  Z  =  ±Xtan9,  with  typical  width  0(VT).  Their 
maximum  amplitude  is 

/>  -  (ff  /  a,o*()3^  )3ff  *  exp(»a^3f + -  imoT]  X 

(4.27) 

The  case  OJ^XO*  cannot  be  considered  in  the  same  way  as  previously,  and  requires  a 
special  analysis.  For  this  case,  the  main  contribution  in  the  integral  (4.23)  is  given  by  the 
most  unstable  modes  now  concentrated  at  the  large  p.  The  approximate  form  of  the 
eigenvalue  in  this  region  is 
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aio  =  Oo((Oo)*P'\+-  (4.28) 

where  the  terms  of  the  expansion  are  obtained  from 

|.Ai(z)ife  -cos^  0. 2o  =  ‘^0  ^ 

a,=-a^JDM)-  (4-29) 

Dependence  on  the  frequency  is  shown  in  Fig.  4a  as  the  limiting  curve  as  )3  -> «.  At 
the  low  frequencies  (the  neutral  frequency  for  air  o,  =0.1 16,  for  helium  ru,  =0.303),  such 
limiting  modes  are  stable,  while  at  high  frequencies  the  growth  rate  is  unbounded.  Due  to 
the  absence  of  the  minimum  in  the  distribution  Ima,a(/9)  the  uncquality 

lmai>0 

holds.  For  this  condition,  the  contribution  of  the  region  Iina,  >0.in  the  integral  (4.23)  is 
exponentially  small,  and  the  approximate  representation  (4.28)  can  be  used  in  the  whole 
region  of  integration.  Taking  for  simplicity  the  residual  of  the  shape  function  equal  to  one 
(^-function  shape  oscillator  in  the  physical  space),  the  solution  can  be  written  in  the  form 

p  -  i(0j^J{2,X)  t [D^ia^)  /aj)' 

JiZ,X)  =  +  UhP'^X]dp  /  .  r4.30) 

Here 

D  =  0 / Dq+...  as p  oo.  The  value  of  J  at  Z=0is  finite 

J  =  {-nlia,xr\ 

and  as  ZjX  it  decays  exponentially  for  Imo,  >0.Thus 
J  -  (JT /  6a,X)''*  exp(3tZ’'*X’'*(fl,  /  4)"’ ] . 


Therefore,  at  supercritical  frequencies  0)  >  m*,  the  perturbation  is  concentrated  in  a  thin 
region  near  the  x-axis  with  atypical  width  0{\/-/x).  A  qualiutivc  picture  of  the  field 
of  perturbations  is  shown  in  Fig.  8. 4.S 

4.5  OlMCUBslon  end  concluBlone 

The  foregoing  has  shown  that  the  strong  interaction  hypersonic  boundary  layer  subject  to 
surface  cooling  of  the  surface  is  uastablc  with  respect  to  low  pha.se  speed  (low  frequency) 
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SC.I8»E.0»2M 


Fig.  7  Character! itics  of  the  far  downstream  field  from  vibrator,  (a)  Angle  of  ray  of 
maximum  amplification  <p„’ys.  frequency  of  vibrator,  (b)  Growth  rate  on  the  ray 
of  the  maximum  amplification  versus  frequency.  (i):Air,(ii):Heiium. 
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Fig.  8  Schematic  picture  of  perturbation  field  induced  by  time-harmonic  source  on  a 
surface;. 

viscou.s  cigenwaves.  This  does  necessarily  preclude  other  types  of  instabilities.  On  the 
contrary,  in  addition  to  inviscid  inflectional  instability  of  the  transitional  layer  considered 
before  in  [  1 1  ],  another  inflectional  instability  should  occur,  because  for  sufficiently  strong 
cooling  a  generalized  inflection  point  appears  within  the  boundary  layer.  (This  can  easily 
be  concluded  from  analysis  of  the  mean  flow  asymptotics  near  the  wall).  Results  of 
computations  for  air  presented  in  Fig.  9b  (see  the  Appendix)  indicate  that  such  a  point  is 
located  approximately  in  the  middle  of  the  boundary  layer.  Moreover,  as  follows  from  the 
analysis  of  [2],  even  for  the  hot  wall  and  realistic  gas  parameten  {a,  y,OJ)  the  mean  flow 
velocity  tends  to  the  freesta'am  value  from  above,  and  an  additional  inflection  point 
exists  near  the  upper  edge  of  the  boundary  layer.  Our  computations  (sec  the  note  in  the 
end  of  the  Appendix)  show  that  this  point  is  extremely  close  to  the  boundary  layer  edge 
and  must  be  obscured  by  the  transitional  layer  for  realistic  Mach  numbers. 

In  any  case,  the.se  inflectional  instabilities  have  a  typical  frequency  not  less  than  the  ratio 
of  the  frecstrearn  velocity  to  the  thickness  of  the  iHJundary  layer  in  the  strong  interaction 
region.  Typical  frequerKics  associated  with  transition  to  tuibuloKC  are  based  on  the 
boundary  layer  thickness  far  downstream  and  are  much  smaller  than  these  inflectional 
instability  frequencies.  Therefore,  the  viscous  instability  considered  may  have  more 
practical  interest  even  if  its  growth  rate  is  smaller  than  the  inviscid  one.  This  conjecture  is 
supported  by  the  experimental  study  by  Holden  [18,26]  of  the  hypcisonic  boundary  layer 
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on  the  cooled  cone  at  zero  angle  of  attack.  He  observed  the  occurrence  of  highly  vi.sible 
waves  in  the  region  where  the  ratio  of  thickness  of  the  shock  and  boundary  layers  was 


Fig.  9  (a)  Mean  flow  distribution  for  velocity  and  tcmperature.(i) 

(ii)r,o(yj)- 10  (b)  Distribution  of  «Jo(y3)/^3o^)'3)''''^hin  boundary  layer. 
Minimum  corresponds  to  generalized  inflection  point. 
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»  10  (under  strong  interaction  it  is  -  1.7  for  flat  plate  in  air).  These  waves  were  thrcc- 
dimeasional,  had  a  continuously  growing  phase  speed  from  0.5  at  the  start  of  observations 
up  to  0.85  at  the  end  of  transition.  The  effect  of  three-dimensionality  is  believed  to  be 
related  to  the  type  of  instability  considered  here. 

Another  interesting  ob.ser\'ation  was  made  in  [19]  in  the  study  of  hypersonic  transition 
induced  by  a  solitary  roughness.  Therein,  it  was  found  that  if  the  distance  between  the 
nose  and  the  roughne.ss  is  less  than  some  critical  value,  the  wake  behind  the  roughness 
has  a  thread-like  form,  and  bursting  into  developed  transition  occurred  far  downstream. 
Our  analysis  shows  that  if  the  frequency  of  the  source  is  less  than  some  critical  value,  the 
far  downstream  perturbation  field  is  like  a  swallow  tail.  For  the  high-frequency  case,  it 
narrows,  with  the  width  of  the  tail  being  inversely  proportional  to  the  square  root  of  the 
distance  from  the  source.  The  downstream  growth  rate  increases  with  frequency  (cf. 
results  of  calculations  in  Fig.  7  and  qualitative  picture  in  Fig.  8).  Therefore  the  high 
frequency  disturbances  must  dominate.  Accordingly,  the  wake  behind  a  point  source  with 
a  broad  spectrum  of  frequencies  must  concentrate  near  the  central  axis  of  the  source,  in 
agreement  with  the  previously  mentioned  experimental  result. 

The  physical  mechanism  of  instability  obtained  is  close  to  that  analyzed  in  [12]  for 
moderate  supersonic  speeds.  In  this  case,  instability  occurs  for  the  waves  directed  outside 
the  free-stream  Mach  cone  which  arc  therefore  effectively  subsonic  and  qualitatively  obey 
the  same  mechanism  as  Tollmien-Schlichting  waves  in  incompressible  boundary  layers. 
For  the  case  considered  here,  the  disturbances  are  confined  by  the  boundary  layer  and  do 
not  penetrate  outside  it,  so  that  the  interaction  is  realized  between  the  vi.scou.s  .sublayer 
and  the  main  part  of  the  boundary  layer.  Hence,  the  role  of  the  freestream  Mach  number 
is  replaced  by  some  effective  Mach  number  averaged  across  the  boundary  layer.  Since  the 
local  boundary  layer  Mach  number  is  of  the  order  of  one  the  Mach  cone  angle  takes  on 
finite  values  different  from  n  tl  and  indeptndent  of  (he  free-stream  Mach  number. 

The  important  difference  from  the  moderate  supersonic  regime  is  the  presence  of  the 
infinitely  fast  growing  modes  in  the  short-wave  limit.  Therefore,  the  forced  problem  is 
well-posed  on  the  class  of  shape  functions  of  compact  support  in  Fourier  space.  For  a 
timc-harmonic  source  switching  on  at  some  instant  of  time,  the  fast  growing  modes 
contribute  to  a  vorticity  spot  which  originates  at  the  start  of  the  oscillations  and  is 
convected  downstream  with  finite  velocity.  Continuation  to  a  wider  class  of  the  shape 
functions  can  be  made  in  the  region  of  physical  space  pas.scd  by  the  high-frequency 
vorticity  spot  and  containing  the  time  harmonic  component  only.  The  high-frequency  part 
of  the  vorticity  spot  cannot  be  described  within  the  present  framework  but  needs 
con.sidcraiion  in  another  asymptotic  regime. 
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4.7  Appendix.  Undisturbed  boundary  layer. 

In  terms  of  the  physical  variables  the  system  of  equations  controlling  the  mean  boundary 
layer  flow  is 

=PiaOyiUio-^Vy^)  (A.l) 

=  -Po '2y=  ((7'3,)Xr 

-('F/4)r„+PoW,o(y-i)/2r=(r-iW“;o)'+<^'’(21S7';)' 

Pno^jo  -Po‘ 

Here,  'F  is  a  stream  function  introduced  for  more  convenient  numerical  integration,  and 
prime  denotes  differentiation  with  respect  to  .  Boundary  conditions  are 

'F  =  t4^  =  r„  =  0  aty,  =  0  {A.2) 

'F  ->«,Vjo=3y„/4,tt,o  =  l,TM  =  0at  y, 

with  the  location  of  the  upper  boundary  y^  determined  simultaneously  with  the  solution. 
The  pressure  Pg  is  not  a  prescribed  value  but  is  determined  only  after  coupling  the 
problem  (A.l), (A.2)  with  the  problem  for  the  inviscid  shock  layer.  This  layer  is  confined 
between  the  upper  edge  of  the  boundary  layer  y^and  the  position  of  the  shock  wave 
Due  to  self-similarity,  its  profile  is  also  described  by  a  system  of  ordinary 
differential  equations  which  are  the  reduced  classical  ones  of  hypersonic  small 
perturbation  theory,  and  may  be  treated  as  a  locally  uniform.  Boundary  conditions  behind 
the  shock  wave  are  the  Rankinc-Hugoniot  relations  and  the  main-order  matching 
conditions  with  the  boundary  layer  which  are  continuity  of  the  pressure  and  vertical 
velocity. 

For  complete  formulation  of  the  boundary  layer  problem  the  detailed  structure  of  the. 
shock  layer  flow  is  not  necessary.  Only  an  additional  relation  between  the  pressure  and 
the  vertical  velocity  at  the  upper  edge  (the  global  interaction  law)  is  sufficient  for 
obtaining  a  unique  boundary  layer  solution.  For  evaluation  of  the  mean  flow  in  the 
Region  4  lying  at  the  bottom  of  the  shock  layer,  the  location  of  the  shock  wave  is 
required.  Numerical  data  appropriate  for  the.se  both  purposes  were  obtained  in  [ISjfor 
mono-  and  diatomic  ga.ses.  Their  .scaling  formulae  (2.1)  and  .solution  representation 
(2.9),(2. 10)  gives 

Fo{vM(y»))'*  =  (i/2)r(r+i)Fo(0){Vo(0)r*.  (a.3) 
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;'3,/>’3o  =  a'2xy+i){v,(0)r'. 

where  ?(,(0),  ^5(0)  are  their  scaled  pressure  and  vertical  velocity  components  at  the 
bottom  of  the  shock  layer.  Substitution  of  their  numerical  values  gives  the  lesulLs  in  Table 
1.  For  comparison  the  analytical  result  for  Newtonian  approximation  (/->  1)  is  also 
shown. 


Table  1 


molecular  stnicturc 

y 

bssbhhhi 

^31  ^  Yv) 

one-atom 

5/3 

3.098 

2.086 

two- atom 

7/5 

1.990 

1.691 

1+0 

2/3 

1 

Noticeable  differences  in  the  global  interaction  law  for  the  different  values  of  the  specific 
heats  ratio  appear  because  our  primary  scaling  does  not  involve  this  parameter. 

For  integration  of  the  problem  (A.l),  the  boundary  conditions  (A. 2)  must  be  replaced  by 
their  asymptotic  counterparts.  Conditions  at  the  wall  follow  from  matching  with  the 
solution  (3.7)  within  the  buffer  layer 


=  (C(l  +  tw)yj 


(A.4) 


Uy,={CJC,){C,{U(D)y, 

Asymptotics  near  the  upper  edge  were  obtained  f21  using  stream  function  as  the 
independent  variable.  In  physical  variables,  they  take  the  form 

7'„  =  Constj(yj„-y3)"‘'**''  (A.5) 

i/„  =  l  +  Const,(y3„-y, 

where 

Const,  ={p,<T(i  +  £U)(I-  +(y  -  1X1  -  <y)y'']/  4}’'"*"’  (A.6) 

Const,  =(1  -  <u){y  +  (y  -  I)(l  -  £u)l-'(l- (T(y  -  l)(l  -  £u)  /  yr'Const,. 
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Calculations  were  conducted  by  the  Runge-Kutta  method  with  variable  step  of  integration 
for  two  sets  of  the  gas  parameters  a.y.o  corresponding  to  helium  and  air.  Numerical 
values  characterizing  the  mean  flow  profile  are  presented  in  Table  2 


Table  2 


Gas 

a 

y 

(O 

Cl 

C2 

Fo 

>'jo 

0“ 

He 

2/3 

5/3 

.647 

.4294 

.1334 

.4538 

.5103 

1.064 

52.94 

Air 

.72 

7/5 

3/4 

.3895 

.1328 

.2151 

.4384 

.7413 

61.67 

Actually,  these  data  are  not  far  from  the  ones  obtained  in  [3]  for  the  Chapman  viscosity 
law  and  Prandtl  number  equal  to  unity.  The  velocity  and  temperature  distributions  for  air 
are  presented  in  Fig.  9a,  and  distribution  of  the  value  m^(>’3)/Tjj(>’j)  characterizing  the 
location  of  the  generalized  inflection  point  is  plotted  in  Fig.  9b.  Behavior  near  the 


Fig.  10  Comparison  of  calculations  of  mean  flow  velocity  and  temperature  with 
a.symptotics  near  boundary  layer  edge,  (solid:  calculation;dashed:  asymptotics). 

(i):((Mjo.  ~l)10*;(ii)Tj<,10' vs. 
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boundary  layer  edge  and  comparison  with  the  asymptotics  (A.5)  arc  shown  in  Fig.  10  in  an 
enlarged  scale.  Note  that  in  agreement  with  the  [2]  the  velocity  approaches  unity  from 
above  (Const  i>0)  so  that  there  is  a  local  maximum  of  the  latter  within  the  boundary 
layer.  Its  position  is  extremely  close  to  the  upper  edge 

=107  10'’at  vjo-y,  •=6.96- lO"* 

and  is  not  resolved  on  the  large-scale  figure.  The  mean  flow  for  helium  exhibits  the  same 
qualitative  behavior. 
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5.  Unsteady  Newtonian  thin  shock  layers  and  hypersonic  flow 
stability 

S.1!ntroductlon 

AlUiough  the  stability  of  high  speed  flows  has  received  much  attention  in  the  recent 
literature,  major  complicating  aspects  have  not  been  treated  in  a  unified  way.  These 
features  include  the  combined  effects  of  the  finite  shock  displacement  on  the  boundary 
layer,  the  nonpai-allelism  of  the  flow  and  the  vorticity  introduced  by  the  shock  curvature. 
The  relevant  strucnire  of  the  shock  and  boundary  layers  has  been  treated  in  [1]-[9J.  In  [6] 
and  [7],  the  aforementioned  stability  issues  were  discussed  within  the  Hypersonic  Small 
Disturbance  approximation  for  the  inviscid  deck  strongly  interacting  with  the  hypersonic 
boundary  layer.  Equations  of  motion  for  the  mean  and  fluctuating  small  amplitude  flows 
were  analyzed.  Because  of  nonparallelism  in  this  framework,  the  spatial  part  of  the  waves 
cannot  be  treated  by  the  usual  Fourier  decomposition  and  an  initial  value  rather  than 
eigenproblem  for  spatial  stability  is  obtained.  The  initial  value  problem  leads  to  partial 
rather  than  ordinary  differential  equations  that  require  a  numerical  marching  method  for 
their  solution.  Results  indicate  that  the  specific  heat  ratio/  plays  a  major  role  in  the 
stability  of  flow  since  it  controls  the  reflection  of  waves  from  the  shock  and  the  radiation 
of  energy  in  the  shock  layer  whose  thickne.ss  scales  witli  /  - 1 . 

Early  experiments  such  as  those  described  in  [2]  showed  that  for  a  practically  interesting 
class  of  flows,  the  shock  layer  becomes  very  thin  compared  to  the  boundary  layer  near  the 
nose  of  hypersonic  flat  plates.  This  feature  and  the  desire  to  further  understand  the  shock 
and  boundary  layer  structure  encourages  the  use  of  the  Newtonian  approximation  y  1 . 
The  conneaion  with  flow  stability  motivates  the  study  of  this  approximation  in  an 
unsteady  context.  In  this  chapter,  limit  process  expansions  will  be  discus.sed  relevant  to 
unsteady  viscous  interactions  as  a  prelude  to  the  analysis  of  hypersonic  stability  and 
transition.  The  application  of  these  limits  is  an  unsteady  extension  of  the  steady  state 
analysis  of  (3].  Although  the  focus  here  is  the  treatment  of  viscous  interaction,  boundary 
layer  stability,  receptivity  and  transition,  the  results  derived  are  useful  in  inviscid 
hypersonic  unsteady  aerodynamic  methodology  and  load  prediction  as  well. 

S.2  AnatyalB 

Figure  1  schematically  indicates  strong  interaction  flow  near  the  leading  edge  of  a 
hypersonic  body.  The  viscous  boundary  layer  which  is  u.suany  thin,  occupic.s  an 
appreciable  fraction  of  the  distance  between  the  shock  and  body  that  will  be  considered 
without  undue  loss  of  generality  a  flat  plate  in  what  follows.  Accordingly  f(x,f)  =  0,  in 
the  notation  of  Fjg.  1.  The  results  in  this  chapter  will  be  expressed  in  terms  of  the 
boundary  layer  thicknes.s  function  A(j,7)  =  0 ,  which  in  the  interpretation  mentioned  in 
the  Introduction  could  be  the  body  shape  in  an  inviscid  context. 
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The  unsteady  form  of  the  Hypersonic  Small  Disturbance  Theory  (HSDT)  equations  (9J 
are  applicable  and  are  obtained  as  in  [7]  from  limit  process  expansions  of  hatted  variables 
defined  as  quantities  normalized  by  their  freestream  counteiparts,  with  p,T,u,v,^  the 
density,  temperature,  horizontal,  vertical  components  of  the  velocity  vector,  and  viscosity 
respectively.  If  the  freestream  density,  pressure  and  velocity  are  denoted  asf/,p_and 
respectively,  then  a  pressure  coefficient  used  in  these  expansions  is  defined 
a!ip=(p-pJlp_U\ 


tC.t€iM9  90*4 


Fig,  1  Schematic  of  hypersonic  strong  interaction  flow. 

With  thc.se  definitions  and  the  coordinate  system  in  Fig.  1  as  well  the  normalization  of  the 
Cartesian  dimensional  coordinates  I  and  y  to  the  unit  reference  length!  and  the 
reference  time  scale  L/f/  for  the  time  7,  unbarred  dimensionless  normalized  counterparts 
of  these  independent  variables  are  defined.  If  Af.  and  /?_  are  respectively  the  freestream 
Mach  and  Reynolds  numbers,  and  5  i.s  a  characteristic  flow  deflection  angle,  then  the 
expansions  are 


p(x,y.i\M^.R^.r)  =  a[x,y,t.H,rV" 

(1.1) 

f=T+^- 

(1.2) 

p  =  5^p+... 

(1.3) 

M=  1+  • 

(1.4) 

(1.5) 

where  >•  =  y/(  W) .  These  expansions  are  valid  in  the  HSDT  limit 


x,y.t,H  ■» 


1 


are  fixed  as  5  -♦  0 , 


where  H  i.s  the  hypersonic  similarity  parameter. 
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In  [7]  and  elsewhere,  e.g.  [1 J,  [8]  it  was  shown  that  S  =  wheree  =  Re'  0  for 
strong  interaction  flows. 

Substituting  (1.1)  -  (1.6)  into  the  Navier  Stokes  equations  and  shock  conditions  gives  the 
following  HSDT  unsteady  equations  and  boundary  conditions; 


f  ^  5  ^  dav 

[dt'^dxr'^  dy 
( d  B  B^  1  ^ 
[Bt'^Bx'^''Byy~aBy 


(1.8) 

(1.9) 


B  B 

- - +  V~“ 

t  Btc  By 


(1.10) 


In  what  follows,  the  boundary  layer  will  regarded  as  an  effective  body,  since  in  the  strong 
interaction  regime,  the  density  of  the  flow  approaches  zero  as  the  edge  of  the  boundary 
layer  is  approached(8],  fl],  and  [4]  so  that  the  latter  i.s  streamline  in  steady  flow  In 
unsteady  flow,  it  is  plausible  that  it  is  a  streakline.  From  Fjg.  !,  denoting  the  ordinate  of 
this  surface  as  y  =  Sdk{xJ)  then  following  [3]  for  steady  flow,  the  body  coordinates 
introduced  for  later  use  in  the  Newtonian  limit  arc 


(1.11) 
(1.12) 

(1.13) 

(1.14) 

(1.15) 

(1.16) 

The  quantity  v,  represents  the  vertical  velocity  component  at  the  boundary  layer  edge. 
Equation  (1.12)  is  derived  from  the  property  of  this  surface  that  it  is  a  strcaklinc. 
Representing  this  line  in  the  implicit  form  B(i,y,f)  =0,  where 


y'  =  y-A{x,t) 

v'  =  v-v^ 

V,  *v(jr,A)  =  A,+A, 
x'  -X 

O^y*  5G  =  g-A. 


R  =  y-A(r,r)aO,  (1.17) 

the  strcaklinc  property  can  be  expressed  as 

R,+^.VR»0,  (1.18) 

where  q  denotes  the  flow  velocity  vector.  Substitution  fur  q  (1.4)  and  (1..5)  into  (1.18) 
gives  (1.12). 
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Defining  the  convective  operators  as 

(1.19) 

(1.20) 


the  HSDT  equations  written  in  these  coordinates  are 

Z)<T  +  (tJV*)  .  =0 
Ev\D^^  =  a''p^. 

V  o''  ) 

If  subscript  S  denotes  evaluation  on  the  shock,  the  shock  conditions  in  these  coordinates 
are 

(1.24) 

(1.25) 

(1.26) 

Recognizing  that  the  shock  layer  is  0{y  - 1)  as  y  I  ,and  //-♦«>,  even  in  the  unsteady 
case,  the  following  Newtonian  limit  process  and  expansions  are  used: 


[D(G+<J)]' 


(1.21) 

(1.22) 

(1.23) 


rs  ^  ^ 

d  d  d  d  d  d 

dt’'^  dx'^'' dt'^ dy' 


^  H  .  ,  ,  y’  r  A  y-’ 

N»~r,t  ,x  are  fixed  us  A  ■ - - 

X  X  y  +  1 

p(x\y',t')  =  p(x*. )?,?*)+•.• 


V  =  AF  +• '  • 

g(x’.y’)  =  A(x',r*)  +  As^(x*,/*)+... . 


(1.27) 

(1.28) 

(1.29) 

(1.30) 

(1.31) 


where  N  is  the  Newtonian  similar*”  parameter. 

Substitution  of  (1.28)  -  (1.31)  into  (1.21)  -  (1.23)  and  (1.24)  -  (1.25)  as  well  us  retaining 
like  orders  gives  the  follcr.ving  unsteady  Newtonian  equations  and  sh(x;k  conditions; 
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(1.32) 

(1.33) 

(1.34) 

(1.35) 


On  the  shock, 

N 

Vs^v(x,g(x.t),t)  =  Dg—~--DA 

Uix 

(OA)^ 

A^  +  (Z)A)- ' 

Eq.  (1.12)  implies  the  tangency  relation 

v(x*,0,/*)s0.  (1,39) 

The  initial  boundary  value  problem  comprised  by  (1.32)-(1.35),  (1.36)-(1.38)  and  (1.39) 
can  be  simplified  by  an  un.stcady  extension  of  the  Von  Mises  transformation.  Introducing 
an  unsteady  stream  function  y/ 

yf;=d,  ~D\tf  =  dv,  (1.40) 

an  unsteady  extension  of  the  Von  Miscs  mapping  is  considered  in  which 

(x',5f.r')-»(^,7j.(-/).  (141) 

with  characteristic  coordinates  defined  a.s 
e  x'-r* 

2“*.  (1.42) 

the  differential  operators  map  as  follows; 


(1.36) 

(1.37) 

(1.38) 


Ed+av.  =0 


E  eD+v-^ 

+  =0 
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p  =  -v 


(1.55) 


where  p,  =  p{x',t’,0)  is  the  pressure  at  the  boundarj’  layer  edge.  The  quantity  y/^  can  be 
obtained  from  the  second  member  of  (1.43)  specialized  to  the  shock  value  of  the  stream 
function  ,  and  the  shock  relations  (1.50)  -  (1.52).  An  integration  with  the  sissumption 
that  the  shock  intersects  the  boundary  layer  and  body  at  the  nose  gives 

yfj=A 

(1.56) 

Evaluating  (1.55)  on  the  shock  gives  an  uasteady  extension  of  the  Newton-Busemann 
pressure  formula.  This  interaction  law  is 


P,  =  (DA)*  +  A£)^A. 

Evaluation  of  K  at  the  shock  gives 

^  =  N  ■<r{DA)' =  Kix' -t\y/^). 


(1.57) 


(1.58) 


which  with  (1.56)  gives  a  pair  of  relation.^  which  parametrically  defines  the  spatial 
entropy  distribution  at  each  time  instant  along  the  various  instantaneous  streamlines.  The 
mapping  >•-♦  yr  is  obtained  from  the  first  member  of  (1.43),  as  well  as  (1.53),  (1.55), 
and  an  integration  as 


Evaluating  ( 1.59)  at  y/,  gives  the  shock  shape 
g(x\t*)  =  J{x’,t',A). 

The  density  is  obtained  from  the  entropy  integral  ( 1 .33)  and  (1 .55)  as 


(1.39) 


(1.60) 


K{x'-t\yf)- 

Writing  K{x’  -  r*,yr)  ■  Q(rj,y/)  ,  the  vertical  velocity  is  derived  from  (1.55),  the 
boundary  condition  (1 .39),  imd  integration  of  (1.46).  This  gives 


(1.61) 


,  f-j.fgm.jo  I  ■ 

j.  ■ 


(1.62) 
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From  the  aforementioned  attachment  of  the  shock  to  the  body  leading  edge  property 
implied  for  continuum  flow,  evaluation  of  Kat  the  nose  j:  =  0and  the  shock  relations 
(I.50)-n.52)  give  the  density  on  the  boundary  layer  edge  to  be 


N  +  (DA)l, 


(1.63) 


5.3  Computstlons  and  results 

The  mapping  y  given  by  (1.59),  can  be  expres.scd  in  explicit  form  convenient  for 
numerical  calculations.  Infoducing  independent  variables  x  =  x’  and  s  =  x’  -t‘  gives 

dV/(s  =  const.)  =  ^-d!T+^~d/*  =:DA(x,jr-j)ir.  (2.1) 

dx  a 

Substituting  (2.1)  into  fl. 59)  and  accounting  for  (1 .55)  and  (1.58),  gives 
fA^  (x,,/  +  jr,-x)  +  (v]A,(x,,r  +  x, -jr) 

y  -  I  — r - ^ - T - dx,  (2.2) 

Ju  A^(j:,r)-[A,(j:,,/  +  jr,-x)-A(x,r)jA2(x./) 

tfr  =  A<x,.x,+r-x)  (2.3) 

A,(x,0  =  DA(x,r),  Aj(x.r)  =  D’(x,/) .  (2.4) 

The  flow  temperature  r,  pressure p ,and  density  a  at  the  pint  (x, y)  and  the  time t  axe 

f  ~  A\{x,,x,-¥t- x)  + N  (2.5) 

p  =  p,(x./)-Ai(x./)vr  (2.6) 

d  =  plf.  (2.7) 

The  shock  deviation  is  evaluated  by  the  integral  (2.2)  at  the  upper  limit  of  integration  .r  . 
Equations  (2.2)  and  (2.3)  indicate  that  the  stream  function  \p  and  the  temperature  T at  the 
pint  (x,y)  correspnd  to  those  just  behind  the  shock  at  the  cross  .section  x,  where 
y  =  ^(x,./). 

If  the  boundary  layer  edge  A(x,r)  is  known  function,  an  unsteady  flow  field  within  shock 
layer  can  be  easily  calculated  u.sing  Hqs.  (2.2)  and  (2.5).  As  an  example,  consider  the  case 
when  the  boundary  layer  thickness  is  much  less  the  body  thickness.  Then 

A(x,/)»F(x./)  +  f7(6).e=<5,/<5-»0.  (2.8) 


Unsteady  prturbations  poduced  by  the  body  .surface  in  the  zeroth-order  approximation 
with  respect  to  the  small  parameter  6  ,  will  be  considered.  Features  of  the  flow  field  will 
be  discussed  for  the  body  shap  examples 
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A  =  ;r+y4x«-"-^’’*cos(<ur)  (2.9) 

A  =  x^^A — - — sin(ca-£ur).  (2.10) 

X  +  Xq 

Shape  (2.9)  dc-scribe;:  a  wedge  of  a  half-angle  <5  with  a  2-D  hump  placed  at  a  distance 
Xf) flora  the  leading  edge  and  oscillating  with  frequency  O)  and  amplitude  A.  Shape 
(2.10)  de.scribes  a  wedge  with  a  superposed  traveling  wave  modulation  starting  near  the 
point  x^and  propagating  downstream  with  phase  velocity  c  =  ta/or  ;  the  wave  amplitude 
achieves  the  constant  value  i4  as  x  . 

Rgures  2a  and  2b  illustrate  a  flow  field  past  the  steady  hump  specified  by  Eq.  (2.9)  at 
Xo=l  with  amplitude  A  =.15,  for  the  Newtonian  parameter N=0.  For  clarity  in  the  plots, 

the  temperature  profile  f(y)at  cross  seaions  x^  =  (5/12)ii:.  (i  =  l,2,...)are  plotted  with 

constant  shift  in  the  abscissae  origins  Af  =0.5  corresponding  to  the  different  k.  Figure 
2a  show.s  that  the  hump  causes  local  shock  perturbation  damping  in  upstream  and 
downstream  directions.  From  Fig.  2b,  it  can  be  deduced  that  the  curved  portion  of  the 
shock  induces  a  non-uniform  temperature  distribution  across  the  shock  layer  which  can 
be  interpreted  as  an  entropy  wake  convecting  downstream  along  streamlines.  If  viscous 
dissipation  is  weak,  the  entropy  wtike  is  preserved  in  the  far  downstream  flow  field, 
although  the  shock  and  pressure  perturbations  are  infinitesimal. 

Similar  data  arc  shown  in  Figures  3a  and  3b  for  the  hump  oscillating  with  the  frequency 
0)  =5  and  the  amplitude  A  =.02 .  The  shock  perturbation  shape  is  close  to  that  Ibr  the 
steady  ca.se  having  an  instantaneous  body  surface  of  the  same  shape  (c/.  solid  curves  in 
Fig.  2a  and  Fig.  3a).  Such  a  quasi-stcady  behavior  is  due  to  the  infinitely  large  speed  of 
pressure  disturbance  propagation  across  the  shock  layer  according  to  Cq.  (1.34). 
However,  the  temperature  and  density  disturbances  induced  by  the  curved  shock  convect 
downstream  at  frecstream  speed  (see  Eqs.  (1.32)  and  (1 .34)).  Because  of  the  time  lag  and 
interference  between  entropy  perturbations,  the  temperature  profiles  are  quite  different 
from  those  for  the  steady  exse  (compare  Fig.  2b  and  Fig.  3b).  They  oscillate  across  the 
shock  layer,  with  decreasing  wavelength  as  frequency  increases. 

In  both  steady  and  unsteady  cases  the  temperature  profiles  have  inflection  points,  that  can 
generate  an  in  viscid  instability  of  the  entropy  wake  for  disturbances  with  short 
wavelength  and  high  frequency,  (see  [10]).  It  is  assumed  that  a  steady  or  unsteady  hump 
placed  on  the  forelxxly  of  a  h>’pcrsonic  vehicle  may  be  u.scd  as  a  boundary  layer  trip  for 
the  shock  layer  flow  captured  by  the  inlet.  The  inflectional  entropy  layer  may  amplify 
short-scale  disturbances  due  to  inherent  imsubility  mechani.sms.  Such  amplification  can 
be  enhanced  by  interaction  of  disturbances  with  inlet  shock.s,  (see  [11],  and  [12]). 
Ultimately,  this  process  may  lead  to  formation  of  turbulent  wakes.  A  .similar  effect  might 
be  caused  by  blowing/suction  through  the  portion  of  the  forebody  surface.  It  is 
conjectured  that  such  inslabililies  and  tripping  mechanism  described  above  can  be 
useful  for  initiation  or  enhancement  of  hypersonic  scrantjet  combustion  processes. 
Figures  4a  and  4b  illustrate  flow  field  perturbations  caused  by  the  vibration  wave  of 
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Shape  (2.10)  for  (O  =  I0,a  =  5,  A  =.0I  and  Xo  =  l.  It  is  evident  that  the  shock  deviation 
grows  down.streara  and  ultimately  becomes  a  spike,  consistent  with  the  nonlinear  nature 
of  the  disturbance  evolution.  Additionally,  the  temperature  profiles  are  qualitatively 
similar  to  those  over  an  unsteady  hump  (cf.  Fig.  2b  and  Fig.  3b). 

As  shown  in  Fig.  5a,  the  amplitude  of  the  shock  vibration  wave  becomes  infinite  at  a 
finite  distance  from  the  starting  point  .  The  wall  pressure  is  equal  to  zero  (sec  Fig.  5b) 
at  this  location.  Since  this  is  a  nonuniformity,  the  present  asymptotic  model  is  inadequate 
to  describe  the  flow  field  near  it  However,  it  is  clear  that  the  singular  region  is  identified 
as  a  function  of  vibration  wave  characteri.siics.  The  vanishing  wall  prc.ssure  vanishing 
here  can  cause  an  inviscid  separation  from  tnc  body  surface  of  the  thin  shock  layer.  This 
process  resembles  shock  separation  from  the  upper  edge  of  hypersonic  boundary  layer  on 
a  flat  plate  at  zero  angle  of  attack  (see  [13]  signifying  that  even  small  vibration  waves 
propagating  along  the  body  surface  can  initiate  a  global  reconstruction  of  tlic  flow  field. 


5.4  Summary  and  concluahns 

In  this  phase  of  the  re.scan;h,  Newtonian  thin  shock  layer  theory  was  extended  to  unsteady 
flow.  In  earlier  work  fof  steady  strong  interaction  flow,  the  boundary  layer  thickness 
A  -  based  on  the  similarity  of  the  initial  boundary  value  problem  and  integration  of 
the  boundary  layer  continuity  equation.  An  integrodiffereniial  equation  for  A  arises  that 
admits  a  similarity  solution  with  the  scaling.  The  literature  assumes  either  a  Tangent- 
Wedge  or  Newtonian  impact  theory  interaction  law  for  the  pressure  in  the.se  .solutions. 
From  the  present  work,  It  appears  that  the  proper  asymptotically  consistent 
procedure  is  to  include  the  centrifugal  term  in  this  law  (second  term  in  the  right 
band  side  of  the  steady  state  version  of  (1.57)). 

Computations  performed  in  a  framework  of  the  unsteady  Newtonian  (hin  shock  layer 
theory  reveals  the  following  features  of  invi.scid  flow  field: 

•  A  local  stcady/unsteady  hump  induces  stcady/unsteady  entropy  perturbations 
propagating  downstream  with  freestream  velocity.  It  generates  inflection  points  in  the 
instantanesrus  shock  layer  profiles.  According  to  inviscid  stability  theory,  these 
profiles  should  be  unstable  with  respect  to  short  wavelength  disturbances. 
Down.stream  amplification  of  inviscid  instability  may  cause  formation  of  turbulent 
wakes  within  a  shock  layer.  It  is  plausible  that  such  a  tripping  mechanism  can  be 
exploited  in  forebody-inlet  design  to  initiate  or  enhance  combustion  process. 

•  Surface  vibration  waves  of  finite  amplitut'e  cause  strong  .shock  penurbafions  that 
grow  downstream.  The.se  ultimately  form  spikes  due  to  nonlinear  effects.  The  infinite 
spike  occurs  at  a  finite  distance  from  the  initial  vibration  point.  This  singularity  could 
lead  to  .separation  of  a  thin  shock  layer  from  the  body  surface  and  global 
reconstruction  of  the  flow  field  within  the  shock  layer. 
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Fig.  2a.  Shock  deviations  caused  by  steady  bump,  A  *.15,  ato  =  l.tw  =  0,  ^  =  0, 
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Fig.  3a.  Shock  deviation.^  caused  by  unsteady  bump,  A  =.02,  jTq  =  l,ft)  =  5,  =  0. 


Fig.  3b  Temperature  deviations  caused  by  unsteady  bump,  A  =.02,  =  l,eo  =  5,  N  =  0. 
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Fig.  4a.  Shock  deviations  caused  by  traveling  vibration 
wave,  A  =.01,  =  l,fi)  =  10,  =  0,  a  =  5 . 


T 


Fig.  4b.  Tcmfjcrature  deviations  cau.scd  by  traveling  vibration 
wave,  A  =.01,  x,,  =  l,<y  =  10,  =  0,  a  =  5. 
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Fig.  5a.  Shock  deviation.?  caused  by  traveling  vibration 


wave,  A  =.0 1,  x„  =.2,  o)  N  =  0,  a  ==5. 


Fig.  5b.  Temperature  deviations  caused  by  traveling  vibration 
wave,  A  =.01,  .Xg  =.2,0)  =.5,  N  =  0,a  =  5. 
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6.  Similarity  Methods  for  Hypersonic  Transition  Prediction 

S.llntroducthn 

Prediction  of  laminar-turbulent  transition  in  hypersonic  boundary  layers  is  of  critical 
importance  for  aircraft  and  missile  design  [1].  Although  considerable  effort  has  been 
invested  in  developing  theoretical  and  computational  methods  to  predict  high  Mach 
number  transition,  a  strong  need  still  exists  for  approaches  suitable  for  rapid  response 
design  application.  To  fill  this  need,  empirical  criteria  based  on  wind  tunnel  and  flight 
experiments  are  currently  the  main  workhorses  in  engineering  practice.  Much  of  the 
difficulty  in  constructing  a  rapid  design  oriented  method  is  due  to  the  unique  features  of 
the  hypersonic  stability  and  transition  problem  that  include  strong  non  parallelism  and 
shock  effects.  In  spite  of  the  success  of  methods  such  as  [2]  below  Mach  12,  these 
aspects  are  significant  at  higher  Mach  numbers. 

In  this  chapter,  similarity  methods  will  be  used  to  indicate  a  potentially  new  useful . 
procedure  for  preliminary  design  and  fast  tradeoff  studies.  These  methods  have  the 
potential  of  dealing  with  some  of  these  difficulties.  The  theoretical  basis  for  the 
similarity  method  will  be  given  and  an  application  to  hypersonic  cones  provided.  Finally, 
possibilities  for  generalization  of  the  method  to  more  aititrary  shapes  will  be  discussed. 

Since  the  problem  of  disturbance  evolution  contains  the  large  Mach  and  Reynolds 
number  parameters,  it  is  natural  to  use  asymptotic  methods  in  predictive  models. 
Asymptotic  theories  of  the  vorticity  mode  as  well  as  cross  flow  instability  and  Gbrtler 
vortices  as  well  as  [3]  are  recent  examples.  Solutions  from  this  class  of  asymptotic 
models  are  limited  in  their  applicability  because  of  the  restrictive  asymptotic  limit 
processes  used.  Other  asymptotic  approximations  provide  more  general  applicability  but 
are  also  more  difficult  to  solve.  Even  without  solution  however,  such  theories  give 
important  similarity  groups,  if  appropriate  normalizations  and  nondimensionalizations 
are  used  that  reflect  the  salient  phenomenological  scales.  These  groups  provide  the  basis 
for  constructing  similarity  laws  that  not  only  can  be  useful  for  design  of  economical 
experiments,  but  prediction  as  well. 

Such  a  procedure  will  be  used  in  this  paper  to  investigate  similitude  of  boundary  layer 
stability  over  nearly  sharp,  zero-incidence  hypersonic  cones.  A  similarity  nile  accounting 
for  viscous  and  compressibility  effects  is  derived  without  solving  the  initial  boundary 
value  problem.  Predictions  from  the  rale  are  compared  with  linear  stability  computations 
as  well  as  (Reentry  F)  flight  data  obtained  on  slender  cones  [4].  It  will  be  seen  that  this 
rule  gives  an  extremely  cost-effective  and  rapid  extrapolation  method  of  code  or 
experimentally  derived  data  to  other  conditions. 
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6,2  Similarity  analysla 

The  unsteady  disturbance  field  in  a  three-dimensional  hypersonic  boundary  layer  will  be 
considered  with  the  following  assumptions;  (i)  The  fluid  is  a  perfect  gas  with  constant 
Prandtl  number  Pr  and  specific  heat  ratio  y.  Denoting  starred  quantities  as  dimensional, 

the  viscosity-temperature  law  is  T  .  (ii)  The  spatial  scales  of  the  disturbances 

are  of  the  order  of  boundary  layer  displacement  thickness 5';  the  time-scale  isS‘  I U',, 
where  subscript  e  refers  to  the  upper  edge  of  the  boundary  layer,  (iii)  Boundary-layer 
disturbances  weakly  interact  with  the  external  inviscid  flow  and  shock  wave. 


In  the  hypersonic  boundary  layer,  the  pressure  is  of  the  order  of  that  of  the  external 
inviscid  flow,  and  the  temperature  is  of  the  order  of  the  stagnation  temperature 
77  -  (y  -  This  leads  to  asymptotic  scaling  for  the  flow  variables  within  the 


boundary  layer :  The  full  Navier  Stokes  equations  written  in  terms  of  this  scaling  contain 
only  the  lumped  parameter/?  =  wherec*^^ — «  land  Reis  the  Reynolds 


number  based  on  the  displacement  thickness^  Re  =  / p]  / ^i*,.  The  dependence  on  the 
similarity  parameter  R  is  valid  for  strong,  moderate  and  weak  viscou.s-inviscid 
interaction  regimes.  However,  <5*  and,  as  a  result,  the  Reynolds  number  Re ,  are  sensitive 
to  the  interaction  type.  For  example,  the  boundary  layer  on  a  flat  plate  has  the 
displacement  thickness  6‘  -  for  strong  interaction  and^  -  x**'*  for  a  weak  one. 


Outside  the  boundary  layer,  the  disturbances  are  assumed  to  vani.sh  and  weakly  interact 
with  the  external  inviscid  flow  .  At  the  wall,  the  flow  satisfies  no-slip  conditions, 
temperature  is  assumed  to  be  constant.  The  upstream  and  downstream  conditions  as  well 
as  initial  conditions  for.the  disturbances  are  also  assumed  to  depend  weakly  on  Mach  and 
Reynolds  number.  These  additional  provisions  preserve  the  applicability  of  the  R 
similitude  previously  indicated.  It  is  therefore  clear  that  the  similarity  law 


f^MP^rJnPr)  (1) 

holds,  where  /  is  any  of  the  flow  dependent  variables.  7)  is  the  wall  temperature  ratio, 
and  signifies  a  functional  dependence  that  can  be  obtained  from  .solution  of  the  initial 
value  problem. 

This  similarity  cun  be  extended  to  the  linear  stability  problem  if  a  typical  eigenmode 
wavelength.  A'  is  of  the  order  of  the  boundary  layer  thickness.  If  small  disturbances  are 
assumed  in  the  original  form  of  the  Navier  Stokes  equations  before  the  previously 
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mentioned  rescaling  is  used,  the  general  linearized  equations  (LSE)  for  small  fluctuations 
are  obtained.  Within  the  spatial  and  temporal  scaling  implied  by  the  wavelength 
assumption,  the  equations  for  the  fluctuations  depend  again  only  on  the  parameters 
y,Tf,PrziidR,  if  further,  their  initial  and  boundary  conditions  are  consistent  with  this 

assumption.  Therefore,  the  fluctuations  obey  a  similarity  law  such  as  (1). 

As  an  example,  consider  the  two-dimensional  vorticity  mode  (second  mode  according  to 
Mack’s  classification  [5])  in  the  boundary  layer  on  a  flat  plate  or  sharp  cone.  In  this  case, 
the  span  component  of  the  wave  vector,  ^  =  Oand  the  growth  rate  a  -  -Im(a)  (a  =the 
complex  eigenvalue).  If  the  Prandtl  number,  specific  heat  ratio  and  vi.scosity  exponent 

r* 

parameters  Pr,y,fl),  as  well  as  the  wall  temperature  factor  =  — ^are  fixed,  where 

is  the  wall  tcmperanire,  ^  is  the  adiabatic  wall  temperanirc,  the  maximum  growth  rate 
<t(£2),  (Q  =  Cl‘S*  /  V]  the  non  dimensional  frequency)  is  a  function  of  Reynolds  number 
/?e  and  Mach  number  A/, . 

In  Fig.  1,  calculations  for  or.  from  a  linear  stability  code  arc  shown  as  a  function  of  the 
similarity  parameter  R  at  =  7, 8, 10, 12  and  14.  The  following  parameters  were  used 
in  the  computations:  Pr  =  0.72,  y  *  1-4,  0.2  and  0=0.75.  The  stability 

characteristics  have  been  computed  by  solving  the  eigenproblem  for  the  LSE  system  of 
equations  [5].  In  this  procedure,  the  mean  flow  in  the  boundary  layer  was  approximated 
by  the  hypersonic  weak  interaction  theory  self-similar  solution.  Note  that  the  similarity 
given  by  (1)  collapses  the  ajRe,M,)  curves  to  O’.(R)  at  moderate  values  of  /?,  as 
Mach  number  increases.  At  sufficiently  large  R,  the  curves  tend  to  spread,  because  the 
wavelength  of  the  vorticity  mode  X « S'  and  Assumption  (ii)  is  violated.  This  is 
con.sistent  with  the  asymptotic  theory  for  inviscid  disturbances  according  to  which  the 
vorticity  mode  is  concentrated  in  a  thin  transitional  layer  at  the  upper  edge  of  the 
boundary  layer.  However,  the  region  where  similarity  works  well  is  sufficiently  wide  and 
relates  well  to  practical  applications. 


Mach  numbers  Ai, . 
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6.3  Appllcstlon  of  similarity  In  hypsraonlc  transition  prediction 

Hypersonic  cones  are  used  as  a  focus  to  explore  the  use  of  similarity  in  transition 
prediction  herein.  Reference  [4]  provides  transition  data  for  a  slender  cone  at  freestream 
Mach  number  20  during  reentry  from  altitudes  of  approximately  30.480  to  18.288  km, 

corresponding  to  the  region  of  freestream  Reynolds  number  6.56  to  52.5  x  10^  per  meter. 
Experimental  points  representing  transition  onset  Reynolds  numbers  Rc^^  =  (/‘x*  /  v) , 

(v)=  kinematic  viscosity)  as  a  function  of  local  Mach  number  M,  are  shown  by  the 
symbols  in  Fig.  2  In  the  figure,  I  denotes  the  sharp  cone,  n  signifie.s  entropy  layer 
swallowing  by  the  boundary  layer  close  to  the  transition  location  and  III  represents 
swallowing  downstream  of  transition.  The  similarity  relation  (1)  applies  for  Regimes  I 
and  m. 

10* 


10* 

Fig,  2  Comparison  of  present  theory  with  flight  experiment.  Solid  line:  theory, 
symbols:  experiment. 


Because  of  recession  of  the  nose  tip  during  reentry,  the  local  Mach  number  at  the  upper 
edge  of  the  boundary  layer  depends  on  the  nose-radius  history.  The  maximum  M, 
corre.sponds  to  a  sharp  no.se.  As  the  tip  is  blunted,  the  boundary  layer  grows  in  a  higher 
entropy  gas  with  smaller  local  Mach  number.  The  transition  Reynolds  number  depends 
on  the  location  of  the  region  where  the  entropy  layer  is  swallowed  by  the  boundary  layer. 
More  details  are  given  in  Ref.  (6]  of  how  this  physics  is  related  to  the  change  in  the 
monotonic  trend  of  the  Regime  ID  data  shown  in  Fig,  2  .  For  Regime  II,  only  local 
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similarity  rather  than  the  global  similarity  (1)  is  valid  that  is  associated  with  almost 
constant  boundary  layer  edge  flow  quantities.  The  latter  corresponds  to  sharp  cone 
similitude  at  a  given  local  Mach  number.  For  this  case,  the  flat  plate  mean  flow  viscous 
interaction  parameter;^  ■(/-!) ■£/,*/ /v)),  derived  in  Ref.  [7]  is 
related  to  transition  onset.  The  assumption  of  weak  interaction  « 1  is  well  satisfied  by 
experimental  data  of  [4].  From  [7],  and  application  of  the  Blasius  formula  locally,  the 
boundary-layer  displacement  thickness  local  Reynolds  number  Re  and  similarity 
parameter,  Rare 


s'=£ 


Re^Ae  ^  (3) 

using  the  perfect  gas  law,  energy  invariant,  and  the  previously  staled  viscosity- 
temperature  law,  where  A  -  the  constant  of  proportionality  in  this  law. 


In  (2)  and  (3),  A  is  a  weak  function  of  M,.]n  the  hypersonic  limit,  A  is  assumed 
independent  of  M, .  This  approximation  simplifies  the  prediction  of  the  transition  onset 

location  that  can  be  estimated  by  thee^  method.  The  amplification  factor  N  for  the  two- 
dimensional  second  mode  can  be  represented  solely  as  a  function  of  the  similarity 
parameter  Ras 


where  <y  *•  cS’ .  Applying  the  transition  onset  criterion  N=:N„,  where  is  the  critical 
amplification  factor  that  is  assumed  to  be  a  constant  for  Regime  III  data  from  [4],  the 
corresponding  similarity  parameter  R„  can  be  obtained.  According  to  (3), 

Re,j,  =  .  Since  N„  sconst  independent  of  Mach  number  Af,  and  from 


N=:ff„{R),  R„  also  =const,  independent  of  Mach  number  A/,.  Therefore,  the 
similarity  rule 


is  obtained  on  use  of  the  definition  of  e  and  the  fact  that  the  other  parameters  are  fixed. 


Equation  (4),  shewn  in  Fig.  2  as  a  solid  line  for  0.75,  agrees  excellently  with  flight 
data  and  provides  a  useful  prediction  tool  that  can  be  applied  for  other  zero  incidence 
cones.  With  (4),  characteristics  at  new  Mach  numbers  can  be  obtained  from  computed  or 
experimental  points  corresponding  to  a  particular  cone  geometry  and  wall-to-freestream 
lemperanire  ratio. 
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6.4  Discussion,  future  work  snd  conclusions 

Effort  is  in  progress  to  extend  the  applicability  of  the  similarity  method  to  more  general 
shapes  than  cones  and  treat  angle  of  attack  cases.  In  this  connection,  it  should  be  noted 
that  similarity  gives  a  relation  between  the  local  characteristics  of  the  boundary  layer.  The 
simple  form  of  the  similarity  rule  (4)  was  obtained  under  the  assumption  that  the  mean 
flow  parameters  at  the  boundary'-layer  edge  are  constant.  Similarity  appears  also  to  be 
valid  for  transition  due  to  some  local  process  that  is  typical  of  nonlinear  by-pass 
mechanisms.  If  transition  depends  on  the  disturbance  streamwise  history  and  the  mean 
flow  has  a  non  uniform  distribution  in  the  x,  2  plane,  similarity  can  be  used  to  relate 
disturbance  characteri.stics  between  two  streamwise  stations  and  account  for  the  global 
process  leading  to  transition.  From  thee*'  method,  the  growth  rate  (T  can  be  extrapolated 
to  the  new  set  of  conditions  at  each  point  using  similitude.  This  can  be  used  to  determine 
the  amplification  rate  N  to  predict  the  transition  Reynolds  number.  Since  this  approach 
would  be  derived  under  the  other  general  conditions  used  here  it  would  not  only  be 
applicable  for  more  complicated  geometries  but  linear  and  nonlinear  regimes  as  well, 
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7.  Transonic  Wind  Tunnel  Wall  Interference 

T.IIntroducthn 

A  central  problem  in  aerodynamics  is  prediction  of  vehicle  flight  characteristics  from 
wind  tunnel  tests  of  subscale  models.  Major  factors  in  extrapolating  the  tunnel  results  to 
free  flight  arc  Reynolds  number  scaling  as  well  as  wall  and  sting  interference.  Currently, 
the  issue  of  U.S.  wind  tunnel  modemi2ation  is  being  addressed.  A  major  thrust  is  attain¬ 
ment  of  near-flight  Reynolds  numbers.  If  large  models  are  used,  wall  and  sting  interfer¬ 
ence  are  limiting  factors  in  obtaining  a  tunnel  flow  even  qualitatively  resembling  that  of 
flight.  In  particular,  testing  at  transonic  speeds  can  produce  steep  wave  fronts  that  reflect 
back  on  the  model.  Porous  and  other  types  of  ventilated  walls  have  been  developed  with 
the  aim  of  canceling  the.se  rcflc  ions.  In  spite  of  the  advances  made  in  this  technology  as 
well  as  computational  simulations,  much  still  needs  to  be  done  to  understand  the  tradeoffs 
in  sizing  wind  tunnel  models  and  test  sections  to  minimize  wall  and  sting  effects  while 
maximizing  model  Reynolds  numbers. 

References  (l]-[3]  exemplify  early  treatments  of  lifting  flat  plate  wings  in  porous  tunnels 
at  subsonic  speeds.  This  work  did  not  address  important  nonlinearities  as.sociatcd  with  the 
transonic  regime  and  assumed  small  wing  dimensions  compared  to  that  of  the  test  section. 
No  framework  was  developed  for  systematic  improvement  of  the  approximation,  and  to 
our  knowledge,  the  results  were  not  compared  with  experiment.  Computational  treat¬ 
ments  of  transonic  flows  are  discussed  in  [4] -[6]  which  remedy  some  of  these  deficien¬ 
cies  and  deal  with  realistic  airplane  shapes.  The  large  scale  (CFD)  computational  ap¬ 
proach  can  model  geometric  complexity  and  nonlinear  effects  but  is  usually  expensive 
and  time-consuming.  Although  considerable  progress  has  been  made,  much  user  inter¬ 
vention  is  still  required  for  geometric  pre-processing  to  mate  structured  and  even  unstruc¬ 
tured  grids  to  the  3-D  surface  topology.  This  is  an  important  i.ssue  for  the  test  engineer 
who  needs  to  quickly  size  the  model,  make  test  plan  adjustments  during  the  test  and  in¬ 
terpret  the  results  afterwards.  Combined  asymptotic  and  computational  (CAN)  methods 
are  discussed  in  {7]  and  [8]  that  address  this  ne^. 


Fig.  1  Schematic  of  confined  .slender  airplane. 
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These  approaches  reduce  the  size  and  computational  intensity  of  the  nonlinear  problem 
and  provide  frameworks  for  approximation  procedures  that  arc  capable  of  systematic  re¬ 
finement 

As  examples  of  the  CAN  approaches,  [7]  and  [8]  discuss  transonic  blockage  interference 
for  confined  slender  shapes  such  as  those  shown  schematically  in  Fig.  1  and  lift  interfer¬ 
ence  for  high  aspect  ratio  wings. 

In  this  section,  improvements  of  .^ur  previous  theories  will  be  discussed  in  four  major  ar¬ 
eas  delineated  in  the  contract  St-tement  of  Work  (SOW).  These  are: 

1.  Treatment  of  lift  interference 

2.  Extension  of  large  wall-height  blockage  interference  to  moderate  height  case 

3.  Porous  wall  effects 

4.  Non-circular  wind  tunnel  test  sections 

7.2  Utt  Inteiiennce  and  peroua  wait  affecta 

The  thrust  of  thi.s  activity  is  to  develop  a  systematic  asymptotic  framework  for  computa¬ 
tion  of  lift  corrections  due  to  the  interaction  of  a  slender  model  with  walls.  Because  of  the 
resemblance  of  the  asymptotic  developments,  to  those  for  transonic  flow,  the  subsonic 
case  was  considered  for  convenience.  As  indicated  later,  this  approach  actually  seemed  to 
provide  good  comparison  with  experiment  at  near-sonic  speeds.  The  treatment  was  gen¬ 
eral  enough  to  address  Item  3. 


Fig.  2  Schematic  of  slender  wing. 
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Initial  developments  under  this  contract  are  described  in  [9].  There,  the  free  field  and 
wind  tunnel  problem  of  the  incompressible  flow  over  a  flat  wing  of  arbitrary  planform  in 
a  circular  wind  tunnel  test  section  was  outlined.  An  in-depth  continuation  of  that  intro- 
ductoiy  treatment  will  be  summarized  in  what  follows,.  Expansions  of  limit  process  type 
have  been  used  to  the  study  the  matching  process  in  greater  detail  than  in  the  preliminaiy 
analysis  of  [9].  In  addition,  the  second  order  inner  approximation  has  now  been  derived. 

Letting  be  the  velocity  potential,  limits  involving  the  semispan  of  the  wing  b  and  the 
angle  attack  a  have  been  considered.  Near  the  wing,  a  limit  process  in  which  b  0  is 
us^.  Refening  to  Fig.  2,  an  inner  limit  is  defined  as 

^^^  =  e  +  flA(jO,(x,y*,z‘)+flbMog-^4)j,(j:,y*,z')  + 

U  (1) 

where  the  inner  limit  is 

a  =  tana,  A  =  a/b, y*  s y /b,z*  s z/b, fixed  asfl,b 0  (2) 

In  (2),  the  characteristic  wing  chord  is  fixed  while  the  semispan  b  and  angle  of  attack  a 
tend  to  zero  at  the  same  rate.  Near  the  wing,  cross  flow  gradients  dominate  and  these  pa¬ 
rameters  give  the  characteristic  lateral  scale  of  the  flow  which  is  b.  Equation  (1)  is  an 
inner  expansion  for  the  velocity  potential  0  in  terms  of  approximating  perturbation  po¬ 
tentials^,,  ,(v=  order  of  the  approximation).  It  contains  the  “switchback”  term  tpj,  and 
the  indicated  gauge  functions  in  anticipation  of  matching. 

The  problems  for  the  9,,  are  obtained  by  .substituting  the  asymptotic  developments  into 

the  exact  problem  for  .  For  this  exact  problem,  the  equation  of  motion  for  the  velocity 
potential  is 

V^<ti  =  0.  (3) 

Bernoulli’s  Equation  is 

=  (4) 

The  equation  of  the  body  is 

B  =  yi-axssO.  (5) 

On  this  surface,  the  tangency  b<}undary  condition  is 

qVB  =  Q.  (6) 

Darcy’s  law 

AP-RV,  (7) 


is  used  to  model  the  porous  walls,  where. 


7-3 


SC71062.FTR 


AP  =  Pressure  difference  across  wall 
R  =  Eow  resistance  constant 
V  =  Transverse  velocity. 

The  dominant  orders  in  (1)  solve  harmonic  boundary  value  problems  In  the  cross  plane 
perpendicular  to  the  frcestream  and  the  higher  orders  solve  Poisson  problems.  In  particu¬ 
lar,  sub.stitution  of  the  inner  expansion  into  the  exact  problci.i,  gives  the  following  prob¬ 
lems: 


=<P,yy  +<P,,v 
(p^.(x,-Aa,z')  =  ~l 


(8a) 


-  0 


(p2^.(x,-Ax,z‘)  =  -<p^^(x-Ax,z')  (8c) 

Problems  (8)  require  far  field  boundary  conditions  for  r‘ » +2'*  -»  «o .  These  come 
from  matching. 

The  dominant  inner  approximation  from  (8a)  provides  a  first  estimate  for  the  flow  and 
pressure  field  of  the  wing.  The  leading  edge  square  root  singularities  dominate  this  flow 
field  which  is  the  stagnating  flow  on  a  finite  flat  plate.  From  a  Joukowski  mapping  to  the 
circle  plane  or  the  Circle  Theorem,  the  solution  of  the  problem  (8a)  is 

qf,  =  -I Re|[cT*  -  - ct|  ct  a  2*  +/(>’’  +  Ax)  (9) 

in  which  the  frcestream  velocity  has  been  normalized  to  unity  through  the  non- 
dimensionalizations  in  (I).  This  solution  has  the  proper  far  field  (downwash  at  infinity) 
related  to  matching  with  outer  solution. 

A  similar  procedure  gives  the  following  for  solution  for  (p^2  which  is 


^1  = 


-Rc^(jj')' 


This  satisfies  the  surface  boundary  condition  in  (8b). 


(10) 


The  far  field  for  (p^  as  r  oe-is  also  needed  for  the  matching.  This  can  be  obtained  from 
particular  solutions  of  the  Poi.s.son  equation  in  (8c).  This  gives 
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^2  =“~(«')^(r‘cos0)logr* +A2(jc)3'*+  ”as  r-^o». 

Refinements  of  the  unconfined  flow  field  and  wall  interactions  come  from  coupling  with 
the  far  field.  An  outer  expansion  involving  an  0(1)  transverse  length  scale  ns  b-^0 
gives  a  semi-infinite  line  doublet  for  this  part  of  the  flow.  In  particular,  for  an  outer  limit 

x,y,z.A  fixed  as 

the  appropriate  outer  expansion  is 

■^~  =  j:  +  aA^^(x,>,2)+— .  (11) 


The  quantity  ^  solves  the  following  problem; 


s^(x) cosd 


as  r->0> 


^ll^(x,/j.0)+)3^^(x,/f,0)  =  O  J 

where  +  z^,  tan5  =  z/y  and  the  last  relation  is  the  approximate  form  of  Darcy’s 

law. 


■  +  wall  correction  function 


The  solution  of  (12)  can  be  written  as 

d  =  -  - I  — r  +  wall  correction  function  (13) 

Jo  V(^-0  +'• 

where  the  first  term  is  a  line  doublet  distribution  of  strength  D{^)  that  satisfies  the  sec¬ 
ond  relation  of  (12).  The  wall  correction  function  is  needed  to  satisfy  Darcy’s  law  in  (12). 
In  (13),  the  Kutta  condition  at  the  effective  trailing  edge  provides  the  appropriate  con¬ 
tinuation  of  the  line  doublet  to  downstream  infinity  and  its  convergent  integral  represen¬ 
tation.  Furthermore,  (13)  gives  upstream  influence  not  present  in  the  inner  solutions. 

The  inner  expansion  of  (13)  can  be  obtained  from  the  r-^0  expansion  of  its  Fourier 
transform.  TMs  gives 


D(,x)  COS0 
2/c  r 


)S0  D"(x)fr,  rl  ^  r  ^ 

—  -  — log  r  -  cosa  -  G  (x)  -  cos  a  -H  my+-  • 


as  r  ^  0 
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where 


my  s  COS0 


D{Xj)r 

2h 


P\K,q  +  K,)(!,-qIo)-^q^I,K, 
P\ql,-l,)^+q^l,^ 


.  ^(x-4)  . 

sin - ;; - aq 

h 


-  COS01 


jLf'k- 

2jrfi‘ 


0  i3’(9/o-/i)'+9^V 


cos- 


dq\ 


A  =  1  for  ^  ?t  1,  (includes  porous  and  closed  walls) 
=  -I  for  p  =  0, (free  jet),  (degenerate  limit  P-*0) 


G(x)  3  J^Z)'(4)sgn(j:  -  ^)  log  2\x  - 


D{x)s^-~a 


(16) 

where  h  is  the  wall  height  in  units  of  the  body  length  and  the  term  my  in  (16)  is  the  wall 
interaction  effect  due  to  the  imaging  of  the  doublet  in  the  walls.  As  an  alternate  approach 
for  the  first  term  of  (13),  the  integral  can  be  directly  expanded  for  r  — >0.  This  delicate 
procedure  is  described  in  (lOJ.  In  the  Fourier  integral  method  used,  the  solution  naturally 
decomposes  into  a  free  field  (no  walls  present  part)  and  a  wall  interaction  portion  as  indi¬ 
cated  in  (19).  Special  limiting  processes  of  the  singular  integrals  were  developed  to  han¬ 
dle  zero  and  infinite  porosity,  corresponding  to  solid  wall  and  free  jet  cases  respectively. 
Equation  (19)  agrees  with  the  results  from  [l]-[3]  in  the  limit  of  vani.shing  chord  to  tunnel 
radius  ratio.  In  future  work,  (19)  will  be  used  to  compute  the  camber  effect  associated 
with  non-vanishing  chord.  To  our  knowledge,  this  study  has  not  been  made  and  we  be¬ 
lieve  it  is  an  important  factor  entering  the  comparison  of  our  results  with  experiment  to  be 
discussed. 


The  inner  and  outer  solutions  match  directly  to  the  orders  to  be  considered  here.*  This 
can  be  shown  by  expressing  each  in  terms  of  an  intermediate  variable  r, ,  For  this  pur¬ 
pose,  an  intermediate  limit 

fixed  as  h -» 0 ,  (17) 

is  used  to  compare  the  inner  and  outer  representations  in  an  “overlap  domain”  =  0(1) 
in  which  both  expansions  are  mutually  valid.  Note  in  the  intermediate  limit 


‘Other  applications  exist  where  intermediate  expansions  is  required. 
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flf,,  Tl(b) 

r*  =  “7 — >«»,/•  =  Tjr„  — »  0,  ~T~  — >  ®®>  0 . 

b  '  b 

The  matching  process  consists  of  writing  inner  and  outer  expansions  in  terms  of  the  in¬ 
termediate  variable  and  comparing  similar  terms  to  determine  unknown  elements.  Car¬ 
rying  out  this  process,  the  two  expansions  appear  as  follow's  in  intermediate  variables 


Inner  =  - cos0  +  (7 

r  b 


+  ■^cosG(log  rir„  +  log  +  /ijfxlcosen  •  •  • 


o  Outer 


D"(x){ 

'nr. 

2ft  1 

.2 

2;r  m,  2;r  2 


f  G"(x))  . 

—  cos  04-  ni - - —  mr^cos0-r- 


Based  on  (18).  the  following  relations  hold 
Jab^']  D(x)  s'^(x) 


<f) 

0(ab^  log  rj); 


0(ab^T]): 

0(ab^i]logb): 


In 

2ft 


D'\x)  1  G"{x)  _  {ss')'  C'\x) 


In  A  2 
switchback  term 


+mc^k2(x) 


=><pj,  =i(ij')'y*  asr*->« 

Note  that  in  (19)  the  essential  wall  interference  effect  is  the  additive  term  m  affecting 
bjix)  in  the  0(ah'T])  matching  relationship.  Another  viewpoint  is  from  (8b),  the  solution 
of  which  consists  of  a  superposition  of  its  homogenous  solution  (response  to  homogenous 
equation  of  motion  and  boundary  conditions)  and  the  effects  of  the  forcing  terms  in  the 
equation  of  motion  and  boundary  conditions.  The  homogeneous  solution  is  non-lrivial 
because  of  the  downwash  far  field  associated  with  the  line  doublet  imaging  in  the  porous 
walls.  Another  interesting  point  is  the  surprising  appearance  of  switchback  terms.  These 
normally  are  as.sociated  with  tran.sonic  flow.  They  ari.se  in  this  subsonic  flow  from  loga¬ 
rithmic  dements  in  the  expansions. 


The  second  order  inner  refinement  of  the  free  field  lift  and  pressure  distribution  comes 
from  the  problem  (8c)  (complemented  by  a  far  field  condition)  for  the  inner  potential  <Pj. 
Techniques  u.sed  to  solve  this  problem  involved  matching  considerations  and  exclusion 
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principles  disallowing  aperiodic  logarithmic  parts  and  removal  of  unallowed  leading  edge 
singularities  to  maintain  finite  forces.'  The  inhomogeneiiy  in  the  equation  of  motion  was 
treated  by  a  particular  solution  <p^  satisfying 

where  <j  is  the  complex  conjugate  of  ff .  Other  solution  components  are  superposed  to 
satisfy  the  far  field  boundary  condition  involving  the  free  field  and  the  wall  correction  to 
the  downwash.  Our  analysis  resembles  that  of  Wang  [1 IJ  in  some  respects  who  in  con¬ 
trast  to  our  effort  did  not  consider  wall  interference  which  is  the  central  thrust  of  our 
work.  Problems  such  as  (8c  )  arc  typical  of  the  transonic  case. 


Fig.  3  Comparison  of  lift  interfcreiKc  theory  with  TsACI  experiment. 


'More  details  will  be  in  our  paper  to  be  presented  at  the  1*'  AIAA  Theoretical  Fluid  Dy¬ 
namics  Meeting,  June  l‘J96,  New  Orleans,  La. 
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Integration  of  the  pressures  on  the  wing  give  the  following  expressions  for  the  lift  L  For 
the  free  field. 


L  . 

— T  =  tana  f log-tana fp+fr  tanaf.+- 
pU  h 


(20) 


where 


f,,  (21) 

=f'  [<p,,]^z‘ =((':!/ 2X1  + log2)-4fC"  +J'‘re  T  A'] 

The  dominant  tenn  for  the  lift  f.  and  pressure  distribution  agrees  with  Jones’  theory 
given  in  [12]. 

As  an  experimental  validation  of  the  lift  interference  theory.  Fig.  3  compares  lift  versus 
angle  of  attack  predicted  by  the  our  asymptotic  theory*  with  transonic  tests  of  a  wing-body 
combination  at  TsAGI  in  Moscow  and  reported  in  19).  It  i.s  striking  that  the  incompressi¬ 
ble  theory  agrees  so  well  with  the  experiment  for  the  high  transonic  Mach  num¬ 
bers  Af  =  .99  and  1.02  tested.  Plausibility  of  this  finding  is  related  to  the  elasticity  of 
slender  wing  theory  to  not-so-slender  planforms  as  M  1  as  discussed  in  [ISJ-fH). 
Approximations  of  wall  interaction  integrals  give  the  porou.s  wall  corrections  for  wall 
openness  factors  /  =  2  and  10%  indicated  in  the  figure.  It  show.s  that  the  experimental 
trend  with  Increasing  /  is  captured  by  the  lift  interference  theory  for  vanishing  chord  to 
tunnel  radius  ratio.  However,  the  compaiison  with  the  data  shows  an  increasing  slope 
with  incidence  not  captured  by  the  first  order  theory.  Preliminary'  indications  are  that  the 
free  field  second  order  effect  shows  a  reduction  in  lift  slope  is  counter  to  experimental 
evidence.  It  is  likely  that  the  reverse  trend  is  due  to  leading  edge  viscous  separation  and 
vortex  formation  as  well  the  need  to  account  for  the  finite  chord  of  the  wing.  A  refine¬ 
ment  accomplished  is  an  estimate  of  the  effect  of  a  vortex  at  the  wing-fuselage  juncture 
occuning  at  higher  angles  of  attack.  An  oil  flow  visualization  of  this  phenomenon  from 
the  Pi’s  Russian  TsAGI  experiment  reported  in  (15)  is  shown  in  Fig  4.  Results  from  a 
preliminary  model  based  on  conical  invariance  of  the  vortex  field  is  shown  in  Fig.  5.  The 
improvement  in  agreement  is  striking  and  suggestive  of  the  importance  of  modeling  dis¬ 
crete  vortex  effects.  In  spite  of  these,  the  wall  interaction  theory  .shows  promise  of  model¬ 
ing  relative  trends.  As  in  the  blockage  theory  work  for  wall  height  of  the  order  of  the 
body  length  to  be  discu.ssed,  estimation  of  the  absolute  levels  can  be  improved  independ¬ 
ently  of  the  interference  estimations  using  vortex  dynamic  and  leading  edge  separation 
approaches  such  :is  that  just  mentioned. 


The  main  point  of  this  effort  is  that  it  provides  a  natural  launching  pad  for  extension  of 
the  theory  to  nonlinear  transonic  flow,  siccounting  for  higher  approximations,  thickness, 
viscous  interactions  and  finite  chord  to  tunnel  height  as  well  as  sy.stematic  higher  order 
refinement.  With  the  exception  of  switchback  terms  and  gauge  functions,  the  inner  prob¬ 
lems  for  the  transonic  case  arc  expected  to  resemble  those  associated  with  the  incom¬ 
pressible  asymptotic  theory.  However,  the  outer  expansions  w'ill  solve  the  three- 
dimensional  Karman-Gudcrlcy  instead  of  Laplace  equation  in  the  dominant  approxima- 
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tion,  and  forced  versions  in  the  higher  orders.  However,  a  major  .simplification  of  the 
usual  lifting  surface  (tran.sonic  small  disturbance  theory)  numerical  problem  is  anticipated 
since  the  angular  variation  can  be  separated  out  by  matching  with  the  inner  multipole 
structure. 


Fig.  4.  Wing-body-juncturc  voncx  formation  in  TsAGI  wind  tunnel. 
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Fig.5  Comparison  of  vortex  improved  theory  with  experiment. 


7.3  Extension  of  large  wall-height  blockage  Interference  to  moderate  wall 
height  case 

Major  thrusts  of  this  work  in  this  pha.se  were 

1.  Validating  a  tran.sonic  .small  disturbance  tnodcl  for  the  flow  in  a  wind  tunnel  agaiast 
experimental  data 

2.  Validating  the  equivalence  rule  for  tran.sonic  wall  interference. 
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Referring  to  Fig.  1,  the  wails  or  pressure  interface  boundary  where  pressures  are  assumed 
to  be  specified  from  experiment  are  at  r  =  A ,  where  h  is  assumed  constant  for  conven¬ 
ience  herein  (circular  test  section).*  Defining  W  =  5  A ,  the  case  if)H  =  0(1)  was  consid¬ 
ered  in  the  research  reviewed  here,  in  contrast  to  the  less  practical  situation  («)  in  [7] 
where  H  -» As  will  become  clear,  the  asymptotic  solution  of  the  Full  Potential  formu¬ 
lation,  for  Case  (i)  leads  to  two,  rather  than  three  decks  associated  with  (ii),  i.e.,  no  wall 
layer  is  required,  the  confined  flow  consisting  only  of  a  nearly  axisymmctric  “outer”  re¬ 
gion  and  a  cross  flow  gradient-dominated  inner  core  which  is  the  near  field  of  the  body. 
In  an  inner  limit  in  which 

r*  a r/5 .  3 (l - Mi)/S \  Am ajS  fixed  as5  -> 0,  (22) 

where  a  is  the  angle  of  attack,  S  is  the  maximum  thickness  ratio  of  the  equivalent  of 
body  of  revolution  and  A/,  is  the  freestream  or  tunnel  Mach  number  which  will  corre¬ 
spond  for  convenience  to  the  flight  Mach  number.  The  inner  expansion  (near  field)  of  the 
slender  airplane  model  B  =  r-6F(,x,d)  =  0  is 


,ix,r,d\M^,8,a) 


=  x  +  {7S^  Iog5 ) s; (x) -(-5  Vi +  ■ 


where  S^  (x)  is  a  source  strength  determined  by  matching  with  the  outer  solution.  The 
inner  problem  is 

=0  (24a) 

ff'M 

where  n  denotes  the  normal  direction  to  anx=const.  cross  sectional  contour  of  the  sur- 
faceBsO.  Further  details  are  in  [13,17]Application  of  Green’s  theorem  to  (24)  shows  that 


^.=-logr 


r- 
},  1 


A'W,  . 

— - — logr  +' 


i.e.,  the  asymmetric  body  appears  as  a  line  source  in  its  far  field  of  inten¬ 
sity  A'<x)  where  A(x)  is  the  cross  sectional  area  at  the  .station  x  of  the  body  B=:0.  Equa¬ 
tion  (23)  determines  3^(x)  in  (23)  as  will  be  seen  in  the  discussion  on  matching  that  fol¬ 
lows. 


The  outer  limit  (with  a  slight  change  in  notation  from  Section  2)  is 
P  mSr,  Km{\-  *,  A  h  a/5  fixed  as  5  -4  0. 

For  (26),  the  appropriate  outer  expansion  is 


=  jr  +  5’^(x,P,0  i:,A)  +  ' 


•All  lengths  are  in  units  of  the  body  length. 
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Substituting  (27)  into  the  Full  Potential  equation  gives  the  Karman-Guderley  (KG)  equa¬ 
tion 

(Ar-(r+l)«u)A„+N„),=0  (28) 

Appropriate  to  the  confined  case  considered  in  the  research,  the  approximate  Darcy  law 
porous  wall  boundary  condition  for  (28)  is 

^,(x./f.0)  +  A^,(x.//.ff)  =  O.  (29) 

where  A  is  the  Darcy  porosity  constant.  If  experimental  pressure  data  is  available  on  the 
interface  r  =  // ,  then 

^(jc.H.9)  =  /(x,0)  (30) 

Appropriate  inflow  and  outflow  conditions  as  well  as  an  inner  boundary  conditions  arc 
required  to  complete  the  formulation  for  (28)  for  subsonic  far  field.  The  first  of  these  are 
discussed  in  [7].  Gauss’  theorem  shows  a  uniform  outflow  is  required  as  x  ->  I"  for  the 
solid  wall,  ( A  in  (29))  which  is  linearly  proportional  to  the  body  base  area. 


A.S  in  the  unconfined  case,  the  inner  and  outer  expansions  can  be  matched  to  dominant 
orders.  However,  intermediate  expansions  are  required  for  the  higher  orders.  This  was 
demonstrated  for  the  unconfined  case  in  1 13].  Matching  is  performed  in  the  same  manner 
as  in  Section  2  using  intermediate  variables.  The  appropriate  intermediate  limit  is 


n(S) 


fixed  as  5— >0. 


(31) 


Under  (31), 


r  =• 


-»«»,?=  Snr. 


q(5)  -»0. 


d  ’  ”>  5 

The  inner  expansion  written  in  intermediate  variables  is 


(32) 


U 


K 


,8;x  +-^(logr?r„  -log(5)|+-  (33) 


1 


where  is  the  portion  of  <p^  that  is  regular  at  r’  «  •• .  Terms  3  and  4  arise  from  (23). 

Term  1  in  (33)  is  a  switchback  term  that  is  needed  for  matching  with  the  outer  expansion. 
In  anticipation  of  matching,  thex  derivative  terms  are  assumed  small  and  the  solution  of 
the  remaining  equation 

«0, 

U 

^  sS(.*)log?’  +  g(.r)  +  -”as  7-»0. 

This  implies  that  the  outer  expansion  in  intermediate  variables  is 
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+  5  ^(s(j:)(log  nr^  +  log  5)  +  ^(;c)  +  •••)  +  •■• . 

1'  r  y 

Matching  of  Term  2'  in  (34)  with  Term  l+Term  4  in  (33)  implies  that 
S(x)  =  5.(;c) 

2}t 


(34) 


(35) 


Also,  matching  of  Term  3'  in  (34)  with  Term  (2)  in  (33)  gives 
(p;'^'{r *.0yc  )=  ^(x)  as r’ ->«»  . 

Finally,  Term  I’  in  (34)  matches  with  Term  3  in  (33).  From  this  matching  and  (35),  the 
appropriate  boundary  condition  to  complete  the  formulation  for  (28)  and  (29)  or  (30)  i.s 
4'(x) 


limF^=- 

?-*o  27C 


(36) 


The  pressure  on  the  body  can  be  obtained  from  substitution  of  the  inner  expansion  into 
the  exact  formula  obtained  from  the  energy  invariant 
a'  al 


where  a  is  the  local  speed  of  sound,  U  is  the  frecstream  velocity  and  the  subscript  -  re¬ 
fers  to  freestream  conditions.  From  (37)  and  isentropy,  the  pressure  coefficient,  is 

-  1 


P-P, 

2 


2 

yMl 


1+— iWi  1- 


ilF. 


V' 


n. 

j 


(38) 


where  P  is  the  pressure,  p  is  the  density  and  ^  ^ 
Sub.stitution  of  (23)  into  (38)  gives 


log5+{(p,, 


25^  “  lit 

for  the  pressure  on  the  body,  r*  =  F{0',x)  ,  where  jV,  tp, 


(39) 


.  A  crucial 


step  for  our  final  conclusion  about  wall  interference  comes  from  the  interpretation  of  (23) 
and  the  matching  between  (33)  and  (34)  and  the  decomposition 

<p,  =  g(x)  +  <p,.  (40) 

12 


Term  1  ua  component  determined  only  by  the  solution  of  the  outer  problem  (28),  (29)  nr 
(30)  and  (36).  From  (36),  it  is  influenced  only  by  the  cross  sectional  area  and  not  the 
cross  sectional  contour  of  the  body  at  the  streamwise  station  x  as  well  as  the  wall 
boundary  conditions.  By  contrast.  Term  2  depends  solely  on  the  inner  boundary  value 
problem  (24)  and  is  independent  of  the  wall  boundary  conditions.  These  statements  lead 
directly  to  the  extension  of  our  transonic  area  rule  for  wall  interference  (JARWI) 
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from  //— >oo*  to  H  =0{\),  a  major  impact  of  this  phase  of  the  research.  From  this 
generalization,  more  practical  situations  than  those  for  H  -*=<>  can  be  considered  in 
which  the  model  distance  from  the  walls  is  of  the  order  of  its  length.  These  are  typical  of 
transonic  testing.  It  should  be  noted  that  angle  of  attack  effects  are  higher  order  for  this 
A  =  t9(l)case  as  contrasted  to  A  -4  »  ca.scs  where  they  will  interact  with  the  near  field  in 
the  dominant  orders  through  line  doublct-wall-imaging/rcflection-induced  downwash. 


Fig.  6  Wing-body  tested  in  TsAGIT- 1 28  wind  tunnel. 

7.4  Validations  of  theoretical  and  computational  simulations  for  moderate 
wall  height  case 

Experiments  in  TsAGIs  T-128  wind  tunnel  in  Moscow-.  Russia,  described  in  (9],  have 
been  performed  to  validate  the  theoretical  work  conducted  under  the  contract.  Figure6 
shows  one  of  the  wing-body  configuration.s  tested.  Results  for  pressure.?  over  the  equiva¬ 
lent  body  of  revolution  (EBR)  for  thi.s  wing  body  are  shown  in  Figure  7  that  compares  the 
combined  a.symptotic  and  numerical  method  developed  in  the  contract  with  the  TsAGI 
experiments  for  the  H=0(l)  ca-se  discussed  in  the  previous  .section. 

The  code  is  quite  efficient,  requiring  only  a  minute  of  execution  time  on  a  VAX  3100 
work  station  and  only  100  iteralions  to  obtain  the  2000  iteration  fully  converged  solution. 
Fig.  7  shows  excellent  agreement  between  the  theory  and  experiment.  To  achieve  this 
fidelity,  it  was  important  to  accurately  simulate  the  sting  model  support.  Thi.s  dement  wa.s 
necessary  to  capture  the  proper  rccompression  process  to  ambient  levels.  Additional  vali¬ 
dations  discussed  in  (9]  are  that  the  shock  position  estimates  from  [161  agree  well  with 
the  TsAGI  measurements.  Work  continues  on  specially  designed  experiments  to  adjust 


'  Enunciated  in  [7,8] 
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ihe  level  of  interference  by  altering  the  wall  porosity.  This  will  provide  a  useful  database 
for  comparison  with  the  ^  =0(1)  theory. 


Fig.  7  Comparison  of  present  theory  with  TsAGI  experiment. 


In  1995,  a  larger  aspect  ratio  model  was  tested  in  the  TsAGI  T'128  wind  tunnel  program. 
A  comparison  of  the  larger  aspect  ratio  wing-body  with  smaller  version  is  shown  in  Fig. 
8.  Figures  9  and  10  compare  the  drag  rise  of  the  smaller  and  larger  aspect  ratio  wing  body 
WBl  and  WB2  respectively  with  their  equivalent  bodies  EBl  and  EB2  for  two  different 
wall  porosities.  These  are  expressed  in  terms  of  the  wall  openness  aiea  ratio,  / ,  which  is 
the  area  of  the  wall  perforations  as  a  percentage  of  the  test  section  cross  section  area. 
Values  of  /  were  2  and  10%  for  this  study.  Figures  9  and  10  are,  to  the  Pi’s  knowledge, 
the  first  experimental  confirmation  of  the  transonic  area  rule  for  wall  interference 
(TARWI)  previously  discus.seJ,  Le.,  if 

then 

where  Cp  is  the  drag  coefficient ,  subscripts  WB  and  EBR  denote  the  wing-body  and  its 
equivalent  body  respectively.  Figure  10  shows  the  TARWI  (41)  has  surprising  robustness, 
although  the  aspect  ratio  of  WB2  is  considerably  larger  than  WB2,  the  TARWI  (41) 
still  hold.s  near  =  1 .  This  is  related  to  the  elasticity  of  slender  body  theory  to  not-so- 
slender  shapes  near  sonic  speeds.  It  is  associated  with  the  coefHcient  of  the  x  derivative 
terms  in  the  KG  equation  (28)  being  proportional  to  Af  - 1  as  A/.  -» 1 ,  where 
is  the  local  Mach  number.  Thus,  although  the  cross  flow  gradients  are  no  longer  0(1/6) 
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but  0(1)  for  not-so-slender  shapes  such  as  WB2,  the  x  derivative  terms  are  still  higher 
order.  Accordingly,  the  near  field  remains  harmonic  in  cross  flow  planes  as  in  the  classi¬ 
cal  slender  body  theory.  Other  robu.stncss  of  (41)  should  also  noted.  Although  (41)  is 
applicable  toH  =  0(1) ,  the  nominal  H  is  closer  to  0.1  for  the  Figs.  9  and  10  cases.  This 
is  related  to  the  larger  length  body  associated  with  inclusion  of  the  sting  in  its  definition 
for  the  computational  modeling.  However,  if  a  large  part  of  the  sting  is  at  nearly  ambient 
conditions,  H  =  0(1)  rather  than  the  nominal  H  =  <7(1) . 

As  a  final  comment,  new  effort  in  generalizing  the  experimental  validations  of  the  theory 
to  more  realistic  configurations  has  started  with  1996  wind  tunnel  testing  of  the  Tupolev 
Tu-144  supersonic  transport  in  the  TsAGI  T-128  program.*  One  aspect  of  the  work  will 
be  investigation  of  the  inlets.  The  model  is  modular  in  the  sense  that  it  can  be  te.stcd  with 
inlet  on  and  off  it. 


Fig,  8  Wing-body  configurations  tested. 


*  A  re-engined,  refurbished  version  of  this  aircraft  is  being  prepared  for  special  U.S.-FSU 
flight  test  program  to  start  in  April  1 996. 


(Wing<body)2 


(EBR)2 


Mach  Number 

Fig.  10  Comparison  of  wave  drag  for  wing-body  WB2  and  its  equivalent  body  EBR2. 
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7.S  Non-eireular  wind  tunnel  eectlone 

Many  wind  tunnel  test  sections  are  non-circular.  Typical  U.S.  installations  have  octagonal 
and  rectangular  test  sections.  The  PFs  testing  in  the  T-128  wind  tunnel  strongly  suggests 
that  these  non-circular  cross  sections  have  only  a  mild  influence  on  the  axisymmetric  far 
field  of  a  slender  model  tested  at  transonic  Mach  numbers.  This  observation  motivated 
the  phase  of  the  research  to  be  described. 


If  the  flow  near  the  walls  is  subsonic  which  is  the  case  in  transonic  flow  with  a  subsonic 
freestream,  it  is  reasonable  to  expect  rapid  elliptic  decay  of  the  disturbances  inward  to¬ 
ward  the  centerline  of  the  wind  tunnel.  This  contrasts  to  other  cases  such  as  slightly  su¬ 
personic  frecstieams,  near  choking  and  supersonic  bubbles  of  the  unconfined  flows  pene¬ 
trating  the  walls.  To  explore  this  hypothesis,  the  flow  inside  a  test  section  that  is  a  slight 
perturbation  of  circular  section  will  be  treated.  Figure  1 1  shows  an  example  of  such  as 
perturbation  which  is  an  octagon.  For  generality,  the  following  wall  shape . 

R^h+egm  (42) 

in  which  e«  1  and  his  a  constant.  For  specific  shapes,  it  is  po5.sibIe  to  get  an  approxi¬ 
mate  numerical  order  of  magnitude  for  £  which  can  be  written  as 


Values  of  £  for  square,  hexagonal  and  octagonal  test  sections  are  given  in  Table  I. 


Table  I 


These  are  based  on  the  following  relations  for  an  n  -sided  polygon: 

,  k  =  1,2,3,..., n, 

4  4 

In  the  first  quadrant. 


(44a) 

(44b) 


h  =  O,l,2,3,...,n/4,Oi0^ff/2. 

Using  the  polar  coordinates  previously  introduced  and  referring  to  Fig.  1 1,  as  well  as  the 
outer  limit  (26)  and  expansion  (27)  as  well  dropping  the  subscript  unity  notation  ir^(p^ ,  the 


equation  for  the  perturbation  potential  in  the  outer  region  is 

(if-(r  +  lV>„++f-‘(f<y),+rX=0.  (45) 
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For  convenience,  a  free  jet  boundary  condition,  corresponding  to  A  -> «»  in  (29)  is  con¬ 
sidered.  Accordingly,  the  exact  boundary  condition 
C^(x,Rm  =  Q 

implies 

Since  R  is  independent  of  x , 

=  constant.  (46) 

The  constant  can  be  assumed  to  be  zero  without  loss  of  generality.  CoTre.sponding  to  (42), 
the  perturbation  potential  0  can  be  decomposed  into  the  axially  symmetric  outer  part  cor¬ 
responding  to  £  =0  and  the  secondary  perturbation  associated  with  the  deviation  of  the 
walls  from  a  circular  cross  section.  Thus, 

((l(x,r,0)=i(po(x,r)  +  e(^^(x,r,9)^...  (47) 

Substituting  into  (45)  gives 

which  implies  for  the  rospective  orders 

O(l):{/ir-(y+l)4.}(lb„+7K-),=0  (^8a) 

(9(£):{A-(y  +  l)0o>,„-(y  +  l)^>,A„+7(j^^,,),+-^<>.«,=O  (48b) 

Equation  (48b)  is  the  so-called  variational  transonic  small  disturbance  equation  consid¬ 

ered  in  [7,8]  in  connection  with  other  applications. 


Substituting  (47)  into  the  boundary  conditions  near  the  tunnel  centerline  gives 
Ijinr0,.  =0, 

f-<0 

where  i4(x)  is  the  cro.ss  sectional  area.  (49b) 

On  the  (non-circular)  cylindrical  walls, 

=  ^  +  (50) 

a  Taylor  scries  expansion  for  the  free  jet  condition  (46)  gives 

<p{x,H+£g-¥...,e)  =  (po{x,H-k-eg+...)+<l\{x,H+eg+...,e)+... 

=  ^  (x,  H} + (x,  7/)  +  ^  (x,  if )+...)+.. . 


Equating  like  orders  gives 
(7(l):^!b(x,W)  =  0 
£?(£):^,(x,f/)  =  -g(0)^,(x,f/). 


(51a) 

(5Ib) 
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A  Fourier  decomposition  to  reduce  the  three-dimensional  Transonic  Small  Disturbance 
(TSD)  problem  “PI”. (48b),  (49b)  and  (51b)  for  a  wall  perturbation  from  cylindrical  to 
one  in  two-dimensions  is 

<l>i=fi<l\S^,r)co&n9.  (52) 

ji«0 

This  decomposition  exploits  the  fact  that  the  only  way  that  asymmetry  is  introduced  into 
the  perturbation  problem  is  tlirough  the  multiplicative  factor  g(6)  in  (51b).  Note  also  that 
the  assumption  of  small  perturbations  allows  the  boundary  conditions  to  be  transferred 
from  the  pemirbed  surface  to  the  simpler  cylindrical  test  section’s.  This  is  essential  to  the 
reducticn  of  the  dimensionality  of  the  problem.  Equation  (52)  is  a  factorization  that  re¬ 


duces  the  problem  PI  to  the  form 

{k-(y  +  1)<I)c,}Ku  +  =0  (51  a) 

(5>b) 

4^(x,/f)  =  4,(x.//)A,  (52c) 

A,  =  7  gie) cosnede,  (n  >  0)  (52d) 

A,«»7jJ^(0)cos^0  (52c) 


By  (40)  and  the  italicized  paragraph  following  it,  the  function  g(x)  which  controls  the 
drag  and  pressure  distribution  i.s  the  only  part  of  the  dominant  near  field  that  interacts 
with  the  dominant  outer  solution.  Since  their  is  no  6  dependence  in  this  portion,  the  only 
solution  of  interest  of  (51a)  is  that  corresponding  to  n=0.  Thus  the  effect  of  the  higher 
harmonics  A,  for  «  >  0  arc  negligible  to  this  order.  Effectively,  the  angular  dependence  is 
“averaged  out”.  This  is  another  kind  of  area  rule  for  the  effect  of  slightly  asymmetric 
wall  sections. 

To  quantifythis  effect,  the  mild  transonic  case  corresponding  to  large  K  in  (51a)  was 
considered.  To  simplify  the  analysis,  the  problem  is  reduced  to  a  harmonic 
(incompressible)  one  by  scaling  out  fiT  by  an  x  transformation  and  noting  that  the  second 
and  third  terms  in  (51a)  arc  negligible  in  this  approximation.  The  x  transformation  is 

I  =  x/V^  (53) 

This  is  equivalent  to  another  procedure  that  relates  the  KG  to  the  Prandtl-Glauert  (PG) 
equation  from  the  definition  of  K  in  (26)  and  the  fact  that  the  appropriate  outer  variable 
for  subsonic  flow  is  r rather  than  7  in  the  KG  regime.  This  gives  the  reduced  PG  equa¬ 
tion 

(54) 

where  with  some  redundancy  in  notation  used  in  a  previous  section,  the  transformation 
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X^xlfi 


in  which 


I- Ml 


is  used. 

The  boundary  conditions  (51b)  and  (52c)  are  unaffected  by  the  large  A" approximation. 
These  relations  and  (54)  constitutes  the  problem  PT  which  can  be  solved  by  the  expo¬ 
nential  Fourier  transform  pair 

S  =  (55a) 

The  subsidiary  equation  or  Fourier  transform  of  (54)  is 

1  dY  (n}  ,  V- 

-7?^+“-  +  ^^.=0-  (55a) 

dr  r  dr  ) 

The  transformed  boundary  conditions  are 

limrd.  =0  (56b) 

r-tO  " 

(56c) 


e^dk. 


Equations  (56)  lead  to  the  solution 

A.  =  ff'T7S#l  «7) 

2JCj^I,{kh)  dr 

where  /,  denotes  the  modified  Bessel  function  of  the  second  kind  of  n**  order.  Equation 
(57)  can  be  further  simplified  from  the  solution  for  the  transform  of  a  similar  boundary 
value  problem  to  PT  for  ^g.  This  gives 

^  ^  A,(kl»)/.(fe/i)-t-  K,{kh)I,m\  T' 

'^,.^"2rrl  I^{kh)  J "  2nA/o(iW»)'  ^ 

where  the  Wronskian  relation  is  used  to  get  the  third  member  of  the  equality.  The  differ¬ 
ence  of  the  non-circular  and  circular  cross  section  wind  tunnel  pressures  on  a  body  is 
given  by  the  formula 

where  the  g,  (j  =  0,1)  are  the  <7(1)  terms  of  the  inner  (r  0)  expansion  of  the  outer  so¬ 
lutions  for  and  ^  respectively,  and  5  is  the  maximum  thickness  ratio  of  the  body. 

The  0,  (jr,0)  for  ^  0  do  not  contribute  in  (59)  since 


/o(0)*l. 

/,(0)  =  0,  «>0. 

Finally,  from  the  Convolution  Theorem,  (57)  and  (58),  the  change  in  the  body  pressures 
due  to  wall  asymmetry  AC^,  is 
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To  illustrate  (60),  a  parabolic  arc  body  of  revolution  inside  a  square  cross  section  tunnel 
is  Cu.isidercd.  For  this  case 

r*(jc)  =  &(! -x),Q^  x^  \ , 
where  is  the  body  radius.  Accordingly, 

A  =  TCrl  =iitS^[x^  -2x^  +x*)  (61) 


3tr' 

tan— 

8 


Figure  12  shows  the  effect  of  increasing  Mach  number  on  the  normalized  correction  of 
the  pressures  from  a  circular  to  a  square  test  section,  AC^ ,  where, 

AC 

Ag,*-- (63) 


when  hsl  for  a  parabolic  arc  body,  i.e.,  the  tunnel  average  radius  is  equal  to  the  body 
length.  Note  that  although  the  body  is  in  the  interval  0^  ;t  <  1,  wall  asymmetry  influ¬ 
ences  the  flow  considerably  upstream  of  the  body  nose.  Moreover,  the  largest  effects  ap¬ 
pear  at  the  nose  and  toil  of  the  body  and  the  correction  increases  with  Mach  number  as 
expected.  Another  observation  is  the  rapid  upstream  and  downstream  decay  of  the  effect. 
This  is  consistent  with  the  flow  ellipticity.  Lastly  and  most  important  is  the  smallness  of 
the  effect  which  is  in  sharp  contrast  with  the  results  for  h  =.5  which  shows  a  dramatic 
ten-fold  increase  with  merely  halving  the  wall  height.  Thi.s  effect  is  brought  out  in  Figs. 
13  and  14  for  A/„=0  and  .7  respectively.  In  accord  with  expectations,  Fig.  15  shows  that 
compressibility  increases  the  change  in  pressure  associated  with  wall  asymmetry. 

7.8  Summary,  eonelualona  and  racommandatlona 


Development  under  the  contract  effort  and  related  work  of  a  number  of  unique  and  useful 
tools  to  the  wind  tunnel  test  engineer  has  been  described  in  the  previous  sections.  These 
are; 

1.  Extension  of  the  area  rule  for  blockage  interference  to  wall  heights  of  the  order  of  the 
body  length. 

2.  A  systematic  asymptotic  theory  for  lift  interference. 

3.  Simplified  corrections  for  asymmetric  deviations  of  circular  wind  tunnel  sections. 

4.  Emerging  database  for  transonic  wall  interference  theories  from  Russian  tests.  Items  I 
and  2  apply  to  porous  wind  tunnels.  Item  3  can  be  readily  extended  to  .such  sections. 
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Pressure  Correction  for  Square  Walls,  h=  1 


ICtMBUM 


Fig.  12  Pressure  corrections  from  circular  to  square  test  sections,  parabolic 
body,  A  si. 

Pressure  Correction  for  Square  Walls,  h  •  0.5 


Fig.  13  Pressure  corrections  from  circular  to  square  test  sections,  parabolic 
body,  h  3,3. 
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Pressure  Correction  for  Square  Walls,  =  0 


Fig.  14  Pressure  corrections  from  circular  to  square  test  sections,  parabolic 
body,  effect  of  h  for  =0. 

Pressure  Correction  for  Square  Walls,  M^o  *  0.7 


Fig,  13  Pressure  corrections  from  circular  to  square  test  .sections,  parabolic 
body,  effect  of  h  for  «.7, 
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The  ultimate  impact  of  this  work  is  to 

•  Reduce  computational  intensity  of  transonic  wall  interference  estimation. 

•  Help  optimize  mode!  sizing  to  maximize  test  Reynolds  number  while  minimizing 

wall  interference. 

•  Provide  a  quick  means  of  extrapolating  ground  tests  to  free  flight. 

To  enhance  the  utility  of  these  tools  the  following  further  effort  is  recommended; 

•  Apply  Items  1  and  2  to  corrections  to  drag  polars 

•  Extend  Item  2  to  transonic  flow,  moderate  chords,  and  thickness 

•  Extend  Item  3  to  the  moderate  K  case. 

-  Strong  evidence  exists  that  the  n  sO  solution  of  (51a)  is  appropriate  to  the  strongly 
nonlinear  transonic  case.  The  argument  is  similar  to  that  following  (59). 

-  It  is  envisioned  that  S\x)  "spiky”  behavior  of  the  variational  solution  near  shocks 
will  be  the  principal  modification  of  the  solutions  previously  discussed  for  the 
high  subsonic  large  K  ca.se. 

-  Validation  of  the  subsonic  solutions  against  the  exact  eigenfunction  and  elliptic 
function  Green’s  functions  should  be  performed 
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8.  Significance,  Impacts  of  Research,  Recommendations  and  Future 
Work 

Various  aspects  of  hypersonic  stability  and  transition  as  well  transonic  wall  interference 

problems  have  been  ^dressed  in  the  effort  previously  described.  Major  observations  are: 

1.  Solutions  exhibiting  a  special  waveguide  behavior  not  previously  studied  reveal 
disturbance  amplification  and  transition  mechanisms  in  which  acoustic  disturbances 
reflect  against  the  edge  of  the  boundary  layer  as  if  it  were  a  solid  wail.  Similar 
phenomena  occur  in  the  region  between  the  edge  of  the  shock  and  the  boundary  layer 
sonic  line.  This  work  is  an  outgrowth  of  the  previous  effort  of  the  Principal  and  Co- 
Investigators  concerning  nonparallel  and  shock  ejects  in  hypersonic  boundary  layers 
as  well  as  waveguide  energy  trapping  that  are  not  modeled  by  current  state-of>the« 
art  stability  or  transition  pr^iction  codes.  These  mechanisms  'ompete  with 
conventional  processes  such  as  amplification  of  Tollmein  Schlichting  waves  and 
crossflow  inst^ilitics  as  paths  to  transition.  Our  identification  of  this  physics  should 
be  useful  to  designers  and  engineers.  They  also  will  be  helpful  in  pre  and  post 
processing  of  direct-simulations  and  interpretation  of  such  solutions.  It  is 
recommended  that  the  waveguide  solutions  be  extended  to  three-dimensional  flows 
and  that  unit  experiments  be  performed  to  validate  the  results  of  our  theoretical 
analyses  such  as  those  that  have  been  initiated  by  Maslov  and  his  IPTAM  group  in 
Novosibirsk. 

2.  We  have  developed  an  unsteady  Newtonian  shock  layer  theory  to  treat  aeroelastic 
and  unsteady  aerodynamic  effects  as  well  as  finite  amplitude  hypersonic  stability. 
From  this  model  we  have  discovered  inflectional  instabilities  that  can  be  exploited  to 
enhance  scramjet  combustor  burning  processes. 

3.  The  initial  boundary  value  problem  of  hypersonic  threeslimensionai  unsteady 
stability  has  been  considered  within  an  asymptotic  framework  for  the  first  time 
to  our  knowledge.  A  major  contribution  is  the  response  of  a  point  vibrator.  We 
recommend  that  comparable  experiments  be  performed  using  the  IPTAM  spark 
discharge  technique.  Our  work  contrasts  with  previous  two*dlmensional  analyses 
such  as  (Nciland  1970,  Brown  and  Stewartson  1975)  for  the  viscous  .steady  case  or 
(Malrauth  1993,  Blackaby  el  al.,  1993)  for  unsteady  inviscid  flow.  Perturbations 
within  the  framework  of  an  initial  boundary- value  problem  instead  of  the 
conventional  eigenproblem  provide  important  insight  into  the  space-time  evolution  of 
instabilities  and  receptivity.  A  significant  impact  is  that  it  reveals  other  instability 
paths  besides  Tollmein-Schlichting,  cross-flow,  bypass  types  and  those  mentioned  in 
Item!.  The  main  feature  is  global  interaction  between  shock  and  boundary  layers  as 
well  as  the  role  of  strong  wall  cooling  in  such  interactions.  Additionally,  we  have 
shed  further  light  on  tlw  role  of  three-dimensionality  in  connection  with  upstream 
influence  in  hypersonic  strong  interactions  since  the  (Nciland  1970,  Brown  and 
Stewartson  1973)  studies.  With  thi.s  new  knowledge,  it  is  possible  to  dilineate  regions 
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of  influence  from  three-dimensional  perturbation  sources.  For  steady  and  low 
frequency  fluctuations,  the  most  intensive  influence  appears  in  the  downstream  wake 
within  the  downstream  Mach  cone.  At  high  frequencies,  effects  are  concentrated  on 
the  source  axis. 

This  theory  gives  improved  basic  understanding  of  the  interaction  of  separation, 
stability  and  transition.  The  practical  implication  for  hypersonic  vehicle  design  is  t}v>( 
shock  interference  heating  induced  by  different  devices  within  the  boundary  layer  and 
strong  separation  can  be  reduced  by  techniques  suggested  by  our  triple-deck  modeh. 
Moreover,  instabilities  and  other  unsteady  effects  can  be  exploited  to  enhance  mixing 
for  scramjet  propulsion  .such  as  with  vibrating  elements  with  tailored  or  maximized 
receptivity.  In  addition,  our  analyses  could  with  further  development  provide 
engineers  and  designers  manufacturing  and  flow  quality  tolerances  for  surface 
roughness,  such  as  near  leading  edges. 

We  recomnend  that  increased  understanding  of  disturbance  propagation  from  this 
pha.se  of  the  work  be  used  to  improve  computational  schemes.  They  could  suggest 
stable  and  fast  running  algorithms  for  hypersonic  interactive  flows  with  strong 
cllipticity  such  as  shock-boundary  layer  interaaions,  and  even  separated  “self- 
interaction”  flows  behind  a  body  flap  or  spoiler.  Other  recommended  generalizations 
that  appear  within  our  grasp  are  more  complex  flows  involving  chemical  reactions 
and  radiation. 

5.  Our  asymptotic  formulations  give  basic  similitudes  that  we  have  applied  to  collapse 
hypersonic  cone  transition  data  to  universal  curves  which  is  the  basis  of  fast 
preliminary  design-oriented  aerodynamic  and  heat  transfer  prediction  methods.  We 
recommend  that  our  procedure  be  extended  to  handle  bluntncss.  This  is  a  major  gap  in 
the  design  of  reusable  launch  vehicles  such  as  the  X-33. 

6.  From  a  new  stability  analysis  with  application  of  a  Darcy  law  boundary  condition,  we 
have  obtained  strong  evidence  that  absorbing  walls  such  a.s  those  associated  with 
thermal  protection  systems  (TPS)  can  quench  second-mode  inviscid  instabilities 
within  hypersonic  boundary  laycn.  V'e  recommend  that  the  approach  be  generalized 
to  viscous  disturbances.  Ground  and  flight  experimental  validations  of  the  benefits 
should  be  performed  so  that  a  practical  system  can  be  developed. 

7.  At  transonic  speeds,  our  transonic  area  rule  for  wall  interference  (TARWT)  was 
extended  from  cases  where  the  wall  height  parameter  in  units  of  the  model  length 
A  ■  o(l)  to  A  =  0(1) .  We  have  performed  experiments  in  Russia  in  the  TsAGI  T- 1 28 
wind  tunnel  that  validate  this  rule.  It  reduces  the  problem  for  transonic  blockage 
interference  from  three  to  a  two-dimensions.  Wc  have  also  applied  matched 
asymptotics  to  provide  a  coasistent  approximation  scheme  to  estimate  porous  wail  lift 
interference  correctioas.  This  theory  has  been  refined  in  our  effort  to  account  for  the 
influence  of  concentrated  wing-body  Juncture  vortices  that  give  nonlincaritics  in  the 
variation  of  lift  with  incidence.  A  recommended  reflnement  is  to  generalize  the 
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approximations  in  the  same  fashion  as  those  for  blockage,  i.e.,  the  lift  interference 
theory  assumed  that  h  =  oO).  Extension  to  A  =  0(1)  will  give  a  camber  effect 
associated  with  the  variable  tunnel  induced  upwash  field  along  the  model  length.  An 
immediate  output  of  our  wall  interference  model  is  that  it  can  be  used  to  extrapolate 
polars  from  the  wind  tunnel  to  free  flight  or  another  facility.  The  TARWI  allows  this 
estimation  to  be  made  without  the  need  for  a  large  scale  CFD  process  involving 
significant  gridding  pre-processing.  Instead,  only  a  desktop  computation  is  required  to 
obtain  the  nonlinear  blockage  effect.  This  will  help  optimize  model  sizing  to 
maximize  test  Reynolds  number  while  minimizing  wall  interference.  Further 
validations  and  extensions  of  the  theory  are  planned  such  as  our  imminent  tests  of  a 
4%  Tupolev  supersonic  transport  Tu-144  model  in  the  TsAGI  T*128  wind  tunnel  in 
the  near  future.  A  three  view  of  this  configuration  and  the  location  of  the  pressure  taps 
are  shown  in  Figs.  I  and  2. 


Fig.  1  Tu-144  three-view. 
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